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COVERAGE OF JOURNAL 


THE aim of this journal of abstracts is to give complete coverage of papers in the field of statistical 
theory and new contributions to statistical method. Papers which report only applications or examples 
of existing statistical theory and method will not be included. There are approximately two hundred 
and fifty journals published in various parts of the world which are wholly or partly devoted to 
the field of statistical theory and method and which will be brought within the scope of this journal of 
abstracts. A complete list of journals covered is printed in the annual Index Supplement. In the case of the 
following journals, however, being those which are wholly devoted to statistical theory—all contributions, 
whether a paper, note or miscellanea, will be abstracted : 

Annals of Mathematical Statistics 

Biometrika 

Journal, Royal Statistical Society (Series B) 

Bulletin of Mathematical Statistics 

Annals, Institute of Statistical Mathematics 
Within the larger group of journals, which are not wholly devoted to statistical theory and method, there 
are some journals which have the vast majority of their contributions in this field. These journals, there- 
fore, will be abstracted on a virtually complete basis : 

Biometrics 

Metrika 

Metron 

Review, International Statistical Institute 

Technometrics 

Sankhya 
After experience in the publication of this journal of abstracts it may be found desirable to add further 
journals to these lists. In any case readers may be assured that all papers properly to be included in this 
journal of abstracts will be included irrespective of such notification. If any reader of this journal discovers 
a paper which happens to have been overlooked, the General Editor will be pleased to be informed so 
_that the appropriate abstract can be made: always provided that the date of publication is after 1st 
October 1958, when the abstracting for this journal commenced. 

In addition to the ordinary journals, there are two kinds of publication which fall within the scope 
of this journal of abstracts. They are the experiment and other research station reports—which occur 
particularly in the North American region—and the reports of conferences, symposia and seminars. 
Whilst these latter may be included in the book review sections of journals it is unusual for any individual 
contribution to be noted at any length. ‘These publications are, in effect, special collections of papers 
and for this reason the appropriate arrangements will be made for them to be included in this journal. 
By the same token, abstracts of papers given at conferences and reproduced in an appropriate journal will 
- be disregarded until the definitive publication is available. 


FORMAT OF JOURNAL 

The abstracts will be up to 400 words long—the recommendation of UNESCO for the “ long ” 
abstract service: they will be in the English language although the language of the original paper will 
be indicated on the abstract together with the name of the abstractor. In addition, the address of the 
author(s) will be given in sufficient detail to facilitate contact in order to obtain further detail or request 


an offprint. Suitable indexes will be provided annually in a supplement. 
iii 


The scheme of classification has been developed upon lines that will facilitate the transfer to punched 
cards of the code numbers allocated to each abstract: to allow for future development it is suggested that use 
is made of 4-column fields. Each abstract will have two classification numbers: the primary number in 
heavy type to indicate the basic topic of the paper and the secondary number in brackets to take account of 
the most important cross-reference. The sheets of the journal are colour-coded according to the twelve main 
sections of the classification and it should be noted that it is the main section number of the primary 
classification which determines the colour code for each abstract. It is believed that this method of 
colour-coding the pages will provide a distinctive visual aid to the identification of abstracts both when 
the journal is in bound form or dismantled and filed on cards or in binders. The format and simple 
binding allows of the following alternative treatments by users of the journal : 

(a) Leave intact as a shelf-periodical. 

(6) Split and filed in page form according to the main sections of the classification. 

(c) Split and guillotine (single cut) each page ready for: 

(i) pasting onto standard index cards. 
(ii) filing in loose-leaf or other binders—for which the appropriate holes are punched. 


It will also be possible for those users who need a completely referenced file to insert skeleton cards or 
sheets according to the secondary classification number provided on each abstract. 


et 


oes saa 


Q Primary and 
' secondary 
’ classification 
| 


Holes for filing 


into binders 

| 

| 

| 

L 
O Serial number 

Mi Ho abstract 
Read 
a ea ee ere Guillotine guide marks — — — — — — | 
en 

0 Primary and 

| secondary 

; classification 
Holes for filing 

into binders 

| 

| 

i) 

v 

© Serial number 


UA of abstract 


nase! 


JOURNAL OF THE 
ROYAL STATISTICAL SOCIETY 


Series A (GENERAL), four issues a year, 25s. per part, annual subscription £4. 4s. (U.S. $12) post free; SERIES B 
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Volume 1. CORRECTIONS 


1/544 1.5 Frankx Author’s name should read “* Franckx ”’. 

1/566 2.9 Bahadur Third line from end: last word should read “ generalisation ”. 
1/588 3.3 Ruben Name in 6th line should read “‘ Behrens-Fisher ”’. 

1/680 8.0 Mikulski & others Title of paper should read “. . . mean properties . . .”. 


Line 5 should read w. 
Line 15, insert word “if” before w(r). 


1/683 8.0 Oderfeld Character at end of Ist line should read w. 

1/749 11.1 Federigh Author’s name should read “‘ Federighi ”’. 

Author Index No. 4 and in Amend author’s name to “ Federighi ”’. 
Supplement 


Volume 2. CORRECTIONS 


2/121 8.9 Barnard The secondary classification should read (10.2). 
2/131 8.8 Lieberman & Bowker The expression in col. 2, line 5 should read 1—f. 
2/213 10.1 Steinhaus & Urbanik Column 2, line 11 should read Se ol BS 


2/219 11.1 Keuls The journal abstracted is Statist Neerlandica. 
Author Index No. 1 Add: Barnard, G. A. 8.9, 8.9 

Durbin, J. 10.2 

Kiefer, J. 9.0, 9.0. 


EpitoriAL Note: In this and future issues all serial numbers of abstracts, whether at the bottom left-hand corner 
of an item, quoted as a reference or mentioned in a Corrections slip, will show the volume number of this 
Journal in bold face: e.g., 1/462. 


ARNAIZ, V. G. & LASUEN, J. R. (Inst. Nacional de Estadistica, Madrid) 0.6 (—.-) 
Economic significance of the coefficients in input-output analysis—Jn Spanish 
Estadistica Espaiiola (1959) 5, 27-33. 


The economic significance of all the inversion process in 
Leontief’s input-output matrix is explained from a 
didactical point of view. The authors emphasise the 
importance of the matrix J— A, not only as a necessary 
intermediate step between A and (I— A)~', but also as 
a meaningful element in itself; whose interpretation is 
not so obvious because of its two kinds of coefficients. 
A deeper comprehension of each algorithm is urged as 
an instrument for research and its applications. 


(F. Azorin) 


2/225 


BARTON, D. E. & MERRINGTON, Maxine (University College, London) 0.1 (2.5) 
Tables for the solution of the exponential equation—IJn English 
Biometrika (1960) 47, 439-446 (10 references, 2 tables) 


The authors consider the solution of the exponential 
equation exp(—a)+ka = 1 for a, given k. For values 
of k = 0-050(1)...1-000 the value of a is given to seven 
decimal places. For k<0-050 it is enough to assume 
a= 1/k. 

Two statistical contexts in which the solution of this 
equation intervenes are discussed. The first is the 
equation for the maximum likelihood estimator of the 
Poisson parameter, 4, when the observed distribution 
is truncated and the zero group is missing. The second 
is the procedure suggested by Laplace for approximat- 
ing to the difference quotients of zero. Laplace’s 
transformation is compared with those put forward by 
De Moivre (1756), Binet& Szekeres [Ann. Math. Statist. 
(1957) 28, 494-498], Arfwedson [Skand. Aktuartidskr. 
(1951) 34, 121-132] and Hsu [Ann. Math. Statist. (1948) 
19, 273-277]. It is shown that Laplace’s approximation 
is very good indeed and superior over the whole range 
to those recently proposed. 


(Florence N. David) 
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CHIASSINO, G. (University of Bari, Italy) 


0.2 (2.2) 


On the analytic representation of survival curves—In Italian 
Riy. Ital. Econ. Demogr. Statist. (1959) 13, 657-660 (6 references, 1 table) 


The author reconsiders a form first suggested by Perrone 
for the smoothing of the survival curve in a life table 
by which the validity of the analytic graduation is 
extended to ranges larger than those usually obtaining 
with other functions. 

The suggested solution, according to Perrone, is based 
on a convenient transformation of the age variable: the 
modification introduced by the author is shown to 
produce better results with regard to the experience of 
the recent 1950-53 Italian table calculated for the male 
sector of the population. 


(V. Levis) 


loa ke 


HARTLEY, H. O. & LOFTSGARD, L. D. 


(lowa State Univ. and North Dakota State Coll.) 


0.8 (0.7) 


Linear programming with variable restraints—In English 
Iowa State Coll. J. Sci. (1958) 33, 161-172 (2 references, 1 table) 


Linear programming methods are becoming increasingly 
important for analysing economic situations. In 
particular, programming methods have been extensively 
applied to farm management problems, where maximum 
revenue defines optimum solutions. In these problems, 
the magnitude of attainable revenue depends on the 
alternative uses of resources and boundaries of the 
planning situation as represented by available resource 
supplies. This paper deals with the latter stipulation 
concerning resource supplies. ‘ 

A more comprehensive method of programming is 
one that allows variation of resource supplies. That is, 
a method that determines continuous optimum solutions 
when one resource supply is varied within a relevant 
range and other resource supplies are held constant. 
One such method developed by Candler is a modified 
simplex solution for linear programming with variable 
capital restrictions. 
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This paper discusses a different modification of the 
simplex solution with continuous variation of one 
restraint. It is believed that the method proposed here 
has considerable advantages for programming on high- 
speed computers. The procedure is numerically illus- 
trated by use of a sample problem from Candler’s article. 


(H. O. Hartley) 
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MIHOC, G. (University “ C. I. Parhon ”, Bucharest) 0.8 (2.3) 
Some remarks in connection with the application of linear programming—IJn Rumanian 
Revista de Statistica (1959) 12, 13-17 


In this paper the author suggests a method for solving 
the problem of linear programming under the hypothesis 
that the coefficients of the linear form to be extremised 
are normally distributed random variables with given 
means and variances. 

The result obtained shows the statistical character of 
the linear programming method when applied to practical 


problems. 
(M. Iosifescu) 
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NORDBOTITEN, S. (Statistisk Sentralbyra, Oslo, Norway) 0.8 (11.5) 


Linear programming and automatic computers—In English 
Skand. Aktuartidskr. (1959) 42, 61-72 (9 references, 1 chart) 


Two applications of linear programming methods are 
indicated. The first problem treated is to find a plan 
of production which maximises the excess of exports 
over imports. The model employed is of the input- 
output type. The second problem concerns the alloca- 
tion of sampling units in sampling a population of 
manufacturing establishments. The population is sub- 
divided into industry groups. A second classification 
is into districts. Estimates are needed for the total 
production of each industry group and the total pro- 
duction within each district. It is assumed that the 
cost of sampling n,; establishments in the ith group 
within the jth district is of the form a;,-c;,/n;;. The 
minimisation of the cost with restraints placed on the 
allowable variances of the estimated totals then becomes 
a case of linear programming. 

The simplex procedure for solving a linear program- 
ming problem is reviewed. A detailed flow-chart is 
presented for applying the procedure on an automatic 
computer. 


(B. Matérn) 
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SHEPHARD, R. W. (University of California, Berkeley) 


0.8 (-.-) 


Non-linear programming. I. Linear representation of separable concave functions 


in production planning—In German 


Qualitatskontr. Operat. Res. (1960) 5, 131-136 (3 tables, 3 figures) 


The method of non-linear programming as treated in 
this paper is a specific form of concave methods: it is 
supposed that the objective function (the function to be 
maximised or to be minimised) consists of a sum of 
separable concave functions of only one variable each, 
and that the constraints of the problem are linear. It 
is shown that problems of this type can be reduced to 
linear programmes. 

The following example illustrates the above-mentioned 
method. A firm manufactures and sells three different 
commodities produced by only one machine with a 
limited capacity. This results in a linear inequality. 
The profit per unit of each commodity, being constant 
in certain intervals of production, decreases from interval 
to intervai (price differentiation). The number of 
intervals is 2 for commodity A and 3 for Band C. After 

exceeding a fixed limit of production the profit becomes 
zero, consequently the profit function of each commodity 
is a concave function composed of sections linear within 
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SHEPHARD, R. W. 


in production planning—J/n German 


(University of California, Berkeley) 
Non-linear programming. I. Linear representation of separable concave functions 


each interval, the sum of all commodity profit functions 
gives the objective function (total profit function). In 
order to reduce this programme to a linear form, for 
each possible unit profit of each commodity, non- 
negative variables must be introduced being permitted 
to become at most as large as the corresponding produc- 
tion interval. By this means the number of variables 
is increased, but the programme has become linear; as 
only one capacity constraint is assumed, it is solvable 
in an elementary manner without using the simplex 
method. The author points out that the approximation 
to a concave function with continuously decreasing 
derivation by means of a polygon leads to similar 
programmes. 

Hereafter, the above example is extended by increasing 
the number of machines producing the three com- 
modities. The capacity of the four new machines is 
restricted and so there are four linear inequalities as 
capacity restrictions. A matrix (a;;) indicates the 


continued 
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continued 


Qualitatskontr. Operat. Res. (1960) 5, 131-136 (3 tables, 3 figures) 


number of hours used on the ith machine to produce 
one unit of the jth commodity. The procedure of 
making the programme linear is the same as that cited 
above. 

Finally, the author presents some ideas on the applica- 
tions in practice, showing how to remove bottle-necks 
by marginal interpretation. He finishes by calling 
attention to the primary-dual relation and illustrates 
the notion of the shadow prices (opportunity costs) to 
demonstrate how to determine the marginal productivity 
by means of the duality theorem. 


(W. Dinkelbach) 
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SHEPHARD, R. W. 


(Universityfof California, Berkeley) 


0.8 (-.-) 


Non-linear programming. II. The Kuhn-Tucker theorem, quadratic programming and 


the gradient method—Jn German 


Qualitatskontr. Operat. Res. (1960) 5, 144-149 (4 references) 


In part one of this paper [Qualitdtskontr. Operat. Res. 
(1960) 5, 131-136: abstracted in this present journal 
No. 2/231, 0.8] the author is concerned with the linear 
representation of separable concave functions in the 
field of production planning. 

In part two of this paper the Kuhn-Tucker theorem 
is stated and discussed; this includes its application to 
quadratic programming (the algorithm of Wolfe). This 
is followed by an outline of the gradient method. This 
paper contains a bibliography of the original papers on 
those subjects. 


(W. Dinkelbach) 
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STAHLKNECHT, P. 


(Karl-Marx-Universitat, Leipzig) 


0.8 (-.-) 


Linear programming applied to problems of managerial science with changing conditions—Jn German 
Wissensch. Zeit. Karl-Marx-Univ. Leipzig (1958/59) 8, 693-702 


The present article is part of the author’s thesis at the 
Faculty of Science, Karl-Marx-University of Leipzig. 
The author briefly discusses the problems underlying 
linear programming. In the following paper he quotes 
the methods most frequently used for finding optimal 
solutions of linear programmes and describes the 
particular aspects to which he devotes this paper. 

In management science, linear programmes often occur 
for which the bounds of some variables are not constant. 
One is not therefore mainly interested in the solution 
of one fixed linear programme in the usual sense of the 
word, but particularly solutions which show how the 
result is influenced by changes of the bounds. 

In developing a method of constructing such solutions 
the author starts with a special problem; this is treated 
numerically in extenso. The experiences gained hereby 
can be generalised without difficulties, leading to a certain 
procedure applicable in general cases. It consists mainly 
in the repeated performance of some simple mathematical 
operations. The following problems may be treated in 
the manner described by the author: 


p= par ys DXi = Min 
i=1 
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under the conditions 
n 
Meese dD, My 20 kantlaem 
i=1 


OSxSr 


The quantities po, Pi, Vz. V., G= 1, ...2; k =n+l,..., 
m) are given constants, while information concerning 
the bounds r; is not available in advance. 

As an illustration, the well-known blending problem 
is used, in the case where the optimal composition of a 
mixture of m gases is to be found, and whose calorific 
value has to remain in given bounds while the pro- 
portion of sulphur in the mixture must be below a certain 
limit. In addition it is required to produce a certain 
quantity of the mixture. The solution especially has to 
show the influence of certain sorts of gases on the 
optimum composition. 

In a final section practical experiences are reported. 
Only the hand computation is treated; the use of the 
electronic computers is not considered. 
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AUSTIN, T. L., FAGEN, R. E., PENNEY, W. F., & RIORDAN, J. 


1.3 (11.3) 


(Silver Spring, Md.; Hughes Aircraft Co.; Silver Spring, Md.; and Bell Telephone Labs.) 
The number of components in random linear graphs—IJn English 
Ann. Math. Statist. (1959) 30, 747-754 (3 references, 1 table) 


Given n distinct (labelled) points, a random linear graph 
is given by m selections of pairs of points made in- 
dependently and at random, each of the (;) pairs having 
the same chance of selection at each trial, upon joining 
each selected pair by a (labelled) line. A ‘‘ component ”’, 
or “ piece ’’, of such a graph is a collection of k of the 
n points, such that any two of the k points are connected 
by some path (possibly through other of the k points), 
and such that no point other than these k points is 
connected to any of the k points. 

Let Timp denote the number of such graphs having 
n points, m lines, and p parts; let C,m = Tm; denote the 
number of connected graphs having n points and m 
ines: Let 


nmp 


SVCD ahs Nis paps B ns z’, and 
n! m! 


C(x, y) = EEC, ~— 2. 
n! m! 
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BARTON, D. E. & DAVID, Florence N. 
A collector’s problem—In English 


(University College, London) 


The authors show that T(x, y, z) = exp[zC(x, y)], and 
exhibit specific formulae for T,,,) and C,, resulting 
from the equality just given. The results described are 
unwieldy for large nm. Consequently, the authors 
determine by other means formulae for T,, ,- 11 Ge et 
and for T;, », 1 = Cn, », and quote the explicit forms for 
ene tle eNO? 

The average number of “‘ components”? is written 
down in explicit but elaborate form. An alternative 
development is presented allowing certain asymptotic 
approximations to be used. Finally, the approximation 
is shown, compared with actual values, for M,,,, 1 = 3, 
pate KO) 


Note: In expression (13), the factor (”) is missing 


inside the summation. To obtain expression (15) from 
(14), multiplication should be by x"/n!, not x"*1/(n+1)!. 


(J. A. Lechner) 


1.3 (4.7) 


Trab. Estadist. (1959) 10, 75-88 (3 references, 2 tables) 


This paper deals with two aspects of the classical balls- 
in-boxes occupancy problem; first the case when the 
number of balls to be dropped is fixed and the number 
of empty boxes is a random variable, and second the 
case in which the number of empty boxes is fixed and 
the number of balls to be dropped is the variable. 

N identical compartments are supposed and one 
further control box; n balls in all are dropped indepen- 
dently. The probability that n—j balls will fall in the 
control box and j in the N compartments to give k 
empty compartments and N—k occupied is found: the 
probability is also calculated of k empty cells and at 
least n—b balls in the control box. The factorial 
moments of k, conditional on n—j balls in the control 
box, and the conditional factorial moments of k when 
there must be at least n—b balls in the control box is 
found. 

The obverse of the problem, when n becomes the 
random variable, is considered under the sequential 
aspect; it is supposed that balls are dropped one by one 
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independently and at random, and that this process 
continues until there are at least (V—k) occupied boxes 
and b in the control box: or, there are precisely N—k 
occupied boxes and at least b in the control box. The 
probability density function, the moments of n and 
some asymptotic distributions are studied. 

Finally, the exact evaluation of some upper percentage 
points are given and the error involved when these points 
are calculated by means of the asymptotic distributions 
previously mentioned is stated. 


(J. M. Garcia) 


BOROVKOFF, A. A. (Moscow) 


1.5 (10.0) 


Limit theorems on the distributions of maximum of sums of bounded latticed random 


variables. I.—Jn Russian 


Teor. Veroyat. Primen. (1960) 5, 137-171 (22 references) 


tie bounded latticed identically distributed random 
variables €,, ¢,, ..., are considered by the author in 
this paper. The local and integral theorems for the 
first passage time 7, over the barrier x>0 in the random 
wanderings along the straight line with the quantity of 
jump ¢, are studied. The formulas for Pr(y, = 7) and 
Pr(y,.>”) in obvious form are obtained for the full 
*“ spectrum ”’ of values x, beginning with x = o(n) until 
x, is equivalent to the product maximum jump €, by n. 

The theorems for Pr(7,.>) simultaneously are integral 
theorems for the maximum of sums ys Gives 2s 

ksy 


...) 1). The formulas for first moments y, and the 
distribution of the quantity of the first excess over the 
barrier x are also obtained. 


(A. A. Borovkoff) 
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ERDOS, P. & RENYI, A. 
On random graphs. I.—Jn English 


(Math. Inst., Hungarian Academy of Sciences, Budapest) 


1.3 (1.5) 


Publ. Math., Debrecen (1959) 6, 290-297 (8 references) 


Let us consider all possible graphs with n given vertices 
Paros ss P and IN edges. Let us select one of them 
at random so that each of these graphs should have 
the same chance to be selected. Such a random graph 
should be denoted byT,, y. Put V.(n) = [jn logn+cn]. 
Denoting by P,[n, N.(n)] (k = 0, 1, ...) the probability 
that the greatest connected component of T,, N.(n) 
consists of n—k points, the authors show that for 
n—>-+oo this distribution tends to Poisson’s law with 
mean value exp(—2c). For k = 0 this theorem implies 
that the probability that [,, N.() should be connected 
tends for n> + 0 to exp(—e °°). This includes previous 
(unpublished) result of Erdés and Whitney. It is shown 
further that in the above case when T,, N(m) consists 
of a connected component having n—k points, then 
the remaining k points are all isolated with probability 
tending to 1. The authors further determine the limiting 
distribution of the random variable Vn defined in the 
following way: the edges of a graph with n possible 
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vertices are chosen successively so that after each step 
the next edge is chosen randomly from the possible edges 
not yet chosen, until the graph becomes completely 
connected. Vn denotes the number of edges of the 
resulting graph. 


(I. Vincze) 
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HALDANE, J. B. S. (Indian Statist. Inst., Calcutta) 
The addition of random vectors—In English 
Sankhya (1960) 22, 213-220 (3 references) 


The author considers the addition of random vectors 
x, such that x has uniform distribution over spheres 
| x || = for each c. If x and y are two such in- 
dependent vectors then it can be seen that z= x+y 
is also uniformly distributed over spheres. The first 
five cumulants of ||z|| are derived in terms of the 
cumulants of ||x || and ||y |. The author defines 
what are called vectorial cumulants 0, of || x || which 
are functions of the cumulants x, of | x |], with the 
following properties: 


., K, only; 

(ii) x, is a function of 0,, ..., 0, only; 

(iii) 6, x+y ||) = 8. x [)+4,([ » I). 

Explicit expressions for 0, in terms of x, and vice-versa 
have been given. The vectorial cumulants 0, are useful 
in studying the distribution of the resultant sum of any 
number of random vectors x, y, Z, .... First, 0, of the 
resultant is obtained by using the addivity property and 
from them the corresponding x, are computed. 


(i) 0, is a function of K,, .. 
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HODGES, J. L., Jr. & LE CAM, L. 


(Univ. of California, Berkeley) 


1.0 (2.1) 


It has been noted for the gamma distribution that 
0, = 0 for r=2. The author mentions possible applica- 
tions of the distribution of the resultant sum of random 
vectors in some problems in Physics and the study of 
migration in Demography. 


(S. R. S. Varadhan) 


1.5 (2.5) 


The Poisson approximation to the Poisson binomial distribution—In English 
Ann. Math. Statist. (1960) 31, 737-740 (5 references) 


Let X, be a random variable which can be realised as 
either zero or one, and let Pr{X, = 1} = p;, i = 1, 2, 


ees 


n 
n. The random variable S = )) X; is sometimes 
1 


said to have a distribution called “‘ Poisson binomial ”’. 
It is known that S tends to become a Poisson variate 
in the limit (on 2) when Lp; = 4 is fixed and « = max{p;} 
tends to zero. It can also be shown that the absolute 
difference between the distribution of S and that of a 
Poisson variate with parameter Lp; is bounded by 
Ca'/>, where C is independent of n and the p;. This 
paper presents a simple proof that the bound Cx’/° just 
stated may be improved to 3". 

The proof is accomplished in three stages: the authors 
are concerned with the pedagogical as well as the mathe- 
matical aspects of their work. In this connection, their 
theorems are stated as approximation theorems rather 
than limit theorems, with the corresponding limit 
theorems indicated briefly at the conclusion. 


(D. H. Shaffer) 
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MANDL, P. (Math. Inst., Czechoslovak Academy of Sciences, Prague) 


1.5 (10.1) 


On limit properties of probabilities in groups of states of an homogeneous 


Markoff chain—Jn Russian 


Cas. pést. math. (1959) 84, 140-149 (3 references) 


Let us consider a finite homogeneous Markoff chain. 
If T is a class of states of this chain and w eT, then the 
symbol Pr(x, = w| 7”) denotes the probability that the 
chain is in state w at the moment n provided that it 
has been in states of the set T at all moments between 0 
and n inclusive. The author studies limit properties of 
these probabilities for noo. 

The class T is called regular if the matrix formed by 
probabilities of mutual transitions between states of T 
is non-decomposable and aperiodic. It is proved in 
the paper that if T is regular, the corresponding limits 
lim Pr(x, = w| 7”) always exist (for every weT), are 


n~?oo 
positive and do not depend on the initial probability 
distribution. 

Further, the author investigates the case when the 
class is divided into sub-classes T,, ..., TJ; such that, if 
for any two states w,, w, € T the transition probabilities 
from w, to w, and from w, to w, are both positive: 


then, w, and w, belong to the same sub-class. The 
system of these sub-classes is partially ordered by the 
relation 7; <T, if the transition from 7; to T, has positive 
probability; the corresponding transition path being 
composed of states belonging to T only. Let us denote 
by p; the characteristic number of the matrix corres- 
ponding to 7;, with maximum absolute value. Let 


p = max p;, and let k be the greatest number such 
1s<jsl 


that there exist k sub-classes T;, <...<T7,, with character- 
istic numbers p;, = p. 

If all classes T; are regular and the chain begins its 
evolution in a state belonging to a class which is minimal 
with respect to the partial ordering <, then for every 
weéT there exists the limit lim Pr(x, = w| 7”). It is 


positive if, and only if, w belongs to a sub-class T; such 
that there exist k sub-classes T;,<T;,<...<T;,ST; 


Jk = 
with p;, = pi = 1, 2, ..., k). If there exists only one 
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MANDL, P. (Math. Inst., Czechoslovak Academy of Sciences, Prague) 1.5 (10.1) 
On limit properties of probabilities in groups of states of an homogeneous continued 


Markoff chain—Jn Russian 


Cas. pést. math. (1959) 84, 140-149 (3 references) 


such sequence 7,,<T;,<...<7;, with p;, = p(i= 1, 2, 
..., k) then the limits do not depend on the initial 
probability distribution provided that the chain starts 
in a minimal class. 

In these cases the probability Pr(7”) that the chain 
remains in states of T at any moment between 0 and n 
is of such an order that 

lim [Pr(T")/n*~ *p"]>0. 


n~>o 


(F. Zitek) 
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MASLOFF, K. A. (Kharkoff University) 


On certain limit theorems in probability theory—Jn Russian 
Teor. Veroyat. Primen (1960) 5, 54-83 (10 references) 


Let Gn (1 SkSk,; fel 2s abe a. sequence, of a: 
series of independent random variables, ¢(x, y) is any 
function and the random variables (,, , are determined as 


Cs, fe ca, 1> Grains = Pn, k3 Snes k= 152) 2.2K, —1. 


The author finds a sufficient condition for the existence 


of a limit distribution ¢,, ,,(2— 00) and for the form of 
this distribution. 


(K. A. Masloff) 
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NARAYANA, T. V. (Institut H. Poincaré and McGill Univ., Quebec) 
A partial order and its applications to probability theory—Jn English 


Sankhya (1959) 21, 91-98 (4 references) 


An r-partition of n is a set of integers (¢,, ..., ¢,) with 
t,+...+t, =n. The partial order in question is the 
following: an r-partition (¢,, ..., t,) dominates (f1, ..., 
h) if t,+...+¢,20,+...+¢(7 = 1, 2, ..., r—1). The 
author shows, using a geometric argument, that the 
total number of ‘‘ dominations ’’ among r-partitions is 


(7)(-2,) 


One application of this partial order is the following 
problem in probability theory: “it is known that two 
candidates A, B for election will secure m and n votes 
(m>n) respectively, what is the probability that through- 
out the scrutiny A holds a L-lead over B(1 <L<m-—n)?” 
The author also discusses the connection with the 
following game: there are two coins | and 2 with p,, 
p> as the probability of heads occurring when the coins 
are tossed. The first toss is made with coin 1 and each 
succeeding trial is made with coin 1 or 2 according to 
whether the result of the preceding toss was tails or 
heads. 


(R. Ranga Rao) 
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PROKHOVNIK, S. J. (University of New South Wales) 


Pattern variants on a square field—IJn English 
Psychometrika (1959) 24, 329-341 , 


A quantitative approach to the psychology of pattern 
recognition requires knowledge of the number of possible 
variants of any particular pattern. The general solution 
for the number (p/m7) of pattern variants that p counters 
can form on a square network of m? positions obtained 
by elementary group theory. The exact solution is 
given in terms of the different types (symmetric, asym- 
metric, etc.) of patterns possible in an approximate 
formula for the total number of patterns is also developed 

Different considerations apply for odd and even cases, 
both for networks and for number of points forming the 
pattern. Further, symmetrical patterns of p points over 
an m? network are possible only if p is a multiple of 
four and, in the case when m is odd, also if p—1 is a 
multiple of four. 

A table is given showing a selection of numerical 
evaluations of Gg, G,, G,, and G,, and the total number 
of patterns for networks up to 67 positions. It is clear 
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RANGA RAO, R. & VARADARAJAN, V. S. (Indian Statist. Inst., Calcutta) 


A limit theorem for densities—In English 
Sankhya (1960) 22, 261-266 (8 references) 


In this paper the problem of almost everywhere con- 
vergence of the density of the normalised sum of a 
sequence of independent and identically distributed 
random variables to the density of the standard normal 
distribution is studied. In particular the following 
results are proved. 

cine c,, ... be a sequence of independently and 
identically distributed random variables with mean 
zero and variance unity; F,(x) the distribution function 
of the normalised sum x, = (€,+..-+,)//m, and p,(x) 
the density of the absolutely continuous part of F(x): 
If for at least one n, F,(x) is non-singular then p,(x) 
converges almost everywhere to d(x) = (2x) *exp(— 4x). 
This result is valid even if €,, €, .-. are random vectors 
with mean zero and dispersion matrix &, provided (x) 
is replaced by the corresponding normal density with 
mean zero and dispersion matrix &. 

It is pointed out that the above result can be extended 
to the case of convergence to stable laws. Generalising 
4 result of Gnedenko, the following estimate is obtained: 
there exist positive constants ¢;, C2, (20, <1) such 
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1.3) 


from this table that the total number of symmetrical 
patterns is only a small ratio of all patterns possible 
for a particular case, and further, that this ratio decreases 
with increasing m, according to a square root law. 
This small ratio leads to two generalisations of practical 
consequences. First, in any physical problem involving 
the elucidation of the pattern of the structure, it would 
be clearly advantageous to seek any evidence of sym- 
metry within the structure. If symmetry is displayed, 
not only is the number of patterns possible reduced 
significantly but also these can be listed and checked 
systematically under the groups and sub-groups of Gg, 
G,, and G,. Secondly, for p, m>2 greater than the 
total number of patterns (p/m?) is of the order of 


2 
1/8 @ though always greater than this number. 


An appendix giving proofs of the theorems presented 
is included. 


(R. E, Stoltz) 


1.5 (—) 


that L{x :| p,(x)—(x) | Sei/ Jn} S cA", L being 
Lebesgue measure. 

As an application of the above-mentioned result, the 
authors prove the following theorem in the case of 
random vectors taking values in k-dimensional Euclidean 
space E;: 

let the common distribution function F(x) of ¢,, 
€,,... be absolutely continuous, A(x), A(x) be two fixed 
independent linear functions on E,, p,(x | y) the density 
of the conditional distribution of A,(é,) given A2(¢,) 
and $(x | y) the corresponding conditional density of the 
normal distribution. Then limp,(x|y) = $(x|y) 
almost everywhere in x for almost all (Lebesgue) y. In 
particular, 


lim sup| P,(A| y)—P(A| y)1 = 0 
n>ao A 


for almost all y, where the supremum is taken over all 
Borel sets, and P,(A | y), P(A | y) represents the corres- 
ponding conditional distributions. 


(K. R. Parthasarathy) 
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ROZANOFF, Y. A. (Moscow) 1.5 (-.-) 
A central limit theorem for additive random functions—Jn Russian 
Teor. Veroyat. Primen. (1960) 5, 243-246 (4 references) 


The additive random functions H(A) of semi-interval 
A = [s, ft), satisfying the strong mixing condition: 


AEM, BEMA, 
a(t) = sup sup|Pr(AB)—Pr(A)Pr(B)|>0 (t-00) 
t A,B 


are considered. Let a(t) = O(t '~*), 6>2/e 
M | H(Ao)—MH(Ao) |?*° SMo, 
fomall A, = [to, t+ to); 
O0< lim a) < im Se) <0 
|A|72o0 <A JA] >a 
then for | A| = t—s—> 


H(A)—MH(A) 1 x 
De re ee ev2 du. 
eo el 


(Y. A. Rozanoff) 
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ADKE, S. R. (University of Poona) 2.0 (-.-) 
Generalised affinity and a class of distance functions—In English 
Proc. Nat. Inst. Sci. India A (1959) 5, 104-110 (7 references) 


Bhattacharya has introduced the concept of affinity 
between probability laws, which is a measure of diver- 
gence or distance between two statistical populations 
[Bull. Calcutta Math. Soc. (1943) 35, 99-109], providing 
a generalisation of Mahalanobis’ distance between two 
multivariate normal distributions. Some properties of 
Bhattacharya’s affinity and the distance functions 
associated with it have been studied by Adhikari & 
Joshi [Publ. Inst. Statist. Paris (1956) 5, 57-74]. The 
present author generalises the definition of affinity and 
considers some properties of a class of distance functions 
associated with the generalised affinity. 


(B. Raja Rao) 
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BLACKWELL, D. & HODGES, J. L., Jr. (University of California, Berkeley) 2.1 (1.1) 
The probability in the extreme tail of a convolution—/n English 
Ann. Math. Statist. (1959) 30, 1113-1120 (6 references, 2 tables, 1 figure) 


For the case of independent and identically distributed 
lattice variables, the authors give approximations to 
Pr{X,+...+X, = na} and Pr{X,+...4+X,2na} where 
&X,<a<sup X,. It is assumed that the moment 
generating function of X, exists in a neighbourhood of 
the origin. The derivation of the approximations 

involves defining, through the moment generating 
function, an auxiliary sequence of independent and 
identically distributed random variables Y,, ..., Y,, 
taking the same values as X,, so that €Y, = a. From 
the definition it follows that the “tail’’ probabilities 
above are related to ‘‘central”’ probabilities of the 
distribution of Y,+...+Y,. The approximations are 
then provided by an inversion argument in terms of the 
characteristic function of Y, and involve the central 
moments of Y,. Two approximations are given for 
each probability and for each case, the approximations 
have relative errors of order n”~1 andn~”. A conjecture 
is given concerning an analogous approximation in the 


non-lattice case. 


(N. D. Ylvisaker) 
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DABONI, L. (University of Trieste, Italy) 
A property of Poisson distributions—Jn Italian 


LT W#{L 


2.5 (1.1) 


Boll. Un. Mat. Ital. (1959) 14, 318-320 (1 reference) 


The author, with reference to the proof by Lamperti 
[The American Mathematical Monthly (1959)], gives a 
new proof of a theorem according to which if, in a 
series of Bernoulli trials their number N is random, 
then the numbers V4 of successes and Ng = N—N, of 
failures are stochastically independent only when N 
follows a Poisson distribution. 

The argument is based on the fact that the generating 
functions of Poisson distributions are the unique solu- 
tions of a functional equation of the generating function 
of N based on the condition that N, and Ng, are 
stochastically independent. 


(G. Chiassino) 
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EKMAN, G. (Inst. Math. Statist., Univ. Stockholm) 


2.1 (8.7) 


Approximate expressions for the conditional mean and variance over small intervals 


of a continuous distribution—Jn English 


Ann. Math. Statist. (1959) 30, 1131-1134 (2 references) 


The paper presents approximations to the conditional 
mean and conditional variance of an absolutely con- 
tinuous random variable, given that the random variable 
is restricted to a small interval of its range. For this 
purpose the density function f of the random variable 
is assumed to have continuous derivatives up to the 
fourth order. The approximations result from the 
combination of expansions for the first three ‘‘ moments ” 
of the random variable over the interval of restriction, 
in powers of the length of this interval, with the Taylor 
series for log f, applied to the same interval. The work 
stems from the problem of determining optimum 
stratification points for stratified sampling with intent 
to estimate the mean of a continuous random variable. 
Application of the results is made in this paper to the 
determination of stratification points when strata sample 
sizes are to be allocated in proportion to the probabilities 
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of the strata. The corresponding problem of selection 
of optimum stratification points to be used with a 
minimum variance linear estimate of the mean of a 
continuous random variable is discussed by the author 
in a previous paper [Ann. Math. Statist. (1959) 30, 
219-229; abstracted in this journal No. 1/268, 8.1]. 


(J. D. Esary) 


GEFFROY, J. 


2.6 (1.5) 


A contribution to the theory of extreme values—Jn French 
Publ. Inst. Statist. Paris (1959) 8, 36-121, 124-185 (26 references) 


The first chapter of this long and detailed paper is 
devoted to a study of the upper and lower asymptotic 
bounds of the extreme values of a sample. After an 
introduction including a section dealing with notation, 
the author investigates stochastic convergence, upper 
and lower probability bounds for extreme values, almost 
certain upper and lower bounds and almost complete 
_ upper and lower bounds. Twelve theorems are stated 
in this chapter which concludes with a note on an 
unusual class of functions. 

The study of stability in probability of extreme sample 
values forms the subject of the next chapter. The author, 
after stating that this study follows logically from the 
results of his first chapter, proceeds to discuss stability 
of random variables, inversion of a frequency function 
and the stability of the sum of random variables, and a 
theorem by Gnedenko [Amn. Math. (1943) 44, 423] and 
notes that this theorem furnishes a stability condition. 
The latter part of the chapter deals with the study of 
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GHURYE, S. G. & WALLACE, D. L. 


(University of Chicago) 


linked extremes of a sample from the point of view of 
stability and ends with the equivalence in probability 
of the extreme values of a sample. This chapter contains 
many theorems and their corollaries. 

Chapter three dealing with almost complete stability 
of the extreme values of a sample is followed in chapter 
four by a study of almost sure stability. 

In chapter five the author gives some consideration 
to the idea of limiting independence of two random 
variables and its application to the study of the stability 
of the “middle”? of the sample and range. In the 
following chapter a detailed study is made of certain 
limit laws and the paper concludes with an investigation 
of the properties of two-dimensional samples. 

In the two parts of this paper there are nearly eighty 
theorems and their corollaries together with an extensive 
bibliography. 


(Mile. Gervaise) 


2.1 (2.9) 


A convolutive class of monotone likelihood ratio families—In English 
Ann. Math. Statist. (1959) 30, 1158-1164 (5 references) 


If f(-, 0) and g(-, 0) define two one-parameter families 
of densities on the real line or the integers, which 
have the monotone likelihood-ratio property, the 
family defined by convolution (f*g) will not generally 
have the property. A sufficient condition for the con- 
volution to have the monotone likelihood-ratio property 
is that for every 0, the one-parameter families with 
parameter € defined by A(x, ¢)=f(x—¢, 0) and 
k(x, €) = g(x—é, 0) also have the property. The 
condition is shown not to be necessary. 

A corollary and application is: let Y be a real random 
variable and, conditionally on Y, let X,, ..., X, be in- 
dependent Bernoulli random variables. Assume that 
for each i, Pr{x; = 1| Y} is a nondecreasing function 
of Y. Then, with Y as a parameter, the family of 
distributions of X, has the monotone likelihood-ratio 
property for each i, and so does that of S = X,+...+%,. 
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Whatever the distribution of Y, the conditional distribu- 
tion function of Y given S = a lies to the left of that 
COGS — Dat. b a: 

Some generalisations of the results to very special 
multivariate distributions are given. 


(D. L. Wallace) — 


GUMBEL, E. J. (Columbia University, New York) 


2.0 (2.1) 


Multivariate distributions with given margins and analytical examples—In English 
Bull. Int. Statist. Inst. (1960) 37, III 363-373 (6 references) 


The aim of this paper is to derive a general and simple 
system which will connect any two given distributions 
F(x) and G(y) to a bivariate distribution H(x, y), such 
that the given distributions are marginal. This system 
is obtained by a generalisation of the construction due 
to Morgenstern [Mitt. Math. Statist. Wiirzburg (1956) 
8, 234-235]. The result, which is also generalised to 
any finite number of marginal distribution, is 
A(x, y) = F(x)G(y)1+a[1 — F(x) ][1 -GO)], 

where @ is any real number with —1<a<1. There is 
an infinity of bivariate distributions, such that the two 
given ones are marginal, which is certainly not exhausted 
by the given infinite system. Fréchet has shown, for 
example, that H,(x, y) = max[F(x)+G(y)—1, 0] and 
H,(x, y) = min[F(x), G(y)] are bivariate distributions 
for which, for all x and y, H,(x, y)S A(x, y)< A(x, y) 
for any A(x, y) having F(x) and G(y) as marginal 
distributions. 


2/255 


GURLAND, J. (Iowa State University, Ames) 


Included in the system is a bivariate distribution with 
uniform marginal distributions, for which the regression 
curves are linear, and a bivariate distribution with 
normal marginal distributions where the regression 
curves are not linear. In the system the coefficient of 
correlation varies only from —1/3 to,1/3: the limits 
are attained for uniform marginal distributions. Hence 
a correlation coefficient 0-30 means strong correlation 
for a bivariate distribution belonging to this system. 

Some further interesting examples of the unexpected 
behaviour of distributions are given. 


(J. T. Runnenburg) 


2.5 (4.0) 


Some applications of the negative binomial and other contagious distributions—In English 
Amer. J. Publ. Health (1959) 49, 1388-1399 (39 references) 


Much of the experimentation in the biological and 
medical sciences pertains to statistical distributions 
which are far from normal, for example the negative 
binomial and other contagious frequency distributions. 
The binomial and Poisson distributions are described, 
with an example, and their relation to the negative 
binomial discussed. If the Poisson distribution is com- 
pounded with a Gamma distribution, that is if the 
Poisson parameter / itself follows a Gamma distribution, 
the resulting distribution is negative binomial; this 
can be referred to as “‘ apparent contagion’”’. A model 
_ based on true contagion also yields the negative binomial 
as a limiting case, and so it is apparently impossible 
from the data to distinguish between the two models. 
Another mathematical model that leads to the negative 
binomial distribution arises when the Poisson distribu- 
tion is generalised through a logarithmic distribution; 
for example, if the number of bacterial colonies in a 
field follows a Poisson distribution and the number of 
bacteria per colony a logarithmic distribution, then the 
distribution of bacteria ineach field is a negative binomial. 
The method of moments and the method of frequencies 
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are described for estimation and fitting. Other com- 
pound and generalised distributions are discussed, and, 
along with the negative binomial, fitted to some medical 
data. 

One example concerns the distribution of dental caries 
in twelve-year-old children. A Pascal distribution 
compounded with a Poisson distribution (which is of 
the same form as a Poisson distribution generalised 
through a Pascal distribution) is fitted by a simple 
combination of the method of moments and method 
of frequencies; the fit is about as close as that accom- 
plished by the negative binomial distribution when the 
method of maximum likelihood is used: Another 
example concerns the sickness distribution of shunters 
over a period of time. Here the negative binomial 
gives a poor fit even when determined by maximum 
likelihood. Using the method of moments the best 
fit is obtained with a Neyman Type A distribution 
(which arises as a Poisson distribution either com- 
pounded with or generalised through another Poisson 
distribution). 


(R. C. Elston) 


SSS 


KHAMIS, S. H. (Amer. Univ., Beirut and Food & Agric. Org., Rome) 


2.7 (0.3) 


Incomplete gamma functions expansions of statistical distribution functions—JIn English 
Bull. Int. Statist. Inst. (1960) 37, IIL 385-396 (14 references). Discussion abstracted No. 2/258 


Incomplete gamma function series expansions of density 
and cumulative distribution functions are developed by 
the author in this paper. The coefficients of the successive 
terms depend upon the moments of the unknown 
distribution function. Finite difference methods are 
developed for the computation of the coefficients which 


result in very simple computational schemes. The 
inversion of the resulting series is also discussed. 
iaeeparticular Gf 90) = [o"/T(n))x*~ 1. e°™,.. for 


positive n, x and a, and if O,(x) = \ 0,(x)dx, the two 
10) 

equivalent series expansions for a density function f(x) 

given by 


a A,A"0,(x) and C 9) as B,A®,+, = 1(X), 


and for the corresponding cumulative distribution 
function given by 


De A,A’®,(x) and Mes B®, +(x) 
r=0 r=0 
are obtained: where the A, are easily computed from 


the leading differences of the ratios of the moments 
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KHAMIS, S. H. (Report of discussion on a paper by) 


Incomplete gamma functions expansions of statistical distribution functions—In English 


Bull. Int. Statist. Inst. (1960) 37, I 102-103 


Those taking part were Hotelling, Neyman and Martin. 

Hotelling observed that this kind of expansion was 
particularly suitable for distributions of positive definite 
quadratic forms. This series of expansions offered 
possible help in studying effects of moderate deviations 
from standard assumptions when standard methods 
were applied. 

Neyman said that the method of deducing these 
formulae could be reduced to the minimisation of a 
certain sum of squares with convincingly devised 
weights. This method suffered from the disadvantage 
that on frequent occasions it led to approximation 
curves having negative ordinates. 

The author replied that he agreed with Neyman as 
to the difficulties regarding negative ordinates; this 
occurred in all series expansions. The only method he 
knows which, once successfully applied, would lead to 
a non-decreasing probability function was the method 
developed by Hartley & Khamis [Biometrika (1947) 34, 
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of the unknown distribution to the corresponding 
moments of 9,(x) and the B, from the A,. 

In the above, the difference operators A” apply to the 
subscripts of 6 and ®. The finite difference properties 
are used in conjunction with the inverse Taylor Theorem 
for inverting the above expansions. 

The following are two by-products of the methods 
used in the paper: 

(i) The Laguerre polynomial L,(x) defined by 

(—d/dx)"[x"0,020)] = L,(%)8,(x) 
is shown to be related to the differences of 0,(x) and 
®,(x) by the identities 

ress I(n+r) ASO (os) i I(n+r) N’0,(X) 

: Tin)  — A®,_ (x) Tin) 9,(x) 

(ii) The finite difference properties of 0,(x) and 
®,(x) are useful in the tabulation of these two functions. 
The method developed by the author for this purpose | 
has been used by Hartley & Pearson, E. S. in computing 
a new y’ table [Biometrika (1950) 37, 313-325}. | 

The paper contains a number of illustrations of the — 
application of the expansion and inversion techniques — 
developed therein. 


(S. H. Khamis) 


2.7 (0.3): fi 


340-351]. This method had an arbitrary parameter | 
and the resulting approximation was satisfactory. 
However, in respect of series expansions no solution had | 
yet been found. 


(W. R. Buckland) 


KHATRI, C. G. & SHAH, B. K. (University of Baroda, India) 
Method of fitting certain curves of force of mortality over the whole range—IJn English 


J. Univ. Baroda (1959) 8, 33-40 (3 references) 


Two new curves which will adequately fit data on force 
mortality , are discussed. They are 


@) yy = A+B dj +Bd3; 
Gi) y, = A+Hx+B,d{+B,d3. 
; The method of fitting known as internal least-squares, 
introduced by Hartley [Biometrika (1948) 35, 32-45] is 
adopted. It is assumed that the values of y, are given 
for values of x in an arithmetic progression. For both 
curves, separate equations for estimation of the constants 
are obtained for the three cases: 
(i) x taking values 0, 1, 2, ..., ; 
(ii) x taking values —m, —(m—1), ..., 1, 0, 1, 2, 
(aan cand 
(ili) x taking values —m+4, —(m—1)+H..., -4, 4, 
2 1 
29 2299 mM—>. 
‘Various particular cases are discussed. The method is 
illustrated by an example. 


(T. V. Hanumantha Rao) 
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LORD, F. M. (Educational Testing Service, Princeton, N.J.) 


An approach to mental test theory—In English 


Psychometrika (1959) 24, 283-302 (37 references, 2 figures) 


Five different true-score models are presented. In each 
an error of measurement is defined as a difference 
between true-score and observed score, and has an 
expected value of zero in all circumstances. The models 
presented then take further assumptions or definitions 


to show which additional inferences about true-scores 


can be made. Models discussed are: 
the matched-forms, 
the rationally equivalent forms, 
the item sampling, 
the Gaussian errors of measurement, 
the binomial error. 

The author discusses the various models and the 
relationships existing between models. He further 
points out that all five of the true-score models yield 
good results when applied appropriately. The need is 
suggested for empirical studies comparing the results 
obtained under different models to compare their various 
advantages and disadvantages. 


(R. E. Stoltz) 
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2.1 (0.2) 


2.1 (2.8) | 


McFADDEN, J. A. (Purdue University, Indiana) 


2307) 


Two expansions for the quadrivariate normal integral—In English 


Biometrika (1960) 47, 325-333 (24 references) 


If X,, X,, X¥3 and X, are normally distributed correlated 
variables with zero means, the author is concerned to 
find P,, the probability that X¥,, X,, X; and X, are 
simultaneously positive, in two special cases. For this 
purpose the generalised tetrachoric series for P, [see, 
for example, Kendall, M. G., J. R. Statist. Soc. (1945) 
108, 93-141 and Kendall, M. G. & Stuart, Advanced 
Theory of Statistics (1958), London: Griffin] are not 
well suited for computation. The two series expansions 
for P, given in this paper do not suffer this disadvantage. 

The first case treated is the stationary Markoff case; 
that is where the inverse of the correlation matrix has 
zero elements except on the main diagonal and im- 
mediately adjacent to it. A series expansion is given 
which, when the correlation coefficients are positive, 
gives a good approximation to P,. 

The second case is when p,3 = ~i4 + Prog = 0. A 
series expansion is obtained which gives reasonable 
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RAJA RAO, B. (University of Poona) 


results provided none of the correlations is large. The 
connection between the two series is discussed, and it is 
shown that the solutions are connected by a simple 
relationship. 


(Florence N. David) 


2.0 (-.-) 


Properties of the invariant I,, (m-odd) for distributions admitting sufficient statistics—In English 


Sankhya (1959) 21, 355-362 (6 references) 


If f and f’ are the density functions for two alternative 
probability distributions of a random variable X then 
the expression 


Le 2 [leer ax 


has been defined by Jeffreys [Proc. Roy. Soc. A (1946) 
- 186, 453-461 and Theory of Probability (1948) Oxford: 
Clarendon Press] as a measure of the distance between 
the two probability distributions. In the case where / 
and f’ are densities for two members of a Koopman 
family of distributions, Huzurbazar [Biometrika (1955) 
42, 533-537] has given the exact form of J,, for even 
values of m as an explicit function of the parameters. 

The author shows that for odd values of m, [,, can 
be expressed as 


lk (fim —fiimy"dx —2 [owrarnyras. 
a ea 


If f and f’ are densities for two members of a Koopman 
family of distributions, the first integral in the above 
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expression is resolved as an explicit function of the 
parameters. 

The author has also shown this to be the case with 
the second integral for the following families: 


(i) Normal distributions with the parameters varying 
separately and simultaneously, 
(ii) Gamma distributions with the parameters varying 
simultaneously, 
(iii) Poisson distributions, and 
(iv) Binomial distributions. 


It is noted that the differential forms of J, for cases (iii) 
and (iv) do not lead to the a priori probability distribu- 
tions of the relevant parameters used by Jeffreys in 
connection with problems of inference regarding these 
parameters. The author has shown however, that the 
mth root of the differential forms of J,, for case (i) 
with the parameters varying separately, do lead to the 
a priori probability distributions of Jeffreys for these 
parameters. 


(P. K. Bhattacharya) 


a ee eee eee. &w«w!wW 


RUBEN, H. (Columbia University, New York) 


WaT. 


Deeb 


Probability content of regions under spherical normal distributions. I.—Jn English 
Ann. Math. Statist. (1960) 31, 598-618 (40 references) 


This is the first of a series of papers attempting to lay 
the groundwork for a relatively well-rounded theory of 
the spherical normal distribution. Many distributional 
problems in the theory of mathematical statistics may 
be regarded as particular instances of the general problem 
of determining the probability content of geometrically 
well-defined regions in Euclidean N-space when the 
underlying distribution is centered spherical normal and 
has unit variance in any direction. Specifically it is 
required that for a definite region R 

e ?**dx 


Pr(R):= (2x) 2” | 
xeR 

in which x’ = (X%, ..., Xy), though integrals of this 
form are rarely capable of being expressed in closed 
form using well-known functions. It is hoped by the 
author that this paper will provide a unifying thread 
and thereby help stimulate further research. 

The opening discussion reviews the methods of express- 
ing the spherical normal distribution. In the remainder 
of the paper the author reviews and discusses a number 
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RUBEN, H. (Columbia University, New York) 


Probability content of regions under spherical normal distributions. I.—Jn English 


of important distributional problems which are formally 
reducible to integral form. The geometric forms which 
are considered are: 


(i) Probability content of a half space, 

(i) Probability contents of centrally and non- 
centrally located hyperspheres, 

(iii) Probability content of a symmetrically and asym- 
metrically located hyperspherical cone, 

(iv) Probability content of a region bounded by a 
variety of revolution of dimensionality N—1 and 
species p, 

(v) Probability contents of symmetrically and asym- 
metically located hyperspherical cylinders, 

(vi) Probability content of a centrally situated 
ellipsoid, 

(vii) Probability content of a regular simplex. 


The “method of sections,’ which will be used 
frequently to deal with integrals in the sequel paper, is 
introduced. This consists of dividing up the definite 
region R by means of a series of parallel and adjoining 


continued 


ee LN 
continued 


Ann. Math. Statist. (1960) 31, 598-618 (40 references) 


(N—1)-flats and of the exploitation of the following 
fundamental property of the spherical normal distribu- 
tion of dimensionality N: The conditional probability 
distribution in any linear subspace of dimensionality 
N-k (k = 1, 2, ..., N—1) is itself spherical normal with 
dimensionality N—k and with variance in any direction 
equal to the variance of the original N-dimensional 

distribution. It follows that the probability content of 
the infinitesimal region intercepted by R between two 
parallel infinitesimally spaced flats requires integration 
over one dimension of a spherical normal distribution 
which has been integrated over N—1 dimensions. If, 
in addition, the section of each cutting flat is a region 
of the same geometrical type as R (for example, R an 
ellipsoid and the section an ellipsoid), with certain 
geometrical relationships being retained by the N—1 
dimensional figure, then the probability content integral 
becomes an integral recurrence relationship. 


(D. A. Roellke) 
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2.1 (1.5) 


A remark on the paper “‘ A moment inequality with an application to the central 


limit theorem ” by C. G. Esseen—In Russian 


Teor. Veroyat. Primen (1960) 5, 125-128 (4 references) 


It is proved that 


lim inf sup ./n | BOjar{ a 26 
n> oo) ail! <0. vx | oO 
0<a<o 


where F(x) is normal distribution function and F,(x) 
is a distribution function of a normalised sum of in- 
dependent identically distributed random variables, 


_ M\|é-Mé,| a, 
(a= —_,- ~~. [he constant: 1/Qx)? cannot be 


improved. 


(B. A. Rogozin) 
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SENTIS, P. 


2.0 (—.-) 


Frequency distributions and some of their applications—in French 
Publ. Inst. Statist. Paris (1958) 7, 17-91, 93-160 (7 references) 


The theory of random variables is well known if one 
possesses a denumerable infinity of random independent 
variables. This paper tries to provide some approach 
to the general case. 

In the first chapter the author supposes that all the 
variables which have the same fixed value can be 
expressed as certain functions of a set of variates for 
whom one can fix arbitrary values. The second chapter 
poses a situation in which the values of the variates 
are known and lend themselves to a statistical study. 

The third and fourth parts are interesting from the 
point of view of the calculation of the properties when 
one is able to define the o functions, additive and 
n-multiplications. It also generalises the concept of the 
series and of infinite products in the aggregate of indices 
which one is unable to count. 

In this present work the author introduces a new 
concept of the multiplication of a variate by means 
of a number, this allows the definition of the integral 
function of independent random values, and then the 
derivation of a Weiner-Lévy process. The generalised 
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infinite products are known to be similar to the expon- 
ential of an integral. That which is defined here enables 
a direct study of the known products of the matrices, 
of a nucleus integral for example. 

It is clear from parts five and six of this paper that 
the preceeding theories have numerous applications in 
the field of the calculation of probabilities and the 
statistic fits the Kolmogoroff equation tolerably well. 
The author illustrates his paper by practical problems 
which occur in certain paper factories. 


(Mile. Gervaise) 


TEICHER, H. (Purdue University, Indiana) 
On the mixture of distributions—In English 


2.8 (1.5) 


Ann. Math. Statist. (1960) 31, 55-73 (20 references) 


Let F = {F} be a family of distribution functions and 

let {u} be a class of measures defined on a Borel field of 

subsets of # such that u(F) =1 for every p of {pu}. 
The distribution function 


H = H(x) = | F(x)du(F) 
F 


is called the p-mixture of ¥. If yu varies over the class 
{u} then one obtains the class {H} of mixtures H. In 
this paper the author considers families F = {F(x; 
1, -++, %m)} Which depend on a finite number of para- 
meters. He treats 


(i) Convergence questions which occur when either 
of the sequences F,, or pw, tend (in the sense of 
weak convergence) to a limit. 

(ii) Mixture of additively closed families. 

Mixtures of additively closed families are shown to 
be closed under convolutions. Conditions which govern 
normality are given for the p-mixtures of normal 
distributions. It is also shown that the mixture, of 
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ZELEN, M. & SEVERO, N. C.. (Nat. Bur. Standards, Washington, D.C.) 
Graphs for bivariate normal probabilities—Jn English 


symmetric stable distributions with fixed exponent 
a(0<a<2) cannot be a symmetric stable distribution 
with exponent a. Similarly, the mixture of normal 
distributions with identical means cannot be normal. 
The author presents a number of complex theorems. 


(E. Lukacs) 


2.3 (11.1) 


Ann. Math. Statist. (1960) 31, 619-624 (5 references, 3 figures) 


In presenting this paper of three charts which will 
enable one to compute easily the bivariate normal 
integral to a maximum error of 1/100, the author states 
that this work was motivated by the recent work by 
Owen & Wisen on this subject and also the fact that 
there has been much activity recently dealing with the 
tabulation of the bivariate normal probability integral: 
see Owen & Wiesen “‘ A method of computing bivariate 
normal probabilities ’’ [Bell System Tech. J. (1959) 38, 
1-20] and the National Bureau of Standards, “ Tables 
of the bivariate normal distribution functions and 
related functions ” [N.B.S., Applied Maths. Series No. 50 
(1959) U.S. Government Printing Office]. 

The authors first discuss notation and useful formulas 
relating to the bivariate normal integral: it is then shown 
(following the work of Owen) that this integral may be 
expressed as a sum of the functions of the two para- 
meters only. This two-parameter function is then 
plotted by exhibiting its level curves in steps of 0-01. 
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The three charts do not cover the complete range of 
parameter values, but simple transformations are 
available for extending the given range. 

The application of the charts, which includes visual 
interpretation, is illustrated by three examples including 
one which involves a non-rectangular region. 


(D. H. Shaffer) 


HARTER, H. L. (Aero. Res. Labs., Dayton, Ohio) 
Expected values of normal order statistics—Jn English 
ARL Tech. Rep. 60-292 (1960), Wright-Patterson Air Force Base. xvi+32 pp. (18 references, 3 tables) 


3.8 (11.1) 


A brief history is given of the development of the theory 
of order statistics since the work of Karl Pearson 
[Biometrika (1902) 1, 390-399], and of past efforts to 
tabulate their expected values for samples from a normal 
population, beginning with a table of the mean range of 
samples of size nm = 2(1)1000, published by Tippett 
[Biometrika (1925) 17, 364-387]. A fuller account is 
given of the method of computation of a new five- 
decimal place table of the expected values of all order 
statistics for samples of size n from a normal population. 
The expected value of the kth largest observation in a 
sample of size m from a standard normal population 
(mean zero and variance one) is given by the equation 
n! 


(n—kK)\(k—1)! 
fi x[3—@(x) 6" *[44+ O(x)]"*b(a)dx, 


E(x, | ns) = 


where (x) =(22)-#e*!? and @(x)= o(x)dx. 


The expected value of the kth smallest observation is 
given by the same expression preceded by a minus sign; 
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HARTER, H. L. (Aero. Res. Labs., Dayton, Ohio) 


Expected values of normal order statistics—In English 


so that for a given value of n it is necessary to compute 
the expected values only for k = 1(1)[n/2]. This was 
done by numerical integration on the Univac Scientific 
(ERA 1103A) computer, for n = 2(1)100 and for values 
of n none of whose prime factors exceeds seven, up 
through n = 400. The resulting table, rounded to five 
decimal places, and accurate to within a unit in the 
last place, is included in the report. 

A discussion is given of a proposal by Blom [Statistical 
Estimates and Transformed Beta-Variables (1958), New 
York: Wiley] for approximating the ith normal order 
statistic (ith smallest normal deviate for a sample of 
size n) by means of the relation 


(xj) = 0-1 (——*_ }, 
n—2a+1 


where O(x)= |  (x)dx, with (x) = {1/(2n)?}{1/e#*}. 
Blom tabulated the value of a required to yield the 
correct value of &(x,) for i = 1(1)[n/2] when n = 2(2) 
10(5)20, and suggested the use of « = 3/8 as a com- 
promise value. Table II of the present report extends 


continued 


3.8 (11.1) 
continued 


ARL Tech. Rep. 60-292 (1960), Wright-Patterson Air Force Base. xvi+32 pp. (18 references, 3 tables) 


these results to include i = 1(1)25(5)100(10)[n/2] for 
n = 25, 50, 100, 200 and 400. When 2 is of this 
magnitude, the value 3/8 for « is too small. Since « is 
quite sensitive to changes in n, no single value of « is 
satisfactory for all values of n; it is possible, however, 
to do a fairly good job of estimating &(x,) by choosing 
one or two compromise values of « for each n. Such 
values are given in Table II for n= 2(2)10(5)25, 50, 
100, 200 and 400, along with regression equations to 
be used for intermediate values of n, and a discussion 
of the errors involved in these approximations. 

A discussion is given of actual and potential uses of 
the tables, foremost of which is transformation to 
standard normal scores of ranked data or of other 
data not conforming to the assumptions underlying the 
analysis of variance. 


H. L. Harter 


HARTER, H. L. & CLEMM, D. S. (Wright Air Development Center, Dayton, Ohio) 


al aT 


3.8 (11.1) 


The probability integrals of the range and of the Studentised range-probability integral, 
percentage points, and moments of the range—Jn English 


WADC Technical Report 58-484, I (1959) Aero. Res. Labs., Wright-Patterson Air Force Base. 


(19 references, 3 tables) 


A description is given of the computation and use of 
tables of the probability integral of the range, percentage 
points of the range, and moments of the range for samples 
from a normal distribution, all most extensive and 
more accurate than previously published tables. The 
introductory material includes a brief summary of the 
development of the theory and the computation of 
tables of the range as far back as the work of Tippett 
[Biometrika (1925), 17, 364-387], also an account of 
applications dating back to the demonstration by E. S. 
Pearson [British Standards Institution No. 600, (1935)] 
that the range is useful in statistical quality control. 
It is stated that the principal reason for computing the 
new table of the probability integral of the range was 
the desire to use it in computing more accurate tables 
of the probability integral and of the percentage points 
of the studentised range (see Volume II of this report; 
abstracted in this journal No. 2/273, 3.8]. The method 
of computation of the tables on the Univac Scientific 
(ERA 1103 and 1103A) computer is discussed, along 
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HARTER, H. L. & CLEMM, D. S. 


(Wright Air Development Center, Dayton, Ohio) 
The probability integrals of the range and of the Studentised range-probability integral, 


xii+ 139 pp. 


with the accuracy of the results. The probability 
integral and the moments of the range were obtained 
by numerical integration, employing the seven-point 
Lagrangian integration formula, and using double- 
precision arithmetic in the case of the moments. The 
percentage points of the range were obtained by inverse 
interpolation in the table of the probability integral. 
In the section on the use of the tables, an example is 
given of the use of percentage points of the range in 
tests of hypotheses concerning the standard deviation 
of a normal population. The table of moments of the 
range, like that of the probability integral, is useful 
primarily in computing other tables. 

The following tables are included: (1) an eight- 
decimal-place table of the probability integral of the 
(standardised) range, W = w/o, at intervals of 0-01 for 
samples of size n = 2(1)20(2)40(10)100; (2) a six- 
decimal-place table of percentage points of the range 
for the same values of n and cumulative probability 
P = 0-0001, 0-0005, 0-001, 0-005, 0:01, 0-025, 0-05, 


continued 


3.8 (11.1) 
continued 


percentage points, and moments of the range—Jn English 
WADC Technical Report 58-484, I (1959) Aero. Res. Labs., Wright-Patterson Air Force Base. xii+139 pp. 


(19 references, 3 tables) 


0-1(0-1)0-9, 0:95, 0-975, 0-99, 0-995, 0-999, 0-9995 and 
0-9999; and (3) a table of moments of the range [mean 
to 10 decimal places (11 significant figures), variance to 
10 decimal places (10 significant figures), skewness to 

8 decimal places (8 significant figures) and elongation 
(sometimes inappropriately called kurtosis) to 7 decimal 
places (8 significant figures)] for samples of size 
me=12(1)100. 


(H. L. Harter) 
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HARTER, H. L., CLEMM, D. S. & GUTHRIE, E. H. (Wright Air Development 


3.8 (11.1) 
Center, Dayton, Ohio) 


The probability integrals of the range and of the Studentised range—probability integral and percentage 
points of the Studentised range; critical values for Duncan’s new multiple range test—In English 


WADC Technical Report 58-484, II (1959) Aero. Res. Labs., Wright-Patterson Air Force Base. xxi+293 pp. 


(15 references, 3 tables) 


A description is given of the computation and use of 
tables of the probability integral of the studentised 
range, percentage points of the studentised range, and 
critical values for Duncan’s new multiple range test, 
all based upon samples from a normal distribution, and 
all more extensive and more accurate than previously 
published tables. The introductory material includes 
a brief summary of the development of the theory and 
the computation of tables of the studentised range 
from the time (1932) that “‘ Student ’’ (W. S. Gosset) 
first proposed the studentisation of the range. It is 
stated that the primary reason for computing the new 
table of the probability integral of the studentised range 
was the desire to use it in computing more extensive 
and more accurate tables of percentage points. More 
accurate tables were needed especially in the case of the 
special percentage points, for protection levels based 
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upon degrees-of-freedom, which are the critical values 
for Duncan’s new multiple-range test, for which 
previously published tables contain sizeable inaccuracies. 
The method of computation of the tables on the Univac 
Scientific (ERA 1103 and 1103A) computer is discussed, 


along with the accuracy of the results. The probability 


integral of the studentised range was obtained by 
numerical integration, employing the seven-point 
Lagrangian integration formula and using a new 8-place 
table of the probability integral of the range [see Volume 


I of this report: abstracted in this journal No. 2/271, 


3.8]. The regular percentage points and the special 
percentage points which are the critical values for 
Duncan’s test were obtained by inverse interpolation 
in the table of the probability integral of the studentised 
range. An example illustrates the use of the tables in 
performing two multiple comparisons tests on observed 


continued 
(Wright Air Development 3.8 (11.1) 
Center, Dayton, Ohio) continued 


The probability integrals of the range and of the Studentised range—probability integral and percentage 
points of the Studentised range; critical values for Duncan’s new multiple range test—Jn English 


WADC Technical Report 58-484, IT (1959) Aero. Res. Labs., Wright-Patterson Air Force Base. xxi+293 pp. 


(15 references, 3 tables) 


data, the regular percentage points being used for the 
Newman-Keuls test and the special percentage points 
for Duncan’s test. 

The following tables are included: 


(i) a table (to six decimal places or six significant 
figures, whichever is less accurate) of the prob- 
ability integral of the studentised range, O = w/s, 
with varying intervals (small enough to make the 
table interpolable) for samples of size n = 2(1) 
20(2) 40(10)100, with degrees of freedom v = 1(1)20, 
24, 30, 40, 60, and 120 for the independent 
estimate s of the population standard deviation; 

(ii) a table (to four decimal places or four significant 
figures, whichever is less accurate) of the per- 
centage points of the studentised range corres- 
ponding to cumulative probabilities P = 0-001, 
0-005, 0-01, 0-025, 0:05, 0-1(0-1)0-9, 0:95, 0-975, 
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0:99, 0-995, and 0-999 for the same values of n 
and y, also for v = 00; and 

(iii) a table (to four decimal places or four significant 
figures, whichever is less accurate) of the critical 
values for Duncan’s new multiple range test at 
significance levels « = 0:1, 0:05, 0-01, 0-05, and 
0-001, corresponding respectively to protection 
levels P = (0-9)"~ !,(0-95)"~ 1,(0-99)"~*, (0:995)"~ 4, 
and (0-999)""1, for the same values of nm and y 
(including vy = 00). 


(H. L. Harter) 


BARTON, D. E., DAVID, Florence N. & O’NEILL, Anne F. (University Coll., London 


32 1(0m) 


and Wheaton Coll., Norton, U.S.A.) 
Some properties of the distribution of the logarithm of non-central F—Jn English 
Biometrika (1960) 47, 417-424 (11 references, 4 tables) 


Fisher & Cornish [Rey. Int. Statist. Inst. (1937) 5, 307- 
321] showed that the distribution of the logarithm of 
central F can be approximated to be means of an 
Edgeworth expansion. 

The authors consider the logarithm of the non- 
central F, where non-central F is defined as the ratio of 
x [fy and x7/f, where y’?is a non-central y? with f, degrees 
of freedom and parameter 4 and ? is an independent y? 
with f, degrees of freedom. This logarithm they denote 
by z* and, following the procedure of Fisher & Cornish, 
they discuss expansions based on the assumption that 
an Edgeworth expansion is adequate. The moments 
of z* are given in terms of the polygamma functions 
plus a simple series expansion which is quickly con- 
vergent. 

Two series expansions are found to be of value. First, 
for A small, the distribution of z* is expanded about the 
appropriate point of the central z. Secondly, the 
Edgeworth expansion for z* is given. Since as many 
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moments of z* can be obtained as desired, the possibility 
of fitting a Johnson S, curve is also considered. 

Tables are given whereby the moments of z* can 
easily be evaluated. Comparisons are made between 
the several approximations suggested: [see also abstracts 
No. 2/278, 3.2 and No. 2/283, 3.2 in this present journal]. 


(Florence N. David) 


3.9 (0.6) 


On conditions for the forms of the type XAX’ to be distributed independently or to obey 


Wishart distribution—In English 


Bull. Calcutta Statist. Ass. (1959) 8, 162-168 (4 references) 


Let X be a p Xn random matrix formed by independently 
distributed column vectors. Let the jth column vector 
have the multivariate normal distribution with mean 
vector p; and dispersion matrix &. What are the 
conditions that a nxn symmetric matrix A should 
satisfy in order that XAX’ should have the noncentral 
Wishart distribution? The author proves that if A is 
a matrix whose rank is greater than p, then a necessary 
and sufficient condition is that A should be idempotent. 

A second problem which the author considers is the 
following: what are the conditions which the two nxn 
symmetric matrices A and B should satisfy in order 
that the forms XAX’ and XBX’ should be distributed 
independently of each other? The author shows that a 
necessary and sufficient condition for this is that AB = 0. 
From this the following corollaries are deduced: 


(i) Let A be an nxn matrix and L an nxt matrix. 
A necessary and sufficient condition that XAX’ 
and XL are independently distributed is that 
AL. =O; 
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Gi) Let" A; =el% 2,.%.,.m) be m, nxn symmetne 
matrices. A necessary and sufficient condition 
that the forms XA;X’ (i= 1, 2, ..., m) be in- 
dependently distributed is that 4,4; = 0(@ # J). 


Using the results mentioned above the author also 
obtains theorems similar to Cochran’s theorem on the 
distribution of quadratic forms. 

Finally the following results concerning the distribu- 
tion of X’AX (as distinct from XAX’) are stated without 
proof. 

(i) Let A be a pXp matrix of rank greater than n. 
A necessary and sufficient condition for X’AX 
to be Wishart is that ALA = A that is, AX should 
be idempotent. 

(ii) A necessary and sufficient condition for X’AX 

and X’BX to be independently distributed is that 
AXB = 0. 


(S. John) 


, Rie CaO RON 
rile ls: 


LAHA, R. G. & LUKACS, E. (Catholic University of America, Washington, D.C.) 3.4 (0.6) 
On certain functions of normal variates which are uncorrelated of high order—JIn English 
Biometrika (1960) 47, 175-176 (5 references) 


Necessary and sufficient conditions are given such that 
an inhomogeneous quadratic form in n independent 
random normal variables shall be independent of a 
linear form in these same variables. Thus if Q = xAx’+ 
bx’ and L = ex’ are the quadratic and linear forms 
respectively where x is a row vector of n independent 
normal variables and A, b and ¢ are matrices of constant 
coefficients, then Q and L are independent if and only 
if cA =0 and cb'=0. If &(O'L) = &(OE(L), 
(ae et a. Ss) O and L. arevsaid to be un- 
correlated to order (r, s). It is shown that if Q and L 
are uncorrelated to order (2,2) then cA = cb’ = 0. 
If Q is homogeneous quadratic, then the conditions for 
independence derived by Laha & Lukacs are weaker 
than those conditions which may be derived from a 
result by Kawada [Amn. Math. Statist. (1950) 21, 614- 
615] which require that Q and L shall be uncorrelated 
to order (2,4). 


(J. G. Saw) 


BOT, 


PEARSON, E. S. (University College, London) 3.2 (-.-) _ 
Editorial Note on non-central F approximations—ZJn English | 
Biometrika (1960) 47, 430-431 (1 reference, 1 table) | 


The author, the Editor of Biometrika, compares and 
contrasts the approximations to the probability integral 
of non-central F suggested by Severo & Zelen [Bio- 
metrika (1960) 47, 411-416: abstracted in this present , 
journal No. 2/283, 3.2] with those put forward by 
Barton, David, F. N. & O’Neill [Biometrika (1960) 47, 
417-429: abstracted in this present journal No. 2/275, 
3.2]: 

He shows that, over a wide range of values of /;, 
f, and 4d, there is little to choose between the two sets of 
approximations. Further, it appears that the results 
obtained by using a Johnson S, curve [Biometrika 
(1949) 36, 149-176] are on the whole better than either 
set. 


(Florence N. David) 
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PRANAB, K. S. (Calcutta University) 
On the moments of sample quantiles—Jn English 


3.8 (1.5) 


Bull. Calcutta Statist. Ass. (1959) 9, 1-19 (3 references) 


This paper deals with the existence of moments of the 
ordered statistics and the convergence of the moments 
of the sample quantiles as the sample size tends to infinity. 

The rth ordered statistic (the rth smallest observation) 
in a sample of size n is denoted by x,,)._ The pth sample 
quantile is defined to be the ordered statistic x (fnp}+1) 
and the pth quantile of the population is denoted by 
€p- Throughout this paper it is assumed that the 
population distribution function F(x) admits of a con- 
tinuous density function f(x). The main results of this 
paper are the following: 

Theorem 1. If the population admits of the dth 
moment with 6>0 then the kth moment of the rth 
ordered statistic exists whenever r satisfies the inequality 
ro SrSn—ryo—1 where dr =k. A consequence of 
Theorem | is that the kth moment of a sample quantile 
will exist in sufficiently large samples. 

It is known that if f(x) is non-zero at €,, then the 
sample pth quantile has an asymptotic normal distribu- 


tion. Concerning the convergence of moments the 
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author states the following theorem: 
Theorem 2. Let 


(i) the population admit of a é6th moment for some 
o>0, 

(ii) f(x) be non-zero at €, and 

(iii) f(x) have a bounded derivative in a neighbourhood 

OL cs: 
then the moments of the standardised sample pth 
quantile converge to those of the normal distributions. 
Theorem 2 is a generalisation of a result due to Chu & 
Hotelling [Ann. Math. Statist. (1955) 26, 593-606]. 

In the concluding section the author indicates how, 
in particular situations, the conditions can be relaxed 
and states, without proving, the following corollary: 

Corollary 1. For any distribution with a continuous 
density function and finite dth moment with 6>0, the 
central moments of the pth sample quantile converge 
to zero provided the population pth quantile is uniquely 
defined. 


(J. Sethuraman) 


3.6 (3.8) 


Variance of the median of small samples from several special populations—Jn English 
J. Amer. Statist. Ass. (1960) 55, 148-150 (4 references, 7 tables) 


The distribution of the median of a sample of size n 
from a population having the density function f(x) is 
asymptotically normal with mean m and variance given 
by the reciprocal of 4n[f(m)]’, where m is the population 
median. This asymptotic variance is compared with 
the exact variance of the median for small samples 
from several special populations. 
_ The populations considered were exponential, normal, 
rectangular, cosine, parabolic and inverted parabolic. 
Samples of size n = 1, 3, 5, 7 were taken for comparison. 
For each population and for each sample size a table of 
exact variance, asymptotic variance and relative error 
is given. 

It is found that for each population the error decreases 
monotonically as the sample size increases. The relative 
error decreases among the six populations as the 
standardised fourth moment increases, so that the 
asymptotic variance is most accurate for the exponential 


population, followed in order by normal, cosine, para- 
bolic, rectangular and inverted parabolic. Finally, it 
is found that the asymptotic variance exceeds the exact 
variance except for the exponential population. 


(S. Krane) 


ROY, A. R. & MOHANTY, S. G. (Indian Council of Agric. Res., New Delhi) 


3.1 (5.6) 


Distribution of y?-analogue for normal population with class intervals 


defined in terms of the sample median—Jn English 


J. Indian Soc. Agric. Statist. (1958) 10, 90-98 (4 references) 


The authors consider independent observations eine oe 
-..» X, from N(@, 1) and the goodness-of-fit chi-square 
statistics 


R, = y (O;—E,)/E; 


where, 0,, 0, and 03 are the frequencies of observations 
in intervals (—0oo, X+a], (X+a, ¥+b), [X+b, «); & 
is the sample median, a and 5b are constants (a<b), 
E, = nn, (i = 1, 2, 3); and 

a 


m, =(2n) * | exp(—41?)dt, 


— 0 


son | 


a 


exp (—4t?)dt, 


Ray On | exp (—4t?)dt. 
b 
It is seen that E,, E, and £3 are expected frequencies 
in the non-stochastic intervals (— oo, +a], (0+a, 6+5), 
[9+b, oo). The authors show that in the limit, as 
noo, R, is distributed in the same way as 4,y7+/,y3 
where y, and y, are independent standard normal 
variables and 4,, 4, are constants depending on a, b 
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but not on 0. For each a, b (a<b), 4,<1 and 1,21, 
the equality sign holding only when a+b = 1. 

The proof runs as follows: consider the intervals 
(—oo, 0+a], (0@+a, 0+b), [0+b, «) and define b,(x) 
(i = 1, 2, 3) as the function which takes the value one 
if x is in the ith interval and otherwise the value zero. 
It is demonstrated that for each i = 1, 2, 3, 


(1/,/n)(O;—E;) = ./n[U; + c(X—0)]+0,(1) 
where U; = (1/n) ) [b(x,)—7;] 


Che (Qn) 7 c= (On)etles Pers oy 

C3 = —(2n) te7 2. 
This makes use of the limiting multivariate normal 
distribution of the median and U-statistics obtained by 
Sukhatme [J. R. Statist. Soc. B (1957) 19, 144-148], 
to establish a similar property of (O;—E;) (i = 1, 2, 3). 
The main bulk of the algebra is in the computation of 
the asymptotic moment-matrix of the variables E?(O;— E,) 
and in deducing the properties of its two non-zero latent 
roots A, and 4). 


(S. K. Mitra) 


3.8 (3.9) 


Coun. Agric. Res., New Delhi) 
Joint asymptotic distribution of U-statistics and order statistics—In English 


Sankhya (1959) 21, 289-298 (3 references) 


It is shown under some mild restrictions that the joint 
distribution of a U-statistic, (see Hoeffding [Ann. Math. 
Statist. (1948) 19, 293-325]), and the a,th order statistic, 
computed from n independent observations on a random 
variable X, both suitably standardised, tends to 
(i) the bivariate normal distribution, if @,/n—p, 
Opals 

(ii) the joint distribution of two independent variables 
one of which is gamma and the other normal, 
if a,2K or n—a,—>K where K is a constant not 
depending on n; 

(iii) the joint distribution of two independent normal 
variables, if a,—0o such that a,/n-—0 or 1. 

The proof requires the following assumptions: 

(a) The U-statistic is generated from a bounded 
symmetric kernel. 

(b) The random variable X has a continuous distribu- 
tion function F(x), and in case (i) it is further 
assumed that 

(c) F has a density which is continuous at the pth 
quartile. 
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An earlier paper by one of the authors, Sukhatme 
[J. R. Statist. Soc. B (1957) 19, 144-148], proves case (i) 
for the special value of p = 1/2. The proof in the 
present paper proceeds on similar lines and consists 
mainly in evaluating the limiting expression for the 
characteristic function of the two standardised statistics. 

It is stated that in proving case (i) the assumption (a) 
could be somewhat relaxed and replaced by the following 
weaker assumptions. 


(al) The U-statistic is generated from a symmetric 
kernel which is bounded on every bounded 
interval of its arguments. 

(a2) The kernel considered as a random variable with 
its arguments replaced by independent observa- 
tions on X, has a finite third moment. 


The results (i)-(iii) are generalised several U-statistics 
and several order statistics: analogous results are true 
for the generalised U-statistics of Lehmann [Ann. Math. 
Statist. (1951) 22, 165-179] computed from independent 
observations on several random variables. 


(S. K. Mitra) 


sos 


SEVERO, N. C. & ZELEN, M. (University of Buffalo and National Bureau of Standards) 


3.2 (3.3) 


Normal approximation to the chi-squared and non-central F probability functions—Jn English 
Biometrika (1960) 47, 411-416 (7 references, 5 tables) 


The authors derive an approximation to the probability 
integral of the x? distribution based on the Wilson- 
Hilferty cube root transformation. Tables are given 
of the various factors which intervene in this approxima- 
tion. 

The non-central F may be defined as (7'?//;)/(y3//s) 
where x" is a non-central y? with f, degrees of freedom 
and parameter A, and x3 is an independent central ide 
with f, degrees of freedom. Abdel-Aty [Biometrika 
(1954) 41, 538-542] suggested that (73°/f,+4)* = Yj 
might be approximately normal and this is the result 
which the authors use. 

Writing Y3 for the denominator of the non-central F 
and c for the cube root of ( f; +4)/f,, the authors remark, 
following the procedure suggested by David & Johnson 
[Biometrika (1951) 38, 43-57], that a probability state- 
ment about non-central F can be reduced to finding the 
probability that L= cY{—(F’)*Y, is less than or equal 
to zero. JL is approximately normal with mean (,) 
and variance (o,) deducable from Abdel-Aty’s work 
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and Pr{L<O} = Pr(—y,/0,). Consequently estimation 
of the non-central F distribution can be thrown back 
to tables of the normal integral. 

Some comparison is given with the exact values of the 
probability integral calculated by Tang [Stat. Res. Mem. 
(1938) 2, 126-149] and with the approximation put 
forward by Patnaik [Biometrika (1949) 36, 202-232]. 
It is clear that the approximation suggested by the 
authors is good over a wide range of values of f,, f, 
and A. 

Further comparison are given by E. S. Pearson 
[Biometrika (1960) 47, 430-431: abstracted in this 
present journal No. 2/278, 3.2]; see also Barton, David 
F. N. & O'Neill [Biometrika (1960) 47, 417-429: 
abstracted in this present journal No. 2/275, 3.2]. 


(Florence N. David) 


3.7 (4.3) 


On the standard error of the Lorentz concentration—ZJn English 


Sankhya (1960) 22, 371-378 (2 references) 


In this paper the standard error of the Lorentz con- 
centration ratio has been estimated on the assumption 
of the lognormality of the distribution of incomes. 
Exact and large sample tests are then suggested for 
testing the equality of several concentration ratios. 

If log x has a normal distribution with mean @ and 
standard deviation A, then the concentration ratio L, 
- defined as twice the area between the egalitarian line 
and the Lorentz curve is given by 


Bl 
Nes 


L is independent of the location parameter 6. Sub- 
stituting 2, the maximum likelihood estimate of 2 in 
the above equation, the maximum likelihood estimate 
L of Lis obtained. ZL has an asymptotic normal distribu- 
tion with mean L and standard error (2mn) * exp(—447) 
where n is the number of observations. Substituting 
A for A, an estimate of the standard error is obtained. 

Usually, however, only grouped data are available. 


en dt. 


L= 24(x//2), $(8) = 


The maximum likelihood estimate of 2 is then tedious 
to compute. A simpler method based on the empirical 
Lorentz curve can be adopted. In the latter case, the 
estimate of the standard error obtained from the above 
expression will be an under-estimate. If the maximum 
likelihood method of estimation for grouped data is 
used, as is the practice in probit analysis, the above will 
give a close lower bound to the true standard error. 

The test for equality of several concentration ratios 
based on independent samples reduces to Bartlett’s test 
for equality of variances. In the special case of testing 
equality of two concentration ratios, the variance ratio 
can be used. If the sample sizes are large, the difference 
log 2,—logi, can be tested by the usual normal 
approximation. 

Finally, an empirical investigation is made to test the 
equality of two Lorentz ratios obtained from two inter- 
penetrating sub-samples of the National Sample Survey, 
India. 


(R. N. Bhattacharyya) 
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TURIN, G. L. (Hughes Research Labs., Culver City, California) 


3.4 (2.1) 


The characteristic function of Hermitian quadratic forms in complex normal 


variables—Jn English 
Biometrika (1960) 47, 199-201 (9 references) 


Wooding [Biometrika (1956) 43, 212-125] has given an 
expression for the multivariate distribution of several 
complex normal variables of a special type often 
encountered in the analysis of random noise. Using 
Wooding’s results the author finds the characteristic 
function of the real quadratic form in these variables 
and thence, by a Fourier transformation on the char- 
acteristic function, the density function of the real 
quadratic forms. 

The distribution is said to arise in, for example, the 
computation of the probability that an ideal binary 
hypothesis-testing communications receiver will err in 
its decision concerning which of two signals was trans- 
mitted through a channel comprising both a random 
transmission medium, such as the ionosphere, and 
additive thermal noise. 


(J. G. Saw) 
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WALLACE, D. L. (University of Chicago) 


3.1 (3.3) 


Bounds on normal approximations to Student’s and the Chi-square 


distributions—IJn English 


Ann. Math. Statist. (1959) 30, 1121-30 (6 references, 2 tables) 


Formulae closely related to 

u(t) = [nlog, (1+¢?/n)]* 

wt?) = [x —n—n log. (x’/n) 
are considered for converting upper tail values of 
Student’s ¢ or chi-square variates with nm degrees of 
freedom to normal deviates. Bounds on the exact 
normal deviate are constructed that are valid for any 
fixed n uniformly in the entire upper tail. Among the 
bounds, the following are the simplest. If x,(¢) is the 
exact normal deviate corresponding to ¢, then for any 
n=1 
u(t) [1 —(0-293/n)] Sx, (1) Sut?) 
u(t) — (0-368/./n) S$ x,(¢). 

If y,(7v’) is the exact normal deviate corresponding to 
y?, then for any n=1, and any y7>n 


w(t?) S In?) S w(W7) + (0-60/,/n). 
Closely related to the bounds for Student’s ¢ is the 
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approximation 


(+) u(t) 

8n+3 

to the normal deviate x,(f). It is claimed to be accurate 
to within 0-02 for t?/n<5 and to be superior to common 
approximations (even to the Paulson approximation). 
A more complex approximation with wider range of 
accuracy is given. 

For chi-square, the Wilson-Hilferty approximation 
is much superior to the bounds as approximations — 
except in the extreme tail and the chief value of the 
approximation is the uniform bound of order ./n on 
the error in the tail which is useful in theoretical work. 


(D. L. Wallace) 


COHEN, A. C., Jr. (University of Georgia, Athens, U.S.A.) 4.3 (2.5) 
Estimating the parameters of a modified Poisson distribution—Jn English 
J. Amer. Statist. Ass. (1960) 55, 139-143 (9 references, 1 table) 


The estimation of the parameter of a Poisson distribution 
is considered for the special case in which there exists 
a (possibly) non-zero probability that any value of 
* one ”’ is reported as “* zero’. The maximum likelihood 
estimate of the Poisson parameter is found as the unique 
positive root of a quadratic equation except for de- 
generate cases. The probability of misreporting ‘‘ one ”’ 
as “ zero” is also readily estimable by maximum likeli- 
hood. Asymptotic variances and covariance of these 
estimates are obtained by inversion of the information 
matrix. As an example, the author alters the well- 
known Bortkiewicz data on deaths from horse kicks in 
the Prussian army. 


(S. Krane) 


2/287 


CORNISH, E. A. (Div. of Math. Statist., CSIRO, Adelaide) 4.2 (3.9) 


The simultaneous fiducial distribution of the parameters of a normal bivariate distribution 
with equal variances—In English 


CSIRO Div. of Math. Statist., Tech. Paper, No. 6 (1960) 1-8 (7 references) 


In his Statistical Methods and Scientific Inference Fisher 
[(1959), Edinburgh: Oliver & Boyd] has derived the 
fiducial distribution of the parameters in a bivariate 
normal population having means €,, €,, variances of, 
o% and correlation coefficient p. This paper is concerned 
with the corresponding fiducial distribution for a bi- 
variate normal population with means ¢,, ¢,, equal 
variances o*, and correlation coefficient p. 

After obtaining the distribution of the sample statistics 
X,, X2, S, r, the author derives the simultaneous fiducial 
distribution of €,, €2, p, o. 


(J. Gani) 
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CZEN, P. & DZAN, D. I. (Inst. Math., Polish Acad. Sci., Wroclaw and Wroctaw Univ.) 


4.0 (8.7) 


On the division of the number of observations by using the ratio-delay method—In English 


Zast. Mat. (1960) 5, 179-194 


If we estimate the ratio of working-time of workers on 
the basis of n observations taken at random moments 
and we use Rubin’s estimator p = (m+ s/n/2) (n+ Jn), 
where m is the number of positive observations, then 
the mean-square error of estimation does not depend 
on the true working-time ratio p and we have for the 
expectation of the square of the error the relation 
E(p—p) = 1/41 + Jn)’. 
Imagine that a factory works in k shifts; that we are 
interested in estimating working-time ratio for each 
shift separately; that the loss caused by error of estimat- 
ing working-time ratio is proportional to the square of 
the error with coefficient c; for the ith shift and that 
we can take nm observations in all. If we then take n; 
observations in the ith shift, n,+,+...+m =n, then 
the expected loss will be 
ci 
1 4(1+./n))? 


&(S) = 


i 


I M= 
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The authors discuss the problem of optimal division 
of the total number of observations between the k shifts 
in order to minimise the expected loss &(S). They find 
an approximate solution of the form bn+a;, where the 
b’s add to unity and are proportional to the square-roots 
of the corresponding c’s. The error of this approximation 
is estimated for k = 2 and k = 3. Some examples are 
calculated and compared with the exact solutions. 


(S. Zubrzycki) 


4.1 (4.4) 


Confidence intervals in estimation of parameters—In Ukrainian 
Dopoy. Acad. Nauk Ukrain. SSR (1959) 355-358 (2 references) 


The author considers the expansion of the confidence 
interval x of the maximum likelihood estimate 0 of the 
parameter 0 by 1/./n powers (where n is the number 
of independent observations of a random variable € 
with a distribution function belonging to the family 
F(x, 0); that is to say x is found in the form 


x X2 
X=Xot + +H... 
pacer 


where x satisfies the correlation 
Pr NG Oe! <x} = 1-4 
Ky 


(0<a<1), where s? denotes the variance of 6. 


(A. I. Dorogovtzev) 
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DUNN, Olive J. (University of California, Los Angeles) 
Confidence intervals for the means of dependent, normally distributed random 


variables—In English 


4.2 (4.8) 


J. Amer. Statist. Ass. (1959) 54, 613-621 (6 references, 4 tables) 


The author considers the problem of constructing k 
simultaneous confidence intervals for the means of k 
correlated normal variates, based on n (k-variate) 
observations. Intervals having bounded confidence are 
given for arbitrary unknown covariances and for vari- 
ances known, unknown, and unknown but equal. The 
mathematical development was given in an earlier paper 
[Ann. Math. Statist. (1958) 29, 1095-1111: abstracted in 
this journal No. 1/45, 4.2]. 

In the present paper the application of the results is 
discussed, some tables are given to facilitate calculations, 
and the methods are illustrated by worked examples. 
The most generally satisfactory solution from the stand- 
point of weakness of assumptions and shortness of 
intervals is one based on a Bonferroni inequality, namely 
Pr(IIA;)=1— = Pr(4,), which allows the joint confidence 
level to be bounded in terms of individual confidence 
levels, independently of the correlations. 
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(Math. Inst., Polish Academy of Sciences, Warsaw) 


Table 615, for use with the Bonferroni solution, gives 
the 1-0-05/2k point of Student’s distribution for 
k = 1(1)10, 15, 20, 50, and for degrees of freedom, 
vy = 5, 10, 15, 20, 24, 30, 40, 60, 120, 00. In the three 
following tables v is the same, and k = 1(1)8. Table 616 
gives the 95 per cent. point of the Studentised maximum 
modulus, is based on work of Pillai & Ramachandran, 
and applies when variances are assumed equal. Table 
619a gives {kFy o5(k, v)}? (F = Snedecor’s variance 
ratio), applies when variances are equal, and is based 
on Scheffé’s simultaneous confidence intervals for linear 
contrasts. Table 619b gives {kvFo.95(k, v—k+1)/ 
(v—k+1)}, applies for arbitrary unknown variances, 
and is based on Hotelling’s 7. 


(R. J. Buehler) 


4.3 (3.8) 


Estimation of distribution parameters by means of statistics from an ordered sample—ZJn Polish 


Zast. Mat. (1960) 5, 149-154 


Let xf<x3S...<x;* be an ordered sample of size n 
from a population with a symmetric distribution. If 
n is even, then denote by x the mean of the values 
Xi}, «++» Xt. and by y the mean of the values X(i2) 41> +++» 
x*_441- In other words, we remove from the ordered 
sample the k—1 largest and the k—1 smallest values 
and then we compute the means of the smaller and 
larger halves of the remaining part of the sample; k 
being the parameter of the procedure. Similarly, if n 
is not even, then x is the mean of the values xj, ..., 
x* 4/2 and y is the mean of the values xf(-1)/21+2> ++» 
oot 

The authors are concerned with using (x+y)/2 as 
the estimate of the mean, and (y—~)/r,,, as the estimate 
of standard deviation in a normal population and they 
determine the values r,, rendering the last estimate 
unbiased and give a table for n = 3(2)9 and k = 1(1)4. 
They also find the coefficients in a linear form in x and y 
so as to get unbiased estimators of a and b for a uniform 
distribution over the segment a<tSb. 


(S. Zubrzycki) 
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HAJEK, J. (Math. Inst., Czechoslovak Academy of Sciences, Prague) 


4.9 (8.2) 


Optimum strategy and other problems in probability sampling—JIn English 


Cas. pést. mat. (1959) 84, 387-424 (13 references) 


We have an universe S consisting of arbitrary elements, 
that is, of numbers 1, ..., N. We select a subset s from 
S by a random experiment, so that every sc S has a 
known probability Pr(s). The selected subset will be 
called a sample, and the system of probabilities {Pr(s); 
sc S} will be called a sampling strategy. Let y,, ..., yy 
be some unknown values associated with the elements 


N 
1, ..., N, the total Y= )’ y; which is to be estimated 
i=1 


on the basis of values y; ascertained on the sample, that is 
for ies. Anestimator Yof the form Y= )\ y,w((s), where 


les 

ws)(ies, s<S) are arbitrary weights, depending not 
only of 7 but also of s, is called a linear estimator. The 
system of weights {w,(s), ies, scS} will be called an 
estimating strategy. A general estimator is of the form 
t = t(s, y), where the observation (s, y) involves the 
knowledge of s and of the values y; fories. 

A linear estimate is called representative with respect 


N 
S Zi; 
i=1 


In the paper it is 


to values z,, ..., Zy if the equation ) z,w,(s) = 


ies 


is fulfilled with probability one. 
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HAJEK, J. (Math. Inst., Czechoslovak Academy of Sciences, Prague) 
Optimum strategy and other problems in probability sampling—Jn English 
Cas. pést. mat. (1959) 84, 387-424 (13 references) 


Cochran and Godambe. Primarily, it is supposed that 
Vi, «+» Yn are realisations of independent random 
variables with mean values My; =n, and variances 


N 
M(y;—n,)? = d;;, and costs equal », cM; where 7; 


is the probability of inclusion of the ith element in the 
sample. It is shown that a strategy is optimum, if the 
linear estimator is firstly, representative with respect to 


values n,, ..., Ny, Secondly, is of the form » y,/m; and 


thirdly, the z,’s are proportional to numbers Jdiile :). 
This result is illustrated by examples and the aspects 
of its application are discussed. 
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shown that the mean-square deviation M( Y— Y)* for 
an Y, which is representative with respect to z,, ..., Zw, 
may be expressed as a weighted sum of the expressions 


\D 
(2 - %) . A further result concerns a use of sufficient 
Z; YE, , 


w] 

statistics in probability sampling: if s’ is a ‘‘ sample ”’ 
which involves not only the information as to which 
elements have been selected but also the order and 
number of times every element has been selected, and 
we consider the observation (s’, y), then the observation 
(s, y), where s has been defined above, is a sufficient 
statistic. 

The paper contains a detailed discussion of linear 
estimates for rejective sampling, consisting of n in- 
dependent drawings of one element with probabilities 
1, ---» &y provided that results where an element is 
selected more than once are rejected. A further method 
of sampling with varying probabilities, called permutation 
sampling, is presented. 

The main results concern optimum sampling and 
estimating strategy formulated in lines of papers by 


continued 


4.9 (8.2) 
continued 


It is then supposed that y,, ..., vy is a realisation of a 
stochastic process with stationary correlation function 
and with stationary coefficients of correlation, and it is 
shown that the best strategy is given by systematic 
sampling with generally unequal probabilities, which 
generalises and strengthens the well-known result by 
Cochran. 


(J. Hajek) 
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LATTER, B. D. H. & ROBERTSON, A. (Institute of Animal Genetics, Edinburgh) 
Experimental design in the estimation of heritability by regression methods—In English 


4.3 (3.5) 


Biometrics (1960) 16, 348-353 (2 references, 1 figure) 


Broadly speaking, there are two different kinds of tech- 
nique for the estimation of heritability: analysis of 
variance and regression coefficients. In some. special 
cases the analysis can be cast in either form. The 
problems of design in analysis of variance technique 
have been discussed recently by Robertson [Biometrics 
(1959) 15, 219-226: abstracted in this journal No. 1/289, 
9.4]. 

This paper is concerned with experimental design in 
regression methods; in particular that of offspring on 
parent. The authors state that the problem is to obtain 
the most accurate estimate of the regression coefficient 
and of the heritability that derives from it for a given 
expenditure of effort. This effort may be perhaps in 
time of measurement or in money expended in the 
rearing of the animals. The assumption made is that 
if n offspring are to be measured in each family then 
the expenditure of effort can be considered as pro- 
portional to n plus a constant, where the constant is 
dependent on the particular circumstances. The 
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LATTER, B. D. H. & ROBERTSON, A. 


(Institute of Animal Genetics, Edinburgh) 
Experimental design in the estimation of heritability by regression methods—In English 


problem and primary concern then is to find the optimum 
value of n under different conditions. 

The authors begin by giving the sampling variance 
of the regression coefficient. From this point they 
discuss the deviation of observed family means about the 
regression line, breaking it up into two parts. The 
section closes with two approximations for optimum n 
when the sampling variance of the regression coefficient 
has a minimum value. 

In the next section the authors discuss regression of 
offspring on single parent value. Three cases are con- 
sidered and approximations for m are given in each. 
The three cases are: 


(i) each offspring being related only through the 
measured parent. The offspring of each measured 
parent form a half-sib group. 


(ii) each measured parent has only one mate, so that 
the offspring are now full sibs. 


continued 


4,3 (3.5) 
continued 


Biometrics (1960) 16, 348-353 (2 references, | figure) 


(iii) all measured parents have the same mate. The 
groups of offspring are then full sibs with a 
half-sib family. 


Following these three cases is a situation where the 
family is a full-sib group. This discussion also gives an 
approximation for optimum 7. 

In the closing discussion the authors point out that 
if there is an environmental cause of similarity among 
members of a family group, the heritability estimate 
obtained from analysis of variance will be worthless, 
but not that obtained by the regression method. In 
this case the optimum family size will be reduced. In 
general, the optimum size by the regression method is 
smaller than that for determination of the heritability 
from the analysis of variance. 


(J. J. Bartko) 
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JILEK, M. & LIKAR, O. (Central Res. Inst. of Food Industr., and 
Res. Inst. of Food Industr. Econ., Prague) 


4.5 (3.2) 


Tolerance limits of the normal distribution with known variance and unknown 


mean—ZIn English 


Aust. J. Statist. (1960) 2, 78-83 (9 references, 3 tables) 


This paper is a further development of previous work 
by the same authors, see ‘‘ Coefficients for the deter- 
mination of one-sided tolerance limits of normal dis- 
tribution” [Ann. Inst. Statist. Math., Tokyo (1959) 11, 
45-48: abstracted in this journal No. 1/401, 4.5]. 

If X,, ..., X, is a random sample of n observations 
from a normal distribution with mean yu and variance 
o”, it is often necessary to find regions R,4, Rp for which 

(i) R, contains at least a given proportion y of the 
distribution with given probability p; 

(11) Rg contains a given proportion y-of the distribu- 

tion on the average. 

The paper gives a survey of solutions for the tolerance 
limits of these regions when yw is unknown but o? is 
known. 

Two tables are given from which one-sided and two- 
sided tolerance limits of R, can be found for given 
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values of y, 6, and n. A third table is provided from 
which one-sided and two-sided tolerance limits of R, 
can be found for given values of y, and n. 


(J. Gani) 


4.5 (10.1) 


A one-sided analog of Kolmogoroff’s inequality—Jn English 
Ann. Math. Statist. (1960) 31, 483-487 (4 references) 


The inequality Pr{X¥>e}<6(X*)/[e?+&(X’)] for a 
random variable with zero mean is well known. The 
author gives a Kolmogoroff-type generalisation for the 
situation in which X,, ..., X, are such that &(X,) = 0, 
eee eee Ge) 10.01 = 23,2. 71), Namely: if 


o? = &(X?) and s, = of then 
1 


Pr{max (X, +... +x)2e1 <S,/(e7+5,): 
isn 
Further, given a positive ¢ and numbers O73 Bons o* one 
can construct random variables X,, ..., X, having the 
properties of the above mentioned theorem and such 
that the above inequality becomes an equality. 

As another generalisation, at least ifn = 2, the author 
shows that if ¢, and ¢, are given positive numbers then 


Pr{X,2e, or X,+ X22} 
S {03 +.07(02/005)7}/{03 + (2/e1)} 


where «, = o7+n17; i= 1, 2, and y, = min(e;, &), 
N> = &. The variables X, and X, are assumed to 
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be as above. Again given 07, 03, €; and e, the equality 
in the last expression can always be achieved by a 
suitable choice of X,, X,. Several generalisations of 
this are given for the case n = 2 but with less restriction 
ony Xa x. 

The author notes that under the original hypotheses 


OnW\ ae ae ka tbee variables) nce (hn be ear) 
j=1 

form a martingale. Correspondingly he states a version 

of the first inequality for a separable martingale 

{Y,; t=0} with 6(Y, =0, &(¥7) =07(1): namely for 

positive ¢ and t 


Pr jsp ze So7(t)/e* +07(t). 


tSt 


This inequality is sharp if o7(-) is continuous on the 
right. 


(R. Blumenthal) 
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MARSHALL, A. W. & OLKIN, I. (Stanford Univ., and Michigan State Univ.) 


4.5 (1.5) 


A one-sided inequality of the Tchebysheff type—In English 
Ann. Math. Statist. (1960) 31, 488-491 (2 references) 


The authors generalise the one-sided Tchebysheff in- 
equality by obtaining an upper bound for Pr{x,>1 or... 
or x,21} where the x’s are random variables with zero 
means, 6xj=o07 and 6&x;x;=07p (i#/j). The 
theorem is that if 1—o7t>0 and k=>o?(k—1)(1+72) then 
Pr{x,21 or... or x,=1} is equal to 


ko*{/[1+(k— Vp +07 —07(k—1) — p)]+(k-1) 1 —p}? 


{k+o7[1+(k—-1)p]}? 


The situation in which one generalises the two-sided 
Tchebysheff inequality by considering the absolute 
values of the x’s has been treated by Olkin & Pratt [Ann. 
Math. Statist. (1958) 29, 226-234]. The argument here 
relies on some of their results. 


(R. Blumenthal) 
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(Econometric Inst., Netherlands School of Economics) 


4.8 (6.6) 


The bias and moments of the general k-class estimators of the parameters in 


simultaneous equations—IJn English 


Econometrica (1959) 27, 575-595 (6 references, 1 table) 


In this paper the author studies the small sample 
properties of the so-called general k-class estimators 
of simultaneous equations. He states that little is 
known about these properties, the only results being 
those of the multi-dimensional confidence regions 
[Anderson, T. W. & Rubin, Ann. Math. Statist. (1949) 
20, 46-63] and those by Theil for the coefficients of a 
structural equation [Proc. Kon. Ned. Akad. Wetensch. A 
(1950) 53, 1397-1412]. 

The first section of the paper is devoted to the analysis 
of the bias, to the order of T~', T being the number of 
observations in a column vector; and the moment 
matrix, to the order of T~’, of the general k-class 
estimators of an over-identified or just-identified single 
equation which is part of a system of simultaneous 
equations. 

The second part of this paper states two resultant 
theorems on bias and on the moment matrix. The 
author obtains the bias of the two-stage least-squares 
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estimator together with an unbiased estimator, to the 
order of T~', as a corollary to the theorem on bias 
(No. 1). The two corollaries of the theorem (No. 2) 
on the moment matrix give the moment matrix of the 
two-stage least-squares estimator and the “‘ best ’’ value 
of k. 

Section three of this paper gives the proofs of the 
theorems stated in section two and finally, the author 
gives an example illustrating the selection of the best 
k-value. A table of alternative point estimates of the 
coefficients in Klein’s model I is given which also shows 
their estimated finite sample bias and standard errors. 


(W. R. Buckland) 
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NEYMAN, J. & SCOTT, Elizabeth L. 


(University of California, Berkeley) 


4.1 (2.2) 


Correction for bias introduced by a transformation of variables—Jn English 
Ann. Math. Statist. (1960) 31, 643-655 (9 references) 


This paper is concerned with deducing minimum variance 
unbiased estimates of the effects of experimental treat- 
ments expressed in the original units when the analysis 
is based on “normally”? transformed data. The 
solution is obtained for a broad category of transforming 
functions. 

Under the basic assumption that the transformation 
used in the analysis of an experiment is faultless, so that 
the transformed variables exactly follow normal distri- 
butions with some postulated means and with the same 
unknown variance, an expression is derived for minimum 
variance unbiased estimates for the treatment effects 
which is valid for a large class of transforming functions. 
The resulting expression is complex. 

A simpler expression is presented for a sub-class of 
this class of transformations. Included in this sub-class 
are the square-root transformation, the logarithmic 
transformation, the angular transformation and the 
hyperbolic sine transformation. The estimates and the 
bias are specifically discussed for these four trans- 
formations. 
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(University of Leeds) 


The authors mention that formulas for correcting the 
bias introduced by the transformation of variables in 
certain particular cases are contained in some of the 
references but that these formulas do not agree with 
theirs. 


(J. W. Wilkinson) 


4.5 (2.7) 


Tables on the bounds of the probability integral when the first four moments 


are given—IJn English 


Biometrika (1960) 47, 399-410 (4 references, 1 table, 1 figure) 


There is a long history of this problem, given the moments 
of a distribution to specify as narrowly as possible the 
probability integral of that distribution. The authors 
are concerned with the case in which four moments are 
given. Tchebycheff’s inequality is commonly applied 
to the case where only the mean and standard deviations 
are available, and the authors note the improvement in 
the bounds for the probability integral obtained by the 
introduction of the third and fourth moments, The case 
where there is some restriction of the range is also 
discussed. 

Using the usual notation, A= f,—f,—1 (the 
boundary of Pearson’s impossible area), the table gives 
the upper and lower bounds to the probability that a 
random variable is less than a quantity « given A, f, 
and f,. The numerical values considered are « = | 10], 
(1)...| 3 |, | 2°5| (0-5)...0:0; for A = 0-0, 0-5, 1(1), ...4, 
Gaseioe loco and. pi. =.0,15.2,.3):4, 8 and various 
values of f,. 


(Florence N. David) 
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WALLACE, D. L. (University of Chicago) 
Conditional confidence level properties—In English 


4.4 (4.1) 


Ann. Math. Statist. (1959) 30, 864-876 (10 references) 


The confidence level (function) of a region-valued 
estimation procedure is the probability that the region 
covers the true parameter and will in general depend 
on the true parameter value. For any subset of the 
sample space, define a conditional confidence level 
(function) as the conditional probability, given that the 
sample lies in the subset, that the region covers the true 
parameter. 

A variety of properties of conditiona! confidence level 
behaviour of region-valued estimation procedures are 
defined. Of these, the most important are the properties 
S,(a) and S,(a). A procedure is S,(«) if there is no 
subset for which the conditional confidence level is 
always (that is for all parameter values) less than a 
or always greater than «. A procedure is $,(«) if there 
is no subset for which the conditional confidence level 
is always less and bounded away from a. These 
properties are irrelevant in general, but are desirable 
properties for a confidence procedure with confidence 
level «. Fisher’s criticism [J. R. Statist. Soc. B (1956) 
18, 56-60] of the Welch solution of the Behrens-Fisher 
problem is equivalent to showing that the level « Welch 
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WALLACE, D. L. (University of Chicago) 
Conditional confidence level properties—Jn English 


procedure does not have property S,(«). Discussion 
of the desirability of these properties, called by different 
names, and examples of their violation have been given 
by Buehler [Ann. Math. Statist. (1959) 30, 845-864]. 
The present paper obtains and illustrates sufficient 
conditions for a procedure to have the properties S,(«) 
or (a). 

A region-valued estimation procedure is said to be 
“level a Bayes’ if it can be generated as follows. For 
some prior distribution on the parameter space, and 
for each sample, the posterior distribution is obtained 
by Bayes theorem. For each sample, a region in the 
parameter space is chosen having posterior probability «. 
A procedure is said to be “‘ level « weak Bayes ”’ if it is 
generated as above, except that in place of a valid prior 
distribution, a non-negative, but not finitely integrable 
function is used formally in Bayes theorem, with the 
requirement that the posterior distribution obtained be, 
for each sample, a valid distribution. 

Principal results are that any “‘level « Bayes”’ pro- 
cedure has property S,(«) and any “‘ level « weak Bayes ”’ 
procedure has property S,(«). Level « Bayes or weak 


continued 


4.4 (4.1) 
continued 


Ann. Math. Statist. (1959) 30, 864-876 (10 references) 


Bayes procedures will generally not have constant 
confidence level «. But if a procedure with confidence 
level « is simultaneously (say) a level « weak Bayes 
procedure for some prior “ distribution’’, then the 
confidence procedure has conditional behaviour no 
worse than the S,(«) property. The weak Bayes result 
is used to show that the property S,(«), but not generally 
S,(«), holds for the usual normal theory confidence 
procedures based on ¢, x’ and F, and also the Pitman 
confidence procedure for location and scale parameters. 
The Behrens-Fisher fiducial procedure has S;(«), but the 
Welch solution does not. 

A property weaker than S,(«) and some associated 
results are given appropriate for such discrete confidence 
problems as those for binomial and Poisson parameters. 


(D. L. Wallace) 


2/304 


e Me 


wore 


.; 


SWAMY, P. S. (University of Poona) 4,3 (2.8) 


Estimating the mean and variance of a normal distribution from singly and doubly 
truncated samples of grouped observations—In English 


Bull. Calcutta Statist. Ass. (1960) 9, 145-156 (9 references) 


In this paper the author derives the maximum likelihood 
equations to estimate the mean and variance of a normal 
distribution from singly and doubly truncated samples 
of grouped observations with known truncation points, 
when the number of unmeasured observations is 

(i) unknown, 

(ii) known separately for each “‘ truncated tail ’’, 

(iii) known jointly for the two “‘ truncated tails ’’. 
These three cases were considered by B. Raja Rao 
[Proc. Nat. Inst. Sci. India A (1958) 24, 366-376] for 
discussing the relative efficiencies of best asymptotic 
normal estimates. The maximum likelihood equations 
are shown to be easily solvable with the aid of 
Gjeddebaek’s tables of Z, and Z, functions. 

Asymptotic variance and covariance matrix of the 
estimates is obtained in each of these cases and an 
investigation is made of the loss of information due to 
truncation. A numerical example is given illustrating 
the practical application of the results. 


(B. R. Rao) 
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BENNETT, B. M. (University of Washington, Seattle) 


5.6 (5.1) 


Note on the power function of the Y,, test in genetics—Jn English 
Skand. Aktuartidskr. (1959) 42, 1-5 (4 references, 1 table) 


The event E is supposed to occur k times in a sequence 
of n independent trials. Let X,, be the rank total 


XxX, a Vat 2yos.. ny, 


where y; = 1 if E occurs, and y; = 0 if E does not 
occur on the ith trial. Haldane & Smith proposed X, 
as a criterion for testing the randomness of the k occur- 
rences of E, more specifically for testing birth-order 
effect [see Ann. Eug. (1948) 14, 117-124]. _X,, is identical 
with the statistic in Wilcoxon’s test, and related in a 
simple way to the statistic used in the Mann-Whitney 
test. 

In a previous paper [Skand. Aktuartidskr. (1956) 39, 
11-18] the author obtained the conditional distribution 
of X, given k as a function of p,, Po, ...; Pas Where p; 
is the probability that E occurs in the ith trial. In 
particular the case 


log [p;/(1—p,)] = t+ G-1)A 
_ was considered. 
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BERTRAM, G. (Techn. Hochsch., Hannover) 


In the present paper the power of the one-sided test 
(with rejection region Y,<X,) is tabulated for selected 
values of k, n, a, and A; the A-values ranging from’ 
—1:0 to —0-05. 

The author also discusses the modifications which 
the distribution of X, undergoes when the number of 
trials is made dependent on the outcome of the trials. 
If for example £ is the birth of a child affected with some 
disease, the appearance of one or more E might cause 
the parents to cease their family. The author gives 
formulae for the distribution of X,, under the following ~ 
two rules for terminating the sequence: 


(i) Terminate if two successive E’s appear, 
(ii) Terminate if E occurs in the first trial. 


The effect of (i) and (ii) on the distribution is also 
illustrated by values of the mean and variance’ of X, 
inthe case? = 6, k = 21-0, 


(B. Matérn) 


5.7 (5.3) 


Sequential analysis for two binomial experiments—In German 
Zeit. angew. Math. Mech. (1960) 40, 185-189 (4 references, 3 tables, 6 figures) 


The standard method in sequential analysis is to 
construct first a sampling plan (a parallel strip in the 
plane). According to the outcome of experiments, 
a certain step-function is sequentially constructed until 
it reaches one of the boundaries. 

In the present note the author reverses this procedure. 
He first draws a step-function according to his observa- 
tions and with fixed sample size. He then determines 
a sampling plan in such a way that the step-function 
ends at one of the boundaries. He finally claims without 
proving, that the inference which can be drawn by a 
correct sequential procedure, will also hold in this 
particular case. 


(W. Vogel) 
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BHATE, D. H. (Delhi School of Economics) 


5.7 (2.5) 


A note on the distribution of successive sums of samples from an exponential 


population—I/n English 


Bull. Calcutta Statist. Ass. (1958) 8, 13-19 (4 references, 3 tables) 


In this note the author solves the following problem: 
let x,, x2, ... be random observations from an exponential 
population. Let X,= x, X;: it is required to obtain 
the distribution of X, given that a;<X;<b,, i = 1, 2, 
..., '—1, a; and b; being linear functions of i with equal 
gradients. 

The author then applies the result to compute the 
probability of arriving at a decision at the rth stage 
when testing sequentially a simple hypothesis regarding 
the parameter of an exponential population against a 
simple alternative. A sequential test of a simple hypo- 
thesis regarding the variance of a normal population 
against a simple alternative is a special case of this 
problem when the mean of the population is known. 


(S. John) 
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COX, D. R. (Birkbeck College, London) 


5.7 (5.1) 


A note on tests of homogeneity applied after sequential sampling—In English 
J. R. Statist. Soc. B (1960) 22, 368-371 (5 references, 1 table) 


Suppose that in order to examine a parameter 0, observ- 
ations (x,, ..., X,) are obtained, the final sample size n 
being determined by a sequential sampling rule, for 
example by a likelihood ratio test with fixed limits. It 
may now be required to use the observations (x,, ..., X,) 
for a different purpose, for example to test a subsidiary 
hypothesis H, that has tacitly been made in applying 
the sequential procedure. 

In general the use of a standard fixed-sample-size 
test is not strictly justified because of the dependence 
of n on the values (x,, ..., X,).. The object of the present 
paper is to point out that when the subsidiary hypothesis 
H, concerns the homogeneity of the observations the 
standard tests are often justified through an appeal to 
permutation distributions. 

Several examples are given. One concerns the 
situation where @ is a binomial probability, and the 
subsidiary hypothesis is that the value of @ for one 
group of observations is the same as for a second group; 
the two groups might, for instance, be males and 
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females. The standard test for the difference between 
two binomial proportions can be used, despite the 
dependence of n on (x, ..., X,), provided that the entry 
of males and females into the trial is random. 


(D. R. Cox) 


DWASS, M. (Northwestern University, Evanston, Illinois) 


5.6 (3.7) 


The distribution of a generalised D,* statistic—In English 
Ann. Math. Statist. (1959) 30, 1024-1028 (6 references) 


Let F,(x) be the empirical cumulative distribution 
function of n independent random variables, each 
_ distributed according to the same continuous cumulative 
distribution function F(x). The paper is concerned 
with the distribution of the random variable, 


Dy) = sup {F,(x) — yF(x)}. 


Elementary considerations show that this distribution 
does not depend on the form of the continuous distribu- 
tion F(x); hence F(x) is taken to be the uniform distribu- 
tion function for the interval (0, 1). 

When » = 1, the variable D,'(1) is the usual one-sided 
goodness-of-fit statistic whose asymptotic distribution 
was first derived by Smirnoff [Usephi Matem. Nauk, 
(1944) 10, 179-206]. The paper gives two expressions 
for 


Pr(D,'(y)<a) = Pr(F,(x)Sat yx, 0S x1). 


The first formula for the case y = 1 agrees with one 
obtained by Birnbaum & Tingey [Ann. Math. Statist. 
(1951) 22, 592-596]. The second formula, the author 
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FISZ, M. (University of Warsaw, Poland) 


feels, would involve fewer computations for numerical 
evaluation. The second formula lends itself to deriving 
asymptotic results as noo and yields some previously 
published facts concerning Poisson processes. 


(D. A. S. Fraser) 


5.3 (5.6) 


On a result by M. Rosenblatt concerning the Von Mises-Smirnoff test—Jn English 
Ann. Math. Statist. (1960) 31, 427-429 (4 references) 


The author points out and corrects an error involving 
the following integral formula. Let S,(¢) and S,(¢) 
denote the empirical distribution functions corresponding 
to two independent random samples x, (k = 1,...,n) and 
y,(j=1, ..., m) from two populations with the same 
continuous distribution function F(t). It is shown that 


| [S,()-S,(0)/?d[(nS ,() + mS,(1))/n+m-—t] 
0 
a [. [S,(t)—t]’d[S,()-1] 
(0) 


ts | * [8,()—tPa[S,() — 1] = 1/6nm. 


In the original paper, the right-hand-side of the 
above equation was claimed to be zero. However, the 
assertion of Rosenblatt’s theorem remains true since he 
was concerned with n, m->oo. Rosenblatt’s original 
paper ‘“‘ Limit theorems associated with variants of the 
von Mises statistic’? was published in Amn. Math. 
Statist. (1952) 23, 617-623. 


(P. L. Meyer) 
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HALPERIN, M. (National Institutes of Health, Bethesda, Maryland) 


5.6 (5.3). 


Extension of the Wilcoxon-Mann-Whitney test to samples censored at the 


same fixed point—In English 


J. Amer. Statist. Ass. (1960) 55, 125-138 (9 references, 15 tables) 


A non-parametric two-sample test for samples censored 
at the same fixed point is obtained by extension of the 
Wilcoxon-Mann-Whitney test. The test is proposed 
for samples of size n on F and size m on G, where F 
and. G are continuous cumulative distributions. The 
observations on F and G are denoted by x and y 
respectively and the point of censoring is denoted T, 
with all x and y in excess of T being censored; there 
being r, censored of the former and r,, censored of the 
latter, or r = r,+r,, censored in all. The null hypothesis 
tested is F(x) = G(x), except possibly in the censored 
portions. 

The statistic employed is denoted by U, and is 
equivalent to the Wilcoxon-Mann-Whitney U statistic 
with an adjustment which treats all censored observations 
as ties. It is shown that under the null hypothesis the 
conditional distribution of U, given r is independent 
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HALPERIN, M. 


same fixed point—IJn English 


(National Institutes of Health, Bethesda, Maryland) 
Extension of the Wilcoxon-Mann-Whitney test to samples censored at the 


of F. This conditional distribution permits the com- 
putation of the exact probability distribution of the 
test statistic. 

The author employs this result to compute lower 
five per cent. and one per cent. points of U, for m and 
n less than or equal to eight and for r less than or equal 
to m+n—1. As a consequence of discreteness the 
significance points are not exact. However, the exact 
probabilities of the nominal five per cent. and one per 
cent. points are given parenthetically; and seven tables 
each are given for these percentage points with tables 
ordered by n and entered marginally by m and r, the 
entries being U,, the integer value of U, nearest the 
appropriate percentage point, and its corresponding 
exact probability level. 

The exact mean and variance of U, under the null 
hypothesis are given and it is shown that the standardised 


continued 


5.6 (5.3) 
continued 


J. Amer. Statist. Ass. (1960) 55, 125-138 (9 references, 15 tables) 


U. is asymptotically normal with zero mean and unit 
variance. Thus a large-sample test is available for 
values of n, m and r outside the tables. A tabular 
comparison of exact and asymptotic probabilities is 
given for n = m = 8 which suggests that for no more 
than 75 per cent. censored observations the asymptotic 
theory is adequate. 

The consistency of the U, test is shown, after the 
proof given by Mann and Whitney for the U test, with 
special consideration for the dependence of the (1+m) 
variables in the censored case. 


(S. Krane) 
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JILEK, M. & LIKAR, O. (Cent. Res. Inst. Food Industry and 


5.6 (3.8) 


Res. Inst. Food. Industr. Econ., Prague) 
Table of random sample sizes needed for obtaining non-parametric tolerance 


regions—In English 
Zastosowania Mat. (1960) 5, 155-160 (1 table) 


A known one-dimensional result of Wilks about order 
statistics is as follows: the probability P that the interval 
between the rth smallest and the sth largest value in a 
random sample of size n taken from a population with 
continuous distribution covers at least the proportion p 
of the population distribution is given by 
P=I1,_,m, n—m+1), 

where [,(a, b) is Pearson’s incomplete Beta function 
andm=r+s. 

This result and its multi-dimensional generalisations, 
showing the applicability of the above-mentioned 
relation in the multi-dimensional case as well as in the 


one-dimensional are the basis for construction of non- 
parametric tolerance regions. 
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KHATRI, C. G. (University of Baroda, India) 


The authors have computed a table based on P = 
I, ,»(m, p—m-+1) which differs from some other known 
tables in that it determines n, if m, p, and P are given. 
The table is computed for m = 1; 2(2)10; p and P = 
0-50, 0-75, 0-90, 0-95, 0-975, 0-99, 0-995, 0-999. 


(S. Zubrzycki) 


5.4 (2.6) 


On testing the equality of parameters in k rectangular populations—In English 
J. Amer. Statist. Ass. (1960) 55, 144-147 (4 references, 2 tables) 


Rider’s statistic for testing equality of ranges of two 
rectangular populations and Murty’s statistic for testing 
equality of ranges of two rectangular populations when 
lower limits are zero are extended to the analogous 
k-sample tests. 

The k-sample test statistic in each case is taken as 
the minimum over all pairs of the corresponding two- 
sample statistic. Thus, for testing equality of ranges 
of k rectangular populations the statistic is upin, the 
minimum of ratios of all pairs of the sample ranges, 
and for testing equality of ranges of k rectangular 
populations when the lower limits are zero the statistic 
iS Umins the minimum of ratios of all pairs of sample 
maximum values. 

Sampling distributions of these statistics are obtained 
under the null hypotheses and the lower five per cent. 
points of the statistics are tabulated for the case in 
which the k samples are of the same size, n. For Umin 
the table has arguments k = 2(1)5; nm = 4(1)10(5)20 
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and four-decimal entries, except for k = 5 in which 
case the entries are given to three decimals. For Unin 
the table has arguments k = 2(1)11; nm = 1(1)10(5)30, 
40, 60, 100, 500, 1000, and four-decimal entries. An 
example is given. 

(Note: a minor erratum is found in section 3(b) 
wherein “ Vmin (OF Umax)” Should be read for “‘ Umin (OF 
Umax) in the first sentence.) 


(S. Krane) 


ROY, J. (Indian Statist. Inst., Calcutta) 


5.8 (11.1) 


Tests for independence and symmetry in multivariate normal populations—In English 
Sankhya (1960) 22, 267-278 (16 references, 3 tables) 


Likelihood ratio criteria for a number of hypotheses 
of interest in multivariate analysis have, under the 
respective null hypotheses, a tth moment of the form 


s 


Il [Gn A | ca [Gn oh b;), | 


where the 5; and c, are constants not involving n (the 
sample size) and satisfying the inequalities c,;>b is 
J = 1,2, ..., 8; (m), =T(@m+0/T(m). Box [Biometrika 
(1949) 36, 317-346], Rao [Bull. Int. Statist. Inst. (1951) 
33, IL 177-180] and Roy [Bull. Int. Statist. Inst. (1951) 


33, II 219-230] provide series expansions in powers of 


(1/n) for distribution functions of such criteria. 
denote the criterion by L and set 


If we 


Me 2 (cj bp, t= 1, 2 m=2r,,¢@=(%1—T)/r;, 


N=n-a, X =Nlog,L 


and 


Q(x) = | 2 UT Gy tat te du, 
Pr(X 2x) = O,(x)+N~*a2[O,+4(x)—Q,(x)] 
correct to terms of order 1/n?. 
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ROY, J. (indian Statist. Inst., Calcutta) 


Tests for independence and symmetry in multivariate normal populations—In English 


The present paper gives tables of r, a and a, by the 
author for the following hypotheses considered in the 
references given against them: 


(1) Specified groups of variables are mutually in- 
dependent [Econometrica (1935) 3, 309]. 


(2) The variates are mutually independent and have 
equal variances [Ann. Math. Statist. (1940) 11, 
204]. 


(3) The variances are equal and the covariances are 
equal [Ann. Math. Statist. (1946) 17, 257]. 


(4) The means are equal, the variances are equal 
and the covariances are equal [Ann. Math. Statist. 
(1946) 17, 257]. 

(5) Specified groups of variables are such that within 
each group the variances are equal and the co- 
variances are equal, and between groups, the 


covariances are equal for each pair of groups 
[Ann. Math. Statist. (1948) 19, 447]. 


(6) Specified groups of variables are such that within 
each group the means are equal, the variances 
are equal and the covariances are equal, and 


continued 


5.8 (11.1) 
continued 


Sankhya (1960) 22, 267-278 (16 references, 3 tables) 


between groups, the covariances are equal for 
each pair of groups [Amn. Math. Statist. (1948) 
19, 447]. 


(7) Specified groups of variables are such that within 
each group the variances are equal and the co- 
variances are equal, and between groups, the 
diagonal covariances are equal and the off-diagonal 
covariances are equal [Ann. Math. Statist. (1948) 
19, 447]. 


(8) Specified groups of variables are such that within 
each group, the means are equal, the variances 
are equal and the covariances are equal, and 
between groups, the diagonal covariances are 
equal and the off-diagonal covariances are equal 
[Ann. Math. Statist. (1948) 19, 447]. 


(S. John) 


2/318 


SAVAGE, I. R. (Univ. of Minnesota) 


5.6 (5.1) 


Contributions to the theory of rank order statistics. The one sample case—In English 
Ann. Math. Statist. (1959) 30, 1018-1023 (4 references) 


The one-sample problem is analysed using techniques 
developed by the author in earlier papers [Ann. Math. 
Statist. (1956) 27, 590-615 and (1957) 28, 968-977]. 
Tests considered for this problem are distribution-free 
and are expressed in terms of the statistic Z = (Z,, ..., 
Zy) having Z; = 1(0) according as the ith smallest value 
in a sample of N is positive (negative). For sampling 
from a distribution with density f(x) an N-dimensional 
integral expression is given for Pr(Z). 

Let a 0 in Z’ be replaced by a 1 yielding Z or let a 
0 and i in Z’ be interchanged (by moving the 1 to the 
right) yielding Z. The author gives conditions under 
which Pr(Z)>Pr(Z’). These conditions are fulfilled by 
the normal distribution. 

The results can be used for finding and for analysing 
tests of the hypothesis that the sample is from a distribu- 
tion symmetrical about zero, against alternatives such 
as slippage to the right. Under the hypothesis each 
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possibility for Z is equally probable. Under alternatives 
inequalities such as Pr(Z)>Pr(Z’) indicate that Z should 
be a part of the rejection region before Z’. Properties 
of the Wilcoxon one-sample signed-rank test that are 
related to these results are considered. 


(D. A. S. Fraser) 


5.6 (3.8) 


Contributions to the theory of rank order statistics. Computation Rules for probabilities 


of rank orders—Jn English 


Ann. Math. Statist. (1960) 31, 519-520 (3 references) 


The author obtained two computational rules for 
evaluating probabilities involving rank orders when the 
null hypothesis is not the true hypothesis. These results 
avoid the usually difficult multiple integrations and are 
obtained for one and two-sample problems. The 
derived rules are ‘‘ back-recursive ’: for example, in 
the two-sample problem, probabilities are computed 
involving samples of size m and n from probabilities 
involving size m+1 and n. These formulae, although 
of only limited analytical value, may be quite helpful 
for Monte Carlo sampling. 


(P. L. Meyer) 
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WEIBULL, M. (Inst. Statistics, University of Lund, Sweden) 


5.6 (3.6) 


Moments of the difference between means in two samples from a finite population: 
applied in connection with a randomisation test—In English 


Skand. Aktuartidskr. (1959) 42, 36-60 (9 references, 4 tables, 2 figures) 


Two samples are chosen at random without replacement 
from a finite population of N individuals. The sizes 
of the samples are n, and ny, respectively. Thus 
ny+n,SN. The author considers the distribution of 


u= (m,—m,)/s /(Q/n, +1/nz) 


where m, and m, are the observed means of the two 
samples and s? is the variance of the finite population. 

The moments and cumulants of u are given in terms 
of the polykays of the finite population. Explicit 
expressions are presented for cumulants up to order 
eight. It is found that uw has a symmetrical distribution 
if the finite population is symmetrical, or if n, = n): 
the author also studies the asymptotic qualities of the 
distribution. The applications of these results concern 
the use of u for testing the significance of the difference 
m,—m,. If ny+n, = N, the test is identical with 
Pitman’s randomisation test for two samples. 

Three types of population are investigated numerically. 
In the first example the population values are 1, 2, ..., 
N; thus w is related to the test-statistic of the Wilcoxon- 
Mann-Whitney test. In the second example the values 
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are 0, 1, 3, 6, ..., N(N—1)/2. For samples from these 
populations the author compares the exact probabilities 
of exceeding certain critical values of u with probabilities 
computed by the normal approximation, and by the 
Gram/Charlier A-series. The approximation by the 
A-series with terms up to and including the fourth order 
is satisfactory for ny, n,=5. The inclusion of terms of 
order five and six does not result in substantially better 
approximations. The normal approximation is not 
satisfactory within the range of n values studied. 

The third finite population has k-statistics coinciding 
with the cumulants of a normal distribution. The 
probability of exceeding critical values is approximated 
by the above-mentioned methods. A comparison is 
made with the corresponding probability in that marginal 
distribution of u which is obtained when the finite 
population is considered as a set of independent observ- 
ations of a normally distributed random variable. 

The author provides a table of approximate critical 
values of u. The table is entered with the coefficients 
of the terms of order three and four in the Gram/ 
Charlier series. 

(B. Matérn) 


5.6 (5.1) 


A generalised Pitman efficiency for nonparametric tests—In English 
Ann. Math. Statist. (1960) 31, 405-414 (7 references, 4 tables) 


To compare two tests of a hypothesis, one may consider 
the different sample sizes necessary for the different 
tests to achieve a specified power against a specified 
alternative. Pitman’s criterion of asymptotic relative 
efficiency, which is at present essentially the only practical 
criterion available for comparing non-parametric tests 
is simply the limit of the ratio of the two sample sizes 
as either the specified power approaches unity, or as 
the specified alternative “‘ approaches” the null hypo- 
thesis. The purpose of this paper is to obtain rates 
of convergence of this ratio to its limit, namely the 
asymptotic relative efficiency, thus allowing comparative 
inferences to be made between tests for smaller sample 
sizes than was previously possible. 

The technique used by the author is the ex- 
pansion of the power functions of the appropriate 
tests by means of Edgeworth’s series, up to terms of 
order O(1/n?), where n is an integer which determines 
the sample sizes. In a one-sample problem, n may be 
taken as the sample size itself. 
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The author works out explicitly the rates of con- 
vergence for the following comparisons: 


(i) For the two-sample problem, the Wilcoxon test 
is compared with both the x test and the ¢ test 
for both one-sided and two-sided normal alter- 
natives. 


(ii) For the one-sample problem; the sign test is 
compared with the X test and the f-test for one- 
sided normal alternatives. Some numerical com- 
putations are also provided. 


(R. Pyke) 


ADAMS, E. “ MESSICK, S. (Univ. California, Berkeley and 


6.4 (-.-) 


Educational Testing Service, Princeton) 
An axiomatic formulation and generalisation of successive intervals scaling—In English 


Psychometrika (1958) 23, 355-368 (27 references) 


A formal set of axioms is presented for the method of 
successive intervals and directly testable consequences 
of the scaling assumptions are derived. Then by a 
systematic modification of basic axioms the scaling 
model is generalised to non-normal stimulus distributions 
of both specified and unspecified form. 

Scaled values obtained in this model, at least under 
certain circumstances, are unique up to a linear trans- 
formation and have two interesting consequences for 
the original successive intervals model based on the 
normality hypothesis. 


(i) If in the errorless case the original model fits, 
then no other successive intervals model which 
assumes a different form for the distribution 
functions will fit. The reason for this is that the 
forms of the distribution functions are determined 
by the values of p, ; lying above the point f,. 
Hence, if the ¢; are determined up to linear 
transformation, so are the curves pg, ;. 
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BARBANCHO, A. G. 
Identification of complete models—IJn Spanish 


Estadistica Espafiola (1960) 6, 5-25 (8 references) 


The author considers the problem of identification of 
relations in complete models; that is the possibility 
of distinction among relations in determined systems of 
equations, from the statistician’s point of view. He gives 
special attention to the so-called multi-equation lineal 
stochastic models. Examples are given to show analyti- 
cally and graphically whether or not the system is 
identifiable. 

Stating the problem in a general form the author 
considers models with random disturbances; classified 
in identifiable, unidentifiable and not exactly identifiable, 
corresponding to the usual terms applying to structural 
models: identifiable, underidentifiable, and overidentifi- 
able. Necessary and sufficient conditions for identifying 
coefficients in structural relations are given; as well as 
practical rules for deciding if a structural equation is 
identifiable on the basis of the given system itself. 


(F. Azorin) 
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(Instituto Nacional de Estadistica, Madrid) 


(ii) Where the normality assumption does not fit the 
data it is theoretically possible to use the present 
generalisation to obtain the scale. Then the 
deviation in the scale values from those obtained 
under a normality requirement can be evaluated. 
This, at least in principle, provides a second kind 
of goodness-of-fit besides the usual least-squares 
regression methods employed where the data 
do not exactly fit the Thurstone model. 


The present axiomatic characterisation of a well- 
established scaling model was attempted because of 
certain advantages which might accrue: Firstly, an 
ease of generalisation that follows from a precise know- 
ledge of formal properties, secondly an ease of making 
comparisons between the properties of different models. 
An important outcome of their analysis is that the 
assumption of normality has directly verifiable con- 
sequences and should not be characterised as an un- 
testable supposition. 


(R. E. Stoltz) 
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BROWN, T. M. (Ottawa, Canada) 


Simplified full maximum likelihood and comparative structural estimates—JIn English 
Econometrica (1960) 27, 638-653 (10 references, 2 tables) 


6.6 (4.8) 


In this paper the author, after commenting on the much 
greater interest now taken both in economic theory and, 
for economic policy considerations, of macroeconomic 
models, states that with the present knowledge the full 
information maximum likelihood estimates with un- 
restricted correlations among contemporaneous disturb- 
ances is the best theoretical method of estimation. It is 
essentially a large sample method. Even with the 
simplification afforded by assuming that the disturbances 
vector has a multivariate normal distribution, the com- 
putations are formidable and bulky. Less difficult 
methods have been developed; for example, the limited 
information maximum likelihood method by Anderson, 
T. W. & Rubin, “ Estimation of the parameters of a 
single stochastic difference equation in a complete 
system’ [Ann. Math. Statist. (1949) 20, 46-63] from a 
suggestion by Girshick. 

The main purpose of this paper is to explore a simplifi- 
cation of the mathematics and the computational layout 
for the full information method. After working through 
the basic theory of estimating the parameters of a com- 
plete economic model by this full information method 
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Simplified full maximum likelihood and comparative structural estimates—In English 


a simplified and condensed routine is developed. This 
layout makes for shorter computations because it reduces 
the direct use of restriction matrices and eliminates the 
redundant and rather bulky computations involving a 
kronecker product. The routine is then applied step- 
by-step to the estimation of the parameters of a known 
numerical model (Model X) which has generated its 
own data. 

Various other methods of estimation are less costly; 
viz., least-squares regression, limited information single 
equations, limited information subsystem and diagonal 
covariance. Some of these methods are tested both on 
the above synthetic data and on actual data of a small 
Canadian model (Model Y) in order to obtain a pre- 
liminary assessment of the efficiency of the various 
methods. Models X and Y are estimated by each 
method in turn and tables are given to illustrate the 
results. The author uses the word “efficiency ”’ to 
mean the inverse function of both bias and computational 
costs: bias is considered relative to the size of the 
corresponding parameter. In order “to appraise the 
relative values of the various methods, computational 


continued 


6.6 (4.8) 
continued 


Econometrica (1960) 27, 638-653 (10 references, 2 tables) 


costs are weighted against the estimates of bias, these 
weights will change with the needs of a particular 
problem and available resources’’. In his conclusion 
the author hopes that continued empirical research will 
provide even better estimates of the efficiency of the 
various methods. He says that until more definite 
results become available, the present results suggest 
that a good research procedure is to use least-squares or 
limited information programmes for probing alter- 
native hypotheses and followed by full information 
maximum likelihood or truncated full maximum likeli- 
hood programmes for making final estimations. 


(W. R. Buckland) 
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6.6 (0.6) 


The computations of generalised classical estimates of coefficients in a 


structural equation—In English 


Econometrica (1959) 27, 72-81 (11 references, 7 tables) 


This paper illustrates the computation of some general- 
ised classical linear estimates of coefficients in a structural 
equation. The author states that he has presented his 
facts in a concise manner suitable for comprehension 
by the non-statistical reader and that his present paper 
is a natural sequel to an earlier one, ‘‘ A generalised 
classical method of linear estimation of coefficients in a 
structural equation’ [Econometrica (i957) 25, 77-83]. 
This present paper is also supplementary to the paper by 
Girshick & Haavelmo, “ Statistical analysis of the 
demand for food: examples of simultaneous estimation 
of structural equations ” [Econometrica (1947) 15, 79- 
110]. This latter paper deals with the computation of 
limited information single equation estimators. 

For the purpose of comparing the limited information 
single equation and the generalised classical estimation 
techniques, where these lead to different results, the 
author applies the generalised classical method to the 
three over-identified structural equations of the Girshick 
Haavelmo model. 
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CURETON, E. E. (University of Tennessee) 
Note on $/@jax—/n English 
Psychometrika (1959) 24, 89-91 (2 references) 


Using conventional methods for computing the fourfold 
correlation coefficient, ¢, the author first points out a 
serious limitation on the use of that coefficient. Since 
¢ is a correlation coefficient, it should be free to vary 
within the limits of +1. These limits, however, can only 
be reached using the ¢ coefficient when certain conditions 
are met, that is to say, when marginal proportions are 
0-5. 

The author proposes a solution to the problem which 
enables the limits of +1 to be always attainable. This 
is accomplished by dividing ¢ by $max, Pmax being the 
maximum value of @¢ consistent with given marginal 
proportions. A formula is given for use when @ is 
positive, and an alternative is given when ¢ is negative. 
In addition, formulas which allow easier computation 
and data handling are also given. 


(D. Schum) 
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The first part of the paper is devoted to an explanation 
of the composition of the Girshick-Haavelmo model, 
a table showing the data used in the example and in- 
structions for the preliminary computations. Section 
four describes the computation of a single structure 
equation: the author states that he is using the forward 
Doolittle method [see Klein, A text book of Econometrics 
(1953) pp. 151-165. Evanston: Row, Peterson] in order 
to carry out the computations. 

In conclusion the author points out that the com- 
putations which he has outlined constitute the essential 
steps in the generalised classical estimation procedure 
and he has used numerical illustrations to illustrate the 
practical applications of this procedure. He makes some 
comparisons between the results of limited information 
single equation and generalised classical estimation 
procedures and discusses at some length how such 
comparisons should be made. The need is stressed for 
further computational research if accurate comparisons 
are to be made. 


(W. R. Buckland) 


6.9 (6.2) 


DUDZINSKI, M. L. (Division of Math. Statist., CSIRO, Canberra) 


6.1 (4.3) 


Estimation of regression slope from tail regions with special reference to the volume 


line—IJn English 


Biometrics (1960) 16, 399-407 (3 references, 4 tables) 


This paper is an extension of the approximate method 
given by M. S. Bartlett [Biometrics (1949) 5, 207-212] 
for estimating the regression line slope from the join 
of means of tail regions for evenly spaced values of the 
predicting variable. The author applies this general 
study to the fitting of a regression line to commercial 
volume on sectional area at breast height in stands of 
even-aged trees. The problem is explained in detail 
with statements justifying the assumptions made. Array 
variances were assumed to increase with an increase in 
the dependent variable. The author proceeds to show 
that the tail-region method for slope estimation is quite 
satisfactory under these conditions. 

Other investigations are cited in which the tail region 
method was used, assuming both constant and changing 
array variances. Tables demonstrating optimum tail 
region sizes and efficiencies for two f-distributions are 
given. 

The author presents data on measurements from 
individual trees on three separate plots. The initial 
task included a least-squares fit to each set of data, 
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omitting the concept of changing array variances: later 
a close observation was made on the squares of deviates 
of actual values from the estimates. This was done 
in order that a general functional relationship could 
be derived determining array variance as a function of 
diameter. No such generalisation was found. Following 
a description of the method for computing weights in 
order that a weighted regression line can be fitted to 
the data, tables concerning the goodness of the weight 
function are given. At this point the slope is estimated, 
using both weighted and unweighted regression, these 
results are compared with the less laborious task of 
estimating by joining tail region means. Tables showing 
results of calculations for all three tree sites are given, 
these results including efficiencies of the slope estimators 
for each of the three methods. 

In concluding, the author emphasises the importance 
of studying the possible heteroscedasticity property in a 
regression study: if this is not taken into consideration, 
the result will produce an unbiased but less efficient 
estimate of the slope. 


(R. H. Myers) 


6.8 (4.3) 


Polychotomous quantal response in biological assay—In English 
Biometrics (1960) 16, 382-398 (12 references, 3 tables) 


A quantal response in biological assay has usually two 
possible outcomes. In a bio-assay, such as that based 
on mortality of the housefly, there are however, three 
possible outcomes; namely alive, moribund and dead. 
More generally in other types of assays, many outcomes 
might be possible; for example varying degrees of 
moribundity besides alive and dead. If a response is 
polychotomous it is more efficient to use this information 
explicitly in analysing the data rather than pool certain 
outcomes to make the response dichotomous. 

In this paper the authors present general results for 
a trichotomous response assay, in which the estimates 
of the parameters are obtained by minimising appro- 
priate chi-square expressions. These chi-square expres- 
sions are generalisations of y? (Normits) and x” (Logits) 
of Berkson [Biometrika (1957) 44, 411-435; J. Amer. 
Statist. Ass. (1955) 50, 529-550] and have the same 
asymptotic properties of consistency and efficiency. 
The results are presented in the context of insecticidal 
assays based on mortality and moribundity of the 
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housefly. A general procedure is also described for 
dealing with any number of outcomes; for example, 
different degrees of moribundity and a general extension 
of Berkson’s minimum Normit and Logit chi-square 
technique is indicated in an appendix. 

The authors also consider the problem of a trichotom- 
ous parallel line assay and include as an illustration the 
analysis of an experiment on Gunthion residue. 


(H. O. Posten) 


< 
‘ 
=~ 
‘ ' 
“ ‘ 
, 7 
; ‘ 
y ’ ‘0 ri 
; 
by Das 
re ; 
" r 1 , 
I 
' 


HOOPER, J. W. (RAND Corp., Santa Monica, California) 
Simultaneous equations and canonical correlation theory—In English 


Econometrica (1959) 27, 245-256 (8 references) 


In this paper the author has a double purpose; firstly, 
to contribute towards the understanding of simultaneous 
equations estimation as a specific part of the more 
general theory of multivariate analysis, and secondly, 
he aims to show by analogy to single-equation systems: 
“in single-equation systems there are regression co- 
efficients which measure systematic relationships between 
the variable to be explained and the explanatory 
variables; and there are also correlation coefficients 
which measure the extent to which these systematic 
relationships do, in fact, explain the fluctuations in the 
variable to be explained.” 

“Regarding simultaneous equations, only the former 
aspect has been developed and it is the second object 
of this paper to develop a generalised correlation co- 
efficient which measures the extent to which the system- 
atic relationships explain the fluctuations in the set of 
ali jointly dependent variables.” 

In the first part of the paper the author uses a topic 
from multivariate analysis, canonical correlation theory, 
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KARON, B. P. & ALEXANDER, I. E. 


(Princeton University, N.J.) 


6.6 (6.2) 


in the development of a generalised correlation co- 
efficient for simultaneous equation systems. He first 
presents canonical correlation theory within the frame- 
work of simultaneous equations and later the generalised 
correlation coefficient is developed from these results. 
The asymptotic variances of this statistic are derived 
and some special problems in dealing with particular 
relation to econometric equation systems are discussed. 
An example is given in which the results obtained in 
the previous section of the paper are applied to an 
econometric model: Tintner’s model of the American 
meat market [Econometrics (1952) New York: Wiley]. 


(W. R. Buckland) 


6.9 (6.5) 


A modification of Kendall’s tau for measuring association in contingency tables—Jn English 


Psychometrika (1958) 23, 379-383 (3 references) 


A coefficient of association t’ is described for a con- 
tingency table containing data classified into two sets 
of ordered categories. Within each of the two sets the 
number of categories or the number of cases in each 
category need not be the same. +t’ = +1 for perfect 
positive association and has an expectation of zero for 
chance association. In many cases t’ also has —1 as a 
lower limit. The limitations of Kendall’s t, and t, and 
Stuart’s t, are discussed, as is the identity of these co- 
efficients to t’ under certain conditions. Computational 
procedure for 1’ is given. 

A test of significance for the hypothesis of no associ- 
ation exists, since, for all t statistics, the test of significance 
is based on the distribution of S and not of Tt. 


(R. E. Stoltz) 
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KASTENBAUM, M. A. (Oak Ridge Nat. Lab., Tennessee) 
A note on the additive partitioning of chi-square in contingency tables—In English 


6.9 (5.2) 


Biometrics (1960) 16, 416-422 (6 references, 2 tables) 


This paper is concerned with the partitioning of the 
asymptotic chi-square statistic involved in the usual 
test of independence of rows and columns in a rxs 
contingency table. The author presents formulae for 
the partition into one-degree-of-freedom sums of 
Squares representing sets of (r—1) (s—1) orthogonal, 
interaction-type contrasts on the cells of the table for 
a broad class of such contrasts. 

The author notes that the overall test of independence 
is equivalent to a test of the hypothesis that s multi- 
nomial samples of r-cells have been drawn from the 
same multinomial population and that in that case 
one would be interested in partitioning the total statistic 
into components which would allow for tests on subsets 
of columns (samples). The formulae are given for the 
partition into three additive components, one represent- 
ing non-homogeneity of the first 7 columns with (r—1) 
(m—1) degrees of freedom, one representing non- 
homogeneity of the last s—m columns with (r—1) 
(s—m—1) degrees of freedom, and one representing 
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the contrast of the first m columns with the last s—m 
columns with r—1 degrees of freedom. The author 
notes that these three components do not necessarily 
provide indepedent chi-square tests, since real non- 
homogeneity as measured by one component may affect 
the distribution and interpretation of another com- 
ponent. 

An example is given of the analysis when the non- 
homogeneity of row-samples is the variation of interest. 


(C. C. Thigpen) 


631-2) 


A synthesis of two factor analyses of intermediate algebra—In English 
Psychometrika (1959) 24, 343-359 (14 references, 18 tables) 


A battery of 18 tests of intermediate algebra and 20 
reference tests was administered to two successive second- 
year algebra classes. Each battery was separately 
factor analysed by Thurstone methods. The two analyses 
were synthesised by the Tucker method. The five 
congruent factors obtained were identified as: verbal 
comprehension, deductive reasoning, algebraic manipul- 
ative skill, number ability, and adaptability to a new 
task. 

The study raised some interesting questions. First, 
what serious consideration can be given to the factor 
loading in a single factor analysis, by the Thurstone 
multiple-factor technique, if in two studies of almost 
identical nature only five of the twelve factors obtained 
were congruent. The author raises the question of 
whether it would be wise in the future to conduct all 
factor analysis studies in two parts, to analyse each part 
separately, and then to test the two parts for congruence. 

The author suggests the possibility that Tucker’s 
technique crams a maximum amount of congruence into 
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the first few factors, and thereby creates factors that are 
psychologically complex and somewhat obscure. He 
feels that this is added to the difficulty of identifying 
the factors with any large degree of assurance, and, too, 
it is added to the resulting high coefficience of correlation 
between the factors after rotation. 


(R. E. Stoltz) 
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LORD, F. M. (Educational Testing Service, Princeton, N.J.) 


6.4 (2.1) 


The joint cumulants of true values and errors of measurement—Jn English 


Ann. Math. Statist. (1959) 30, 1000-1004 


Suppose U fallible measurements Ni ose p Ol. a Char 
acteristic are made on each of a large number of objects. 
These measurements and the true value (="¢)are 
considered random variables with errors, e, = x,—€, 
u=l1,..., U; the errors e and € are termed latent 
variables. 

The cumulants of the U’s are denoted by "Ke; |. ce 
where C, is a non-negative integer referring to the 
variable x,, and the cumulants of the latent variables 
are Kg,,..,By Where the subscript By refers to € and 
the others to the U measurement errors. 

The article gives explicit expressions for all latent- 
variable cumulants K through the fourth order in 
relation to the observed-variable cumulants « for the 
case U = 4. Formulas are derived for determining any 


U 
Pth order cumulant (?- y B.) of the observed 


u=0 
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measurements as a linear function of the Pth order 
cumulants of the latent variables. The case where each 
error of measurement has zero mean and is uncorrelated 
with every product of the remaining latent variables is 
also discussed. 


(B. M. Bennett) 


6.2 (9.2) 


Analysis of variability of statistical relations—Jn Italian 
R.C. Ist. Lombardo Sci. Lett. (1958) 92, 621-630 (2 references) 


The correlation between two normal variables (linearly 
dependent one on another) is supposed to be under 
the influence of associated conditions. These can be 
considered as particular quantitative or qualitative 
determinations of a third variable or characteristic, or 
of a combination of variables and/or characteristics 
associated with the two principal random variables. 

Corresponding to each particular concomitant con- 
dition, Bernoulli samples are taken from the two- 
dimensional normal variable; admitting variations in 
sample size from case to case. In order to ascertain 
whether the variations of the sample statistics r, with 
changing conditions, are significant, and to estimate their 
direction, recourse is made to Fisher’s z transformation 
so as to have an approximately normal sample dis- 
tribution. 

When the concomitant conditions are particular 
attributes of a qualitative characteristics, comparing 
different values of z with a Student’s f-test one can 
ascertain possible significant variations in r. When 
some of those conditions can be considered a particular 
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determination of a discrete variate, the use of the 
generalised ¢ permits the application of the criteria of 
analysis of mean values, introduced by Pompilj and 
to divide the variations of the z’s observed in orthogonal 
comparisons into simple components, each of which 
shows the influence on correlation owing to elementary 
variations of the concomitant conditions. 

The criteria discussed are illustrated through an 
application to data by Eden & Fisher [“‘ Studies in Crop 
Variation, IV,” J. Agric. Sci. (1927) 17, 548-562] where 
the two main normal variables are the weight of the 
seed and that of the straw of oats sowed in 96 plots 
in eight randomised blocks, while the concomitant 
associated conditions are date of fertilising, quality and 
quantity of fertiliser. 


(A. Naddeo) 
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PAWLOWSKI, Z. (Central School of Planning and Statistics, Warsaw) 


6.1 (4.1) 


On the total differential method and its efficiency in the case of linear regression—In English 


Zastosowania Mat. (1960) 5, 97-118 


The author considers a series of random variables Ye 
depending on a discrete time t. Moreover, it is assumed 
that the expected value of Y, is a linear function of A 
real variables x,, ..., x,, whose values at each time 
moment may be observed. The co-efficients of that 
linear function are to be estimated. A further assump- 
tion is, that the residuals Z,, that is the differences 
between Y, and the linear function in the variables, 
form a stationary time-series with zero mean and finite 
variance. 

The author constructs the estimators of the co- 
efficients in question with the aid of least-squares 
method; but applied to the differences Y,,,—Y;,, 
Z,+1—Z,, and of all the other variables involved, instead 
of applying it to the raw values of Y,, Z,, etc., as 
is usually done. Hence the term “total differential 
method ”’. 

After this general introduction it is proved for h = 1 
that the estimators so obtained are unbiased and, under 
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some additional assumptions, consistent. The efficiency 
of total differential estimators and of classical ones is 
compared. The paper ends with a preliminary discussion 
of the possibility of applying the total differential 
method in the case of non-linear regression. 


(S. Zubrzycki) 


6.9 (5.1) 


A note on the use of Mood’s likelihood ratio tests for item analysis involving 2 x 2 tables 


with small samples—IJn English 
Psychometrika (1959) 24, 371-372 (2 references) 


Mood’s likelihood ratio test is generally considered an 
unreliable y* approximation in 2 x2 contingency tables 
containing expected cell frequencies less than five. 
Probability values were computed for 60 such tables as 
part of an item analysis for two thirty-item alternative 
forms of a measure. The rank orders of the items, 
from the best to the worst differentiators as determined 
separately by Mood’s tests and by Fisher’s exact test, 
correlated 0-97 for one form and 0-96 for the other. 
Thus, although assumptions for its use were not met, 
Mood’s tests give a very good indication of the relative 
probability values for the items. This may indicate 
that where ranking is the goal, as in the item analysis 
described in the article, the y? approximation, by Mood’s 
likelihood ratio test is an adequate statistical test, even 
with small samples. 


(R. E. Stoltz) 
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7.4 (9.1) 


On the accuracy of current approximate variances of treatment differences in 
randomised block designs with missing observations—JIn English 


J. Indian Soc. Agric. Statist. (1958) 10, 83-89 (3 references) 


In this paper the author states the two principal methods 
for computing the standard errors of estimated treat- 
ment differences in randomised blocks with missing 
observations, and makes certain comparisons. 

The two methods stated are, firstly, the multiplication 
of the plot variance by the sum of the reciprocals of the 
actual number of replications of the two treatments 
involved and calculation of the square root: the second 
method makes use of an “‘ effective number of replica- 
tions ’’ instead of the actual number of replications. 

The rule suggested by Yates [Emp. J. Exp. Agric. 
(1933) 1, 129] for getting the effective number of replic- 
ations is the following: subtract from the existing 
number of replications half the number of those replicates 
of the other treatment which are in blocks in which 
the first treatment is not missing. Later, Taylor [Nature 
(1948) 162, 262] suggested that replacing the factor 1/2 
in Yates’ method by 1/(k—1) where k is the number of 
treatments, would produce an improvement. 
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The present paper makes certain numerical com” 
parisons of the values of the standard error given by 
the various rules with the exact values which are known 
in the case of comparisons between two treatments only 
one of which is affected. The author concludes that 
Yates’ rule tends to over-estimate the standard error, 
whereas the errors computed by using the actual number 
of replications are generally smaller than the exact 
values. An average of the values arrived at by these 
two rules provides a close approximation to the exact 
value, closer than values given by Taylor’s rule. 


(S. John) 


7.6 (9.2) 


A comparison of the precision of three experimental designs employing a 


concomitant variable—In English 


Psychometrika (1958) 23, 335-353 (21 references, 6 tables) 


Three techniques are commonly employed to capitalise 
on a concomitant variate and improve the precision of 
treatment comparisons: firstly, stratification of the 
experimental samples and use of a factorial design, 
secondly, analysis of covariance, and thirdly analysis of 
variance of difference scores. The purpose of this paper 
is to compare the effectiveness of these alternatives in 
improving experimental precision, to identify the most 
precise design and the conditions under which its 
advantage holds, and to derive, in the case of the 
factorial approach, recommendations as to the optimal 
numbers of levels. 

Data presented by the author tend to justify the 
conclusion that the less stringent assumptions of the 
factorial design more than compensate for the relatively 
small advantage in precision which may obtain for 
covariance. This is especially true in educational and 
psychological research in which the number of degrees 
of freedom for error are usually quite large and relatively 
little is known about the form of the relationship between 
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criterion and control measures. The more general 
applicability of the factorial design becomes apparent 
from the relationship between the homogeneity of 
regression from treatment population to treatment 
population and the phenomena of interaction between 
treatments and levels in the two-factor design. It is 
shown that heterogeneity of regression is equivalent to 
such an interaction, and the presence of either implies 
the other. 

Heterogeneous regression renders the covariance 
technique, as it is typically applied in educational and 
psychological research, somewhat invalid. If the usual 
covariance model is used, the effects would appear to 
be more serious than those of non-normality and hetero- 
geneity of variance are to an analysis of variance. 
Further research may indicate that such violations of 
assumptions is less serious than heterogeneity of variance 
or non-normality; however, no such conclusion seems 
warranted at this time. 


(R. E. Stoltz) 
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7.3 (9.1) 


The single Latin square design in psychological research—Zn English 
Psychometrika (1958) 23, 369-378 (10 references, 4 tables) 


It has been assumed that interactions must be non- 
existant for results to be adequate in using the Latin 
square designs. Gourlay, using a variance component 
analysis, indicates that for a valid application of the 
Latin square techniques, however, interactions do not 
always have to be zero. Furthermore, contrary to 
McNemar’s assertion, he found that too few significant 
F values might result; a negative F-test bias. Gourlay 
investigated this problem with reference to two main 
types of interactions that occur in psychology. Firstly, 
where each individual or unit receives only one of 
several treatments and is represented by one measure- 
ment in the data. In this case interaction is between 
main effects. Secondly, repeated measurements are 
made on the same individuals or groups. In this case 
earlier measurements may interact with those that 
follow. The author feels that a more general and 
instructive procedure would be to determine the com- 
ponents of variance included with each mean-square 
under four conditions: zero, one, two, and three random 
variates. The first condition corresponds to the fixed 
variate model, the last to the random variate model, 
and the others to the mixed model. He states that this 
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procedure would exhaust all possibilities in the single 
Latin square design and would permit an evaluation of 
the behaviour of each test of significance. 

The author demonstrates the possible defects inherent 
in the single Latin square tests of significance when 
interactions occur under the four conditions mentioned 
above. The data presented show that as the number 
of random variates increases the number of unbiased 
F-tests increases likewise until in the three random 
variates model all tests are free of bias. Paralleling this 
trend is an increasing number of valid F-tests as well. 
In the first two models negative bias occurs most 
frequently but small positive bias is possible when the 
triple interaction is present. In the third model all biases 
are negative. As the number of random variates in- 
creases, more tests are insensitive to deviation from the 
assumptions that interactions must be zero. Less bias 
will occur with interactions present inasmuch as some 
tests will be unbiased; however, a fewer number will be 


valid. Contrary to McNemar’s assertions, and in agree-° 


ment with Gourlay, most of the tests have negative bias: 
the probability of a type 2 error increases with these 
Latin square designs for some tests of significance. 

(R. E. Stoltz) 


7.2 (4.2) 


Maximum likelihood estimation of genetic components of variation—Jn English 
Biometrics (1960) 16, 369-381 (6 references, 4 tables) 


In this paper the methods of least-squares, weighted- 
least-squares, maximum likelihood, and one due to 
Nelder [Heredity (1953) 7, 111-119] are applied to the 
problem of estimating the additive, dominance, and en- 
vironmental effects in a genetic population. The ordinary 
least-squares method is unsatisfactory when the observed 
components of variance are correlated or have unequal 
errors, and for this reason Mather [Biometrical Genetics 
(1949). London: Methuen] considered weighted-least- 
squares estimation using weights obtained from replica- 
tions within the experiment. Assuming normality, the 
maximum likelihood estimates can be obtained from the 
weighted-least-squares by an iterative procedure. 

The similarities and differences of the estimates are 
illustrated from results of a breeding experiment on the 
plant Nicotiana rustica, for which four generations 
were available in 1954 and five in the succeeding year. 
The analyses are carried through on plant height for 
1954, and time of opening of first flower for both years. 
For the maximum likelihood method, five iterations 
are calculated although the chi-square goodness of fit 
stabilises at the second iteration. 
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In the analysis on height, the methods agree that 
heritability is high due to additivity of gene action, but 
that dominance and epistasis are unimportant. How- 
ever, they differ in assessing the significance of environ- 
mental effects. 

For the analysis of flowering time, four models are 
considered ranging from fitting 10 parameters allowing 
for seasonal changes in both genetic and environmental 
components, to a model with 4 components not allowing 
for seasonal changes. It is found that a model having 
varying environmental but constant genetic components 
(6 in all) fits the data almost as well as the most general 
model. A more detailed investigation shows that the 
environmental effects decrease by about one-half over 
the two years, but also suggests that the (significant) 
dominance component has decreased somewhat. These 
conclusions are based on the maximum likelihood 
method; if unweighted-least-squares are used the con- 
clusion concerning dominance is that it is not even 
significant in either year, a contradiction due to the 
incorrect weighting of the latter method. 


(G. A. Watterson) 
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LORD, F. M. (Educational Testing Service, Princeton, N.J.) 


7.6 (7.5) 


Large-sample covariance analysis when the control variable is fallible—In English 
J. Amer. Statist. Ass. (1960) 55, 307-321 (9 references, 2 figures, 1 table) 


When the control variable, or covariate is subject to 
errors of measurement, the usual analysis of covariance 
is no longer valid; its use may lead to equivocal con- 
clusions. In particular, it may point to significant 
differences among adjusted-group-means when none in 
fact exist; it may likewise fail to detect such differences 
when they do exist. 

The purpose of this paper is to outline an appropriate 
procedure for dealing with analysis-of-covariance data 
in the simple case where there are two groups, composed 
of individuals on each of whom we have one measure- 
ment on the criterion variable and two duplicate measure- 
ments on a control variable. The hypothesis to be 
tested is that the two regressions are collinear. 

The test developed in this paper is obtained by a 
large-sample approach. The basic step is the construc- 
tion of an asymptotically unbiased estimator of the 
difference between the adjusted group means. Since 
the estimator is a continuous function of moments with 
continuous first derivatives, its distribution is approxi- 
mately normal for large samples. Once the large-sample 
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standard error of the estimator has been obtained, it is 
possible to make use of familiar procedures for testing 
significance. 

The author provides an example; not only as a guide 
for practical computation, but also to illustrate the fact 
that the method explained can demonstrate a significant 
deviation from the null-hypothesis in data where no 
such deviation would be shown by the usual analysis of 
covariance methods. 


(M. Dorff) 


Notes on the analysis of variance of logarithms of variances—Jn English 
WADC Technical Note 59-82 (1959), Aero. Res. Labs., Wright-Patterson Air Force Base. 


iv+47 pp. (15 references, 13 tables) 


One of the assumptions needed in the analysis of variance 
is equality of within-cell variances. Often the validity 
of this assumption is investigated by the use of Bartlett’s 
modification of the Neyman-Pearson L, test. The 
present report is concerned with the investigation of an 
alternative test, the analysis of variance of logarithms 
of variances, which seems to have some more desirable 
properties. 

The technique consists of dividing the observations 
in each group into sub-groups, applying the logarithmic 
transformation to the sub-group variances, and then 
performing an analysis of variance on the logarithms. 
Let us suppose that we are given random samples of 
equal size J, one from each of K normal populations. 
Let y,; be the jth observation in the kth group, j = 1, 
2,....J; k = 1,2, ..., K, where the y,; are NID(y,, 02). 
We wish to test the hypothesis Ho : of = 03 =... = ox. 
We divide the observations in each group into M sub- 
groups each of size AWA = J, M>1, J>1), then 
compute the logarithm of the variance of the observa- 
tions within each sub-group, and perform a simple 
one-way analysis of variance on the logarithms of 
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variances. The ratio of the mean-square between groups 
to the mean-square within groups is denoted by Fy. 

It is shown that, under the null hypothesis, F, has 
asymptotically, as A increases, an F distribution with 
(K—1) and K(M-—1) degrees of freedom. Since the 
exact distribution of F, for small samples is mathe- 
matical intractable, a study is made of the Type I error 
rates if the percentage points of the F distribution are 
used in testing the significance of F;. For K = 2(1)7, 
9, 11; M = 2(2)8, 12; and A = 2(2)8, 12 it is found 
that Pr[F,=Fo.05] is in the range 0-041 to 0-056 and 
that Pr[F,=Fo.091] is in the range 0-007 to 0-015, so 
that no serious errors result from the use of the per- 
centage points of the F distribution. An example of 
the use of the proposed test is given. 

A test in this form would be expected to be “‘ robust ” 
(relatively insensitive to non-normality), which the 
Bartlett test is not, and would also have the advantage 
that we could break out single-degree-of-freedom con- 
trasts. The question of the power of the test is now 
being investigated, and results will be given in a sub- 
sequent report. 

(H. L. Harter) 
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7.5 (2.2) 


Notes on the use of F in the analysis of variance of attributes data and on the 
transformation of measurements data to attributes data—In English 


ARL Technical Note 60-118 (1960), Aero. Res. Labs., Wright-Patterson Air Force Base. 


iv-+42 pp. (11 references, 14 tables) 


In the usual procedures for the analysis of variance, 
significance tests are based on the assumption that the 
data come from populations having normal distributions. 
The recognition that, in practical situations, this assump- 
tion is at best only approximately fulfilled has led to a 
continuing investigation of how lack of normality 
affects the results of the significance tests. Box [Bio- 
metrika (1953) 40, 318-335] and others have shown 
that the usual F tests on means are “ robust ’’, that is 
that they are insensitive to lack of normality. The 
analysis of variance technique therefore remains a potent 
analytical tool even though a basic assumption inherent 
in it is not strictly satisfied. An applied statistician, 
when faced with analysing data that, obviously are 
not from a normal population, may be able to overlook 
this: fact and still expect to get useful results. Cochran 
[Biometrika (1950) 37, 256-266] has suggested that the 
F test might serve as an approximation even when the 
table consists of 0’s and 1’s, since both the F test and 


2/345 


OLDS, E. G. & LEWIS, J. S. 


5. (Carnegie Institute of Technology, Pittsburg) 
Notes on the use of F in the analysis of variance of attributes data and on the 


the conventional y? test assume normality, and it is 
not obvious which is more sensitive to the assumption. 
Batson [National Convention Transactions (1956) Amerian 
Society for Quality Control, 9-23] has described a 
development of the x? test as applied to the analysis of 
attributes data, and has ascribed to Brandt the formula- 
tion of this method, called “factorial chi-square ”’. 
This report is concerned with the comparison of this 
method to the usual analysis of variance as applied to 
attributes data. In the simple case of n replications at 
each of k levels of a single factor, both test statistics, 
called F* and y*, have as their numerator the mean- 
square between levels. F* has as its denominator the 
mean-square within levels, and F* is treated as an F 
variable with (k—1) and k(n—1) degrees of freedom. 
y* has as its denominator pq, where p is the proportion 
of zeros and q = 1—p? is the proportion of ones in the 
experiment, and y* is treated as a y? degrees of freedom 
variable with (k— 1) degrees of freedom. Asymptotically 


continued 
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continued 


transformation of measurements data to attributes data—I/n English 
ARL Technical Note 60-118 (1960), Aero. Res. Labs., Wright-Patterson Air Force Base. 


iv+42 pp. (11 references, 14 tables) 


(as n> 00) the distributions of F* and x* are identical, 
since both converge in distribution to y7/(kK—1). On 
the basis of a complete enumeration of the possible 
values of F* and y* for several pairs of values of n and 
k (nk small), it is conjectured that F* will always be as 
good as, or better than y* if nominal 95 per cent. points 
are to be utilised. 

The problem of finding the best method of transform- 
ing measurements data from a two-factor experiment to 
attributes data is discussed and exemplified. 


(H. L. Harter) 
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7.5 (2.2) 


Some asymptotic results on transformations in the analysis of variance—Jn English 


ARL Technical Note 60-126 (1960), Aero. Res. Labs., Wright-Patterson Air Force Base. 


(13 references) 


Several transformations have been used in the past to 
stabilise the variance in particular problems of the 
analysis of variance. Most of these studies use only 
approximations. The difficulty, or impossibility, of 
formalising the heuristic argument involved here has 
been demonstrated by Curtiss [Ann. Math. Statist. 
(1943) 14, 107-122], who sought to obtain satisfactory 
asymptotic solutions. Recently, for the logarithmic 
transformation, Severo & Olds [Ann. Math. Statist. 
(1956) 27, 670-686] have shown that, if the exact solution 
is considered, different tests applied to the same data 
for the same hypotheses have very distinct power 
properties while all those tests have identical power 
functions asymptotically. Thus the solutions in the 
study of transformations can be broadly classified into 
two groups. The first involves obtaining certain approxi- 
mate solutions, and the second proving certain asymptotic 
results. The present note belongs to the latter group. 
The problems considered, and the results obtained, 
in this note may be stated as follows. The asymptotic 
properties of the power function obtained by Severo 
& Olds are shown to hold, under similar conditions, 
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also for the square-root transformation (variance 
assumed known), so that this type of result is not the 
property of a very special distribution. Clearly, there 
is no reason to expect the same behaviour of different 
tests for any such (stabilising) transformation in the 
non-asymptotic, or exact case. 

Next, the problem of testing the mean of a log-normal 
distribution when the variance is not known is con- 
sidered. In this case a necessary and sufficient condition 
for the usual normal-theory analysis of variance pro- 
cedures to hold asymptotically, as the mean tends to 
infinity, is that the coefficient of variation (of the log- 
normal variate) should tend to zero. 

Finally certain general results, in the spirit of Curtiss’ 
paper, are presented. These deal with the problems of 
normal convergence of certain random functions. 
Sufficient conditions for the applicability of the normal 
theory are presented for a class of distributions depending 
on a finite set of parameters with one parameter large, 
while the others, if any, are relatively small, or are con- 
fined to a fixed bounded set in the parameter space. The 
last section includes some remarks on related problems. 

(H. L. Harter) 


7.7 (5.6) 


A multiple comparison sign test: treatments versus control—In English 
J. Amer. Statist. Ass. (1959) 54, 767-775 (3 references, 4 tables) 


A number of treatments are to be compared with a 
control. The hypothesis to be tested is that the median 
of each of the other k treatments equals the median of 
the control. The test statistic is the vector (r,, ..., 7;), 
where r; is the number of times the sign of the difference, 
ith treatment effect less the control effect, is negative 
in n trials. (It is assumed that the ith treatment effect 
and the control effect are never equal. A simple rule 
of thumb is suggested for use if they are.) If the r; 
denoted positive signs, the distribution would remain 
unchanged. The author demonstrates the combinatorial 
method for working out the distribution of (r,, ..., 7) 
for the case where the null-hypothesis is true. From this 
k-variate distribution one can then obtain the distribu- 
tion of r= minimum r;. This is tabled for k = 2, 
n = 4(1)10 and for k = 3,n = 4(1)7. 

If the alternative hypothesis is that at least one of the 
other treatment medians is greater than the control 
median, then one rejects the null-hypothesis when r 
is less than or equal to the tabular one-tail critical 
value of r. This is tabled for k = 2(1)9, n = 5(1)20, 
and probability in the tail of 0-05, 0-01. The ith treat- 
ment is declared to have a median significantly higher 
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than the control median, if r; is less than or equal to the 
one-tail tabular r. 

If the alternative hypothesis is that at least one of 
the other treatment medians is unequal to the control 
median, then one should redefine r; to correspond to 
the less frequent sign. One then rejects the null- 
hypothesis when r is less than or equal to the tabular 
two-tail critical value of r. This is tabled for k = 2(1)9, 
n = 6(1)20, and probability in the two tails of 0:05, 0-01. 
The ith treatment is declared to have a median signifi- 
cantly different from the control median, if r; is less than 
or equal to the two-tail tabular r. 

In the last section the author conjectures that the 
distribution of (r,, ..., r,) is asymptotically multivariate 
normal and that (minr—y,)/o, is distributed approxi- 
mately as Dunnett’s ¢ with infinitely many degrees of 
freedom. A comparison with their tables gives 
close agreement. The final approximation used is 
r = 1(n—1—1m’*), (Dunnett’s ¢ with degrees of freedom 
equal to infinity). 


(S. W. Nash) 
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ZYSKIND, G. & KEMPTHORNE, O. (Iowa State University, Ames) 
Treatment errors in comparative experiments—In English 
WADC Technical Note 59-19 (1960), Aero. Res. Labs., Wright-Patterson Air Force Base. iv+78 pp. 


7.5 (9.4) 


(45 references, 6 tables, 3 figures) 


The primary concern of the present report is with the 
investigation of a representative number of alternative 
experimental schemes in which explicit allowance is 
made for the fact that application of intended amounts 
of treatment may be subject to error. These repre- 
sentative experimental schemes are, firstly, the completely 
randomised one-factor experiment; secondly, designs in 
generalised randomised blocks with three alternative 
randomisation schemes for assigning treatment attempts 
to the plots within blocks; and lastly, two-factor 
completely randomised designs with two alternative 
randomisation schemes for assigning treatment attempts 
to experimental units. 

In Chapter Two, the authors demonstrate how the 
question of treatments subject to error may be structurally 
incorporated into a broader formulation of experimental 
design. This more general formulation involves notions 
of the explicit act of randomisation, the use of design 
and sampling random variables, balance (in the occur- 
rence of levels of a treatment in the experimental design), 
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(45 references, 6 tables, 3 figures) 


the & expansions are explicitly written down. The 
design and sample random variables are introduced 
and the statistical model then presented. Finally an 
analysis of variance table showing the sources of variation 
and expected values of their corresponding mean 
squares is shown. 

Consideration of the above-described experimental 
schemes which were investigated allow one to abstract 
the following overall trends in the conclusions. If 
treatment errors are operative then repeated independent 
attempts at realising the preassigned amounts of treat- 
ments are to be desired since without them partial or 
total confounding of treatment amounts and treatment 
errors will result. The presence of treatment errors 
tends to decrease the sensitivity of tests of significance 
for treatment effects. As is also the case with inter- 
actions involving experimental units, the presence of 
interactions involving treatment errors tends not only 
to decrease further the sensitivity of tests but also to 
introduce negative bias both into tests of significance for 
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x expansions of the variances of effects, rightmost 
bracket subscripts, and complete ambivalence. The last 
three are new notions which are defined in the report. 
Using these notions, the authors state, without proof 
a theorem expressing the expected mean-square of a 
partial sample-mean as a linear combination of Sigma 
expansions. For further details of this more general 
formulation, reference is made to Error Structures in 
Experimental Designs [Zyskind (1958) unpublished 
Ph.D. thesis, Iowa State Univ. Library, Ames]. As a 
result of this theorem, when certain sets of structural 
conditions are satisfied, the expected values of mean- 
squares in the analysis of variance tables have simple 
and easily specifiable forms. These conditions have 
been verified to hold in all of the experimental schemes 
discussed in the present report. 

In discussion of these experimental schemes the 
following procedure is generally used. The population 
structure is symbolically represented and the population 
identity is given. Applying their general definitions, 
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7.5 (9.4) 
continued 
iv+78 pp. 


treatments and into the estimation of the treatment 
component of variation; but the presence of these inter- 
actions tends to over-estimate the average variance of 
estimates of treatment differences. Since the magnitudes 
of non-additivities are scale-dependent, the results of 
the present report lend further support to the importance 
of seeking a scale which as nearly as possible achieves 
additivity. 


(Mary D. Lum) 
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8.9 (5.7) 


A note on the efficiency of sequential sampling plans based on gauging—ZJn English 
Bull. Calcutta Statist. Ass. (1960) 9, 117-121 (3 references) 


Sampling plans based on gauging are operationally 
simpler, though less efficient, compared to those ‘based 
on variate values. The author demonstrates that in 
tests of hypotheses specifying the mean of a normal 
population with known variance, against one-sided 
alternatives, the loss of efficiency for sequential plans 
based on gauging is not high. 

Observations are drawn one by one and, at each stage 
the total number of observations having values less 
than G,, the number lying between G, and G, and the 
number above G,, are noted. The observations are 
supposed to come from a normal population with mean 
@ and variance unity. The hypothesis to be tested is 
that 0 = 0 against the alternative that 0 = 0,>0. 
The sequential probability ratio test for testing this 
hypothesis is constructed. 

The average sample number varies with the values 
of G, and G,. By numerical computations the author 
finds that the average sample number can be minimised 
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by setting G, = 0)—0-61 and G, = 0,+0°61. It is 
noted by the author that the optimum gauge setting is 
independent of the value of 6,. 

It is of interest to inquire how efficient this procedure 
is when compared to sampling plans based on variate 
values. Taking the ratio of the average sample numbers 
in the two cases as a measure of relative efficiency, the 
author evaluates the efficiency of the sampling procedure 
under consideration for 0;—6 ) = 0:25 (0-25) 1-0 and 
6,;—@) = 0-1. In each case the efficiency is about 
eighty per cent., though there is a tendency for it to 
decrease as 0, — 0, increases. 


(S. John) 


8.1 (4.9) 
and General Analysis Corp.) 


Design and estimation in two-way stratification—Jn English 
J. Amer. Statist. Ass. (1960) 55, 105-124 (7 tables, 11 references) 


One may encounter populations for which there are 
two effective stratifying criteria, both of which are 
desirable in a sample design. However, the number of 
permitted observations may be less than the number of 
strata formed by the usual double stratification tech- 
nique. The authors present a method which will permit 
estimation in these cases. Both biased and unbiased 
estimators are considered. 

It is shown that, if the stratification effects are additive 
in the analysis-of-variance sense, the method is par- 
ticularly effective. It is stated that if the substrata 
sizes (formed from the two-way classification of the 
population according to the stratifying criteria) are 
proportional to the product of the corresponding one- 
way strata sizes, the possible loss in efficiency compared 
to single stratification is trivial. Even in populations in 
which substrata disproportion is great, one can still 
use the method effectively by employing a method of 
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allocating, with certainty, some of the sample observa- 
tions to the substrata. 


Variances of both biased and unbiased estimators are 
given, along with a method for obtaining essentially 
unbiased estimates of the variances. 


(H. O. Hartley) 
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- and normally distributed error. 


CURNOW, R. N. (Agric. Res. Council, Statist. Unit, Univ. Aberdeen) 


8.6 (8.0) 


Consequences of errors of measurement for selection from certain non-normal 


distributions—In English 


Bull. Int. Statist. Inst. (1960) 37, III 291-308 (5 references). Discussion abstracted No. 2/354 


The mathematical problem can be briefly summarised 
as follows: the true values of a certain characteristic 
are distributed over an infinite population: observed 
values of this characteristic are subject to an additive 
A new population is 
formed by selecting a predetermined fraction of the 
original population possessing the highest observed 
values. The distribution of true values within this 
selected fraction is to be studied. 

This type of selection problem occurs frequently in 
plant breeding work. Assuming that the true values 
are initially normally distributed and given certain 
details of the amount of land available and the cost 
involved, optimal procedures have been proposed for 
the conduct of plant breeding programmes. The 
sensitivity of these optimal procedures to departures 
from normality clearly needs investigating. 

Finney, D. J. has provided a method for studying the 
effects of errors on selection from a general distribution 
specified by its cumulants. See “ 
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The consequence of 


(Report of discussion on a paper by) 


selection for a variate subject to errors of measurement”’ 
[Rev. Int. Statist. Inst. (1956) 27, 22-29]. However, 
more direct methods are available when the distribution 
takes one of a certain number of forms. This paper 
presents the methods involving Hermitian probability 
functions and recurrence relationships that are being 
used to study selection with error from rectangular, 
gamma (or y”) and y-distributions. 


(R. N. Curnow) 


8.6 (8.0) 


Consequences of errors of measurement for selection from certain non-normal 


distributions—Jn English 
Bull. Int. Statist. Inst. (1960) 37, I 59-61 


Those taking part in the discussion included Hemelrijk, 
Neyman and Hamaker. 

Hemelrijk commented that the proportion P was 
introduced in the paper as the proportion of varieties 
selected in the sample; but in the mathematical treat- 
ment the quantity n, related to P by Pr{y>n} = P was 
treated as constant. This meant that P was a proportion 
in the population instead of the sample, thus, selection 
in reality referred to y, taking a value greater than 7, 
and was only asymptotically equivalent to selecting the 
fraction P of the largest value in the sample. For large 
samples there would not be much difference in practice, 
but, from a theoretical point of view he would like to 
known whether the relation between P in the sample 
and 7 was an asymptotic one. The author replied that 
the results were asymptotic in that they depended on 
the assumption that the number of varieties was very 
large. He hoped to be able to investigate the importance 
of the assumption. 

Neyman stated that both Mrs P. C. Tang and Nash 
had investigated exactly the same problems: that is, 
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those dealing with the selection of the number of 
varieties and intensity selection. The author replied 
that he had referred to Mrs Tang’s work but that she 
had dealt with slightly different problems. 


Hamaker, while drawing attention to some industrial’ 


applications of the problem hoped that the author would 
pay some attention to these partially solved applications. 
The author replied that while certain industrial problems 
were probably closely analogous to those discussed the 
work done by O. L. Davies [‘‘ The design of screening 
tests in the pharmaceutical industry ’’, Bull. Int. Statist. 
Inst. (1959) 36, III 226-241: abstracted in this journal 
No. 2/149 and No. 2/151 9.4 and “Some statistical 
aspects of the economics of analytical testing ”, Techno- 
metrics (1959) 1, 49-62] was worth mentioning as being 
probably of rather more direct application to the in- 
dustrial problems. 


(W. R. Buckland) 
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HAMAKER, H. C. (Philips Research Labs., Eindhoven) 


Attribute sampling in operation—In English 


8.8 (-.-) 


Bull. Int. Statist. Inst. (1960) 37, If 265-281 (9 references, 4 tables, 4 figures) 


Discussion abstracted No. 2/356 


This paper brings a critical survey of attribute sampling 
techniques as applied in industry to-day. In setting 
up a sampling table we always attempt a compromise 
between a number of conflicting requirements. The 
choice of the most suitable sampling method depends 
on the features that seem most desirable. Opinions differ 
and this explains the variety in sampling tables that 
have been proposed. Various features of these tables 
are discussed in detail such as the acceptable quality 
level concept, the relation between lot size and sample 
size, the use and utility of tightened and reduced 
inspection, etc. 

It is emphasised that the purpose of sampling in- 
spection is not only to separate bad lots from good, 
but also to control the quality of the lots by making 
use of the information provided by the samples. This 
second purpose is too often disregarded. For purposes of 
control it would be advantageous to use a constant 
sample size independent of the lot size. This is illus- 
trated by an example. 
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HAMAKER, H. C. 
Attribute sampling in operation-—In English 
Bull. Int. Statist. Inst. (1960) 37, I 134-138 


Those taking part in the discussion were Hald, Barnard, 
Riordan, Anson, Madeira and Wallis: the author 
replied. 

Hald took the view that, although MIL-STD 105A 
could be regarded as a convenient catalogue of sampling 
plans with corresponding operating characteristic curves, 
as could the author’s system be regarded as a collection 
of single sampling plans, it was not a good thing to 
throw back to the user the difficult problem mentioned 
in the first section of the paper. This was, however, 
the effect of the author’s definition of a sampling plan 
as being merely the sample size and the acceptance 
number. Instead of a move to greater simplification 
of plans it was probably better in the long run to try 
and make them more realistic: for example, with 
respect to the prior distribution, the costs and also the 
mechanism of statistical feed-back. 

Barnard accepted the economic approach proposed 
by Hald and gave a practical plan for using a feed-back 
system to determine whether too much sampling was 
being carried out. Riordan drew attention to the 
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(Report of discussion on a paper by) 


Finally, the author discusses the question of choosing 
the most desirable features of a universal sampling 
standard. 

It is pointed out that normal, tightened, and reduced 
inspection might well be incorporated in a single table 
provided with sliding scales for both lot size and accept- 
able quality level. A modification of the well-known 
Military Standard 105A in this sense is presented in 
detail. . 


(H. C. Hamaker) 


8.8 (—.-) 


distinction between factory data and data from per- 
formance in service. The latter was both difficult to 
collect and to process. A difficulty of using MIL-STD 
105A was the changing inspection work-loads. 

Anson preferred to regard the problem as consisting 
of two parts (i) the discrimination between good and 
bad lots; which was a single-stage decision problem 
when all the information had to be collected in a first 
sample, and (ii) the control quality of lots; which was 
a multi-stage decision problem. 

The author replied at some length on the difficulty 
of pursuing the economic viewpoint in practice and 
gave an example illustrating the difficulty of dealing 
with data from performance in service. 


(W. R. Buckland) 
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-HAMAKER, H. C. (Phillips Research Labs., Eindhoven) 


Qualitative control of the sample—In French 


8.1 (-.-) 


Rey. Statist. Appl. (1960) 8, 5-40 (24 references, 6 tables, 9 figures) 


In this paper the author first states the fundamental 
principles which are the basis for quality control by 
sampling and then discusses efficiency (operating char- 
acteristic) curves. Graphs are given showing com- 
parisons of several efficiency curves in simple sampling 
plans. 

Quality control as applied to double, multiple and 
sequential plans which allow economy of observations 
is also discussed and some comments made on the 
methods of choosing sampling plans for particular 
situations. 

Acceptable Quality Level is also mentioned as another 
factor which influences the choice. A figure in this 
section illustrates the efficiency curves corresponding to 
the simple sampling plan P = 0-95 for p = 2:5 per 
cent. when the acceptance numbers are 0, 1, 2 and 3. 

It is stated that in the well-known Military Standard 
105A [Sampling procedures and tables for inspection by 
attributes, U. S. Government Printing Office (1950)] the 
probability of acceptance corresponding to the Accept- 
able Quality Level covers the range of numbers as far 
as 85 per cent. Also included in this section are a 
study of the important differences in the main objectives 
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of sampling; economic theories which are related to 
the choice of a sample design and a discussion on the 
distributions of the percentage of defective lots in the 
controlled batch. It is stated that there may be no 
precise information on this distribution: figure four 
illustrates various forms that the distribution can take. 
The following paragraphs contain detailed studies with 
illustrative tables of the relationship between the size 
of a batch and the size of a sample. 
In the latter part of this paper, the special character- 
istics of the principle sampling plans are examined: 
(i) US Military Standard 105A 

(ii) Columbia tables 

(iii) Philips SSS 

(iv) German AWF 

(v) German VG 95083 

(vi) Swedish military standard 

(vii) Dodge & Romig tables. 
Special emphasis is placed on the interest to be found in 
the use of a sample of constant size independent of the 
batch size. The author investigates the possibility for an 
international standard of sampling instead of the many 


plans now in existence. 
(Mile. Turlot) 


8.8 (8.3) 


Comparison of the principal acceptance sampling systems with a new German system—Zn French 
Rey. Statist. Appl. (1959) 7, 5-25 (6 references, 9 tables, 6 figures) 


A new sampling system to be called ASQ-UA6 has been 
proposed by the Sampling Plans Sub-Committee of the 
German Association for Statistical Quality Control. 
The author of this paper has discussed this system and 
also compared it critically with the best known of other 
sampling systems; the systems used for his comparisons 
include the Dodge-Romig system, the Phillips SSS, and 
the systems MIL-STD 105 and AWF. All these are 
compared and discussed at considerable length and the 
paper includes numerous tables and graphs with which 
the author illustrates his statements. 

Firstly the author gives an account of the origin and 
the necessity for the new system and, after stating that 
it has much in common with the MIL-STD 105A, he 
proceeds to describe its advantages. 

After extensive comparisons with other systems he 
also compares it with the UDE 0560 Regeln fur Konden- 
satoren. He states that the presentation varies in each 
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of the systems of sampling but that it is easy to regroup 
the tests according to individual preference. Tables are 
also given comparing the VDE 0560 system with the new 
ASQ-UA6 system. 


(Mlle. Gervaise) 


ek: 
Ey 
% 


“ ee < 4 the 5 
OR? eo 


i 


1 ae ak Pi 
Pe aah) aoe a 
SP Riva dhe IE Ce 


he aaa f ie 
Py.” AY ) “a i ee are re 
UBT hat, Sine Wa a ete ae ee 


MATERN, B. (Forest Research Inst., Stockholm) 


8.0 (10.0) 


Spatial variation. Stochastic models and their application to some problems in 
forest surveys and other sampling investigations—In English 


Meddelanden fran Statens skogsforskningsinstitut 49-5 (1960) 144 pages (133 references) 


The theorey of spatial processes has in recent years 
come to the fore in a number of important applications, 
from macrocosmic patterns such as the distribution of 
galaxies in space, to micro-structures such as the surface 
pattern of manufactured products, e.g. photographic 
film, paper, metal; and between these extremes a great 
many phenomena such as the spatial distribution of 
plants or animals in the field, of trees in a forest, the 
pattern of various rock formations on a geological map, 
the variation of tensile strength in a piece of metal. 
Matérn gives a review and development by way of a 
flexible treatment that starts with generalities about 
processes that are stationary and isotropic (invariant to 
translation and rotation), continues by a detailed 
exploration of specific processes at different levels of 
generality, and makes use of statistical data for a 
discussion of the realism of the various approaches. 
The bulk of the paper is devoted to statistical applica- 
tions, including discussions of the efficiency of random 
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Spatial variation. Stochastic models and their application to some problems in 


and systematic designs in plane sampling, the estimation 
of the sampling error from the data of a systematic 
sample, and specific questions connected with large- 
scale forest surveys. 

Most of the processes considered are constructions 
based on the Poisson process, as when centres y, are 
distributed at random over the plane with constant 
probability density, and the plane is divided into poly- 
hedric cells Y; such that all points x of Y; have y; for 
nearest centre. In developing formulae for covariances 
of integrals, the author studies the distribution of 
distances between points chosen at random in regions of 
different kinds. A formula is given for the average 
surface content (per unit volume) of the boundary of a 
random set. This average is shown to be proportional 
to —c’(0) where c(v) is the (supposedly isotropic) 
covariance-function of the corresponding process. When 
dealing with data read from maps the author treats the 
influence of errors of measurements and rounding-off 


continued 


8.0 (10.0) 
continued 


forest surveys and other sampling investigations—Jn English 
Meddelanden fran Statens skogsforskningsinstitut 49-5 (1960), 144 pages (133 references) 


errors on the observations of a realisation of a process, 
the effect of “‘ local integration ”’ (e.g. of fertility in field), 
and the effect of competition between plants, trees, etc. 

As to systematic sampling with fixed intensity, it is 
shown that points in a triangular network give lowest 
variance among the investigated sampling schemes for 
the classes of isotropic processes considered. The travel 
distance connected with different random and systematic 
sampling schemes is discussed. 


(H. O. Wold) 
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STECK, G. P. & OWEN, D. B. (Sandia Corporation, Albuquerque, New Mexico) 8.8 (8.9) ; 
Percentage points for the distribution of outgoing quality—In English | 
J. Amer. Statist. Ass. (1959) 54, 689-694 (5 references, 2 figures) 


When only a few, perhaps L, lots are going to be pur- 
chased, the average outgoing quality limit will be less 
pertinent for the consumer than the percentage points 
of the distribution of outgoing quality. The authors 
derive formulas for these percentage points as functions 
of incoming quality. 

Three cases are considered: the case where the given 
L lots all have the same known quality; the case where 
the given L lots have different but known quality; and 
the case where the L qualities associated with the L 
given lots constitute a random sample from a distribution 
with known mean and variance. The first case is handled 
exactly; the second and third are handled approximately. 


(H. T. David) 
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BOCK, R. D. (University of Chicago) 


Remarks on the test of significance for the method of paired comparisons—ZIn English 


9.5 (7.3) 


Psychometrika (1958) 23, 323-334 (10 references, 1 table, 2 figures) 


A three-component model for comparative judgment 
which allows for individual differences in preference is 
proposed. An implication of the model is that errors 
in the observed proportions due to sampling indivuals 
in paired comparison experiments are correlated. By 
neglecting this correlation, Mosteller’s test for the method 
of paired comparisons tends to accept falsely the good- 
ness-of-fit of the Case V solution. It is shown that 
bounds may be set for the correlation effect which make 
a valid test possible in some cases and provide useful 
standard errors for the estimated affective values. — 

In 1951, Mostelier suggested the use of the arcsine 
transformation for proportions to test whether the 
variance of those discrepancies is in excess of that 
expected from the binomial sampling variability of the 
observed proportions. When this test is applied to 
data from moderate-sized samples, it persistently shows 
that discrepancies are smaller than those expected from 
sampling variability. That is, the fit of the Case V 
model appears to be too good rather than too poor. 

This paper shows that the anomalous behaviour of 
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Mosteller’s test is the result of assuming the sampling 
errors independent when, in general, they are not. It 
is concluded that under Case V assumptions, sampling 
errors for comparisons involving a common object are 
correlated and on certain assumptions, bounds for this 
correlation may be set. With a minor alteration, 
Mosteller’s test may be recast in the form of an analysis 
of variance and the effect of correlation on the variance 
due to departure from the Case V solution may be 
derived. It is apparent that the binomial sampling 
variance used as the error term by Mosteller is in general 
too large and that the test must frequently fail to detect 
departure from internal consistency in the Case V 
solution. 


(R. E. Stoltz) 


9.0 (9.7) 


The analysis of field experiments according to Mitscherlich—Jn German 
Zeit. Acker Pflanzenbau (1959) 109, 255-290 (11 references, 16 tables, 1 figure) 


The paper contains a comprehensive discussion of the 
adjustment methods due to the so-called Mitscherlich 
designs, which are systematic designs having the plots 
in a row with constant order of the treatments in every 
replication 

In the first part the author assumes that the soil 
fertility is constant for all plots. He computes the 
estimates of the difference between two treatment effects 
and of the error variance, without applying any adjust- 
ments and applying the moving average and Mitscher- 
lich’s adjustment method. As one can expect, all the 
estimates are unbiased and the variance of the estimated 
difference may be greater after adjustment than that of 
the corresponding unadjusted values. Therefore, in 
this simple case of a constant soil fertility, the adjustment 
methods are of no value. 

In the second part the author investigates the influence 
of a variation of soil fertility, which he assumes to be a 
parabola of the ordered plot numbers. Comparing the 
different adjustment methods he obtains the following 
results: firstly, without applying any adjustment, the 
estimates of the difference between two treatments and of 
the error variance obtained by analysis of variance 


methods are biased; the bias being dependent on the 
linear and quadratic term of the soil fertility function and 
on the two treatments compared. Secondly, after 
adjustment by the moving average method the estimates 
are unbiased and can be computed by analysis of 
variance applied to the adjusted values. 

By the adjustment method of Mitscherlich and the 
regression method of Rundfeldt, the biases are decreased 
but do not vanish. However, they depend only on the 
quadratic term of the soil function. The author discusses 
the problem of which systematic row-design will adjust 
the soil fertility in the best way if the latter can be 
expressed by a parabola. He develops his “balanced 
arrangement ”’, in which the treatments are so arranged 
that: (1) their average place within the blocks is constant 
and equals the middle point of the whole design, and (2) 
the sum of squares of the place numbers within the whole 
design is constant for all treatments. These conditions 
can only be fulfilled exactly, if the product n(p+ 1) is odd 
(n =number of blocks, p = number of treatments). 
The author gives some examples of such arrangements. 


(B. Schneider) ; 
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DAMARAJU, R. A. (University of Bombay) 


9.1 (—.-) 


A note on the construction of group divisible designs from hyper-graeco-latin cubes 


of the first order—In English 


Bull. Calcutta Statist. Ass. (1959) 9, 67-70 (7 references) 


In this paper the author gives a method of construction 
of group divisible designs from hyper-graeco-latin 
cubes of the first order. Consider s* treatments arranged 
_ Im_s squares each of side s. Suppose a complete set of 
(s*+s—2) mutually orthogonal latin cubes of side s 
exists. Superimpose all the latin cubes on the arrange- 
ment of treatments in the cube formed by s squares. 
For each latin cube the sets of s* treatments corres- 
ponding to each of the letters give us s blocks with s? 
treatments in each. Thus s(s* +s—2) blocks are obtained. 

Treatments in a particular row of all the squares give 
a block of s* treatments. The s rows give rise to s 
blocks. 

Another s block can be obtained from the s columns 
in a similar way, so that all together, from the letters, 
columns and rows, we generate s?(s+1) blocks with s? 
treatments in each. Let us consider s? treatments 
contained in a square as a group. There are s such 
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EISENHART, C. (Statist. Engineering Lab., Nat. Bur. Stand., Washington) 


Some canons of sound experimentation—In English 


groups. It is easily seen that the s?(s+1) blocks as 
derived above provide a group divisible design with the 
following parameters. 


b=s*(s+1), k=s7, v=s?, 


An =S+1 


r= s(s+1), A, =‘s, 


The group divisible design corresponding to s = 3 is 
constructed to illustrate the procedure. 


(K. R. Shah) 


9.0 (9.7) 


Bull. Int. Statist. Inst. (1960) 37, III 339-350 (9 references, 3 tables). Discussion abstracted No. 2/367 


Randomisation, replication, and balanced arrangements 
are presented as fundamental principles of the planning 
of experiments in general. In the case of comparative 
experiments, one needs to add to the above mentioned 
the use of controls and ‘‘ blind-fold techniques”. These 
principles are elucidated and illustrated in terms of the 
1954 Salk poliomyelitis vaccine field trials. This study 
consisted of two parts, one planned and carried out 
in accordance with the above-mentioned principles, and 
the other without regard to these principles. 

The advantages of balanced arrangements in lab- 
oratory work are illustrated in terms of a titration 
experiment devised by Youden. This experiment 
involves only eight observations, but the structure and 
scheduling of these is such that, from the results obtained, 
it is not only possible to determine the concentration 
of the solution titrated, but also: 


(i) to compare the upper and lower halves of the 
burette to determine the effect of repeated fillings 
of the burette, 

(ii) to compare the two analysts, and 

(iii) to determine whether the analytical procedure 
is subject to a systematic error in the sense of a 
plank 7. 


(C. Eisenhart) 


EISENHART, C. (Report of discussion on a paper by) 
Some canons of sound experimentation—Jn English 
Bull. Int. Statist. Inst. (1960) 37, I 57-59. 


Those taking part in the discussion were Hamaker, Box, 
Martin and Wold. 

Hamaker drew attention to the fact that the list of 
canons of sound experimentation was not complete. 
He mentioned three examples of canons which equally 
contributed to the success of an experiment, among those 
he included the choice and the setting of the levels of 
the factors. 

Box said that while the author stated that the experi- 
ments must be designed in advance and not as they 
proceeded he was sure that he did not mean to imply 
that the individual experiments, which were often steps 
in a larger investigation and so became iterative processes, 
should not occur as such. The object of the statistical 
methods in these situations was to make the experimental 
process converge as certainly and as rapidly as possible. 

Wold stated that while he acknowledged that the 
author had made a significant contribution to the 
classification of the principles of controlled experiment- 
ation, some of these principles could be re-stated in 
formal terms and so further clarified. 
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(Phillips Research Labs., Eindhoven) 


: 9.0 (9.7) 


The author thanked the participants and said that as 
the comments were non-controversial and of a con- 
structive nature he did not feel the need for further 
comment. 


(W. R. Buckland) 


9.1 (—.-) 


Some examples relating to the design of experiments—In French 
Rey. Statist. Appl. (1960) 8, 9-27 (8 tables, 9 figures) 


After saying that the design of industrial experiments 
constitute a field where statistical methods have an 
important part to play, the author proceeds to present 
some practical examples to illustrate this statement. 

The first example, dealing with power-press swaged 
joints, states nine factors which have influence 
on the quality of the assembly and shows how, in 
practice, all nine factors cannot easily be inspected 
simultaneously. In these circumstances it is preferable 
to make a series of simple experiments, so that one can 
develop a sequential system of experiments. There are 
various choices as to where these experiments should 
start and these are next discussed by the author; these 
include choices from the point of view of production, 
cost or materials. The author illustrates and describes 
these experiments. 

The second experiment uses fuse links to prove once 
again how much simpler and quicker it is to use sequential 
methods. Average results for the fusion time as a 
function of electric current as well as the maximum 
and minimum are presented graphically. The third 


example deals with a steel foundry which produces 
cyclindrical steel magnets in very great numbers. The 
production difficulties which arise in this particular 
design and which include the punching of holes and 
control of the accurate height of the cylinder are dis- 
cussed and once again a simple series of experiments is 
drawn up. 

Finally, the author states that he has tried to present 
his examples and conclusions without using an excessive 
amount of statistical terms: he shows that statistical 
techniques provide powerful methods for analysing 
complex situations. 


(Mile. Turlot) 


. 
{ 


! 
: 
f 


2/368 ; 


GRIDGEMAN, N. T. (National Research Council, Ottawa) 


The lady tasting tea, and allied topics—Jn English 


9.0 (9.9) 


J. Amer. Statist. Ass. (1959) 54, 776-783 (7 references, 1 figure) 


The well-known discussion of the principles of experi- 
mentation, illustrated by a taste-testing problem, in 
Fisher’s Design of Experiments, is the basis of this 
expository paper. Reappraisal is conducted in the light 
of the work of Neyman and of Wrighton, which is 
reviewed briefly. The author supports Fisher’s notion 
of a hypothetical population of experiments. -He 
observes that the null-hypothesis of ‘‘ no sensory per- 
ception’ may be tested without constructing a model 
for an alternative hypothesis; that is, at least from 
one point of view it is a matter of indifference which 
of the following is true: 


(i) that all ladies are binomially classifiable according 
to whether or not they can detect the particular 
tea difference; 

(ii) that all ladies may be ranked according to their 
degree of sensory perception. 

It is noted that in many investigations of sensory 
perception involving only a single subject a marginally 
detectable stimulus often gives a frequency of correct 
responses less than 100 per cent. but significantly higher 
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The lady tasting tea, and allied topics—Jn English 


than random expectation. One may therefore postulate 
a neural mechanism which responds with probability p,. 
In the tea-tasting experiment, p, is taken to be the 
probability of recognising a cup of either type. Cups 
which are not recognised may still be assigned correctly 
by chance. Based on this model the author gives the 
probability that the subject scores 2R successes when 
confronted with 2N cups, known to comprise N of 
either kind. For N=4 (double tetrad), Pr(R), &(R), 
and Var R are plotted as functions of p,. 

The author distinguishes between “ sensitivity ’’ (as 
used by Fisher in this context) and “efficiency” (in 
terms of power). Alternative experiments are con- 
sidered: 


(i) duplicate double tetrad; 
(ii) quadruplicate double pair; 
(iii) eight pair comparisons. 
If the lady is allowed one mistake in each, then the 
sensitivities are equal. Pair comparisons are however 
most efficient and double tetrad least efficient, uniformly 
in p,. Since efficiency depends on the unknown p,, the 


continued 


9.0 (9.9) 
continued 


J. Amer. Statist. Ass. (1959) 54, 776-783 (7 references, 1 figure) 


author asserts that operationally speaking, sensitivity 
is specifiable but unimportant, whereas efficiency is 
important but unspecifiable ”’. 

Fisher’s assertion that randomisation is “‘ the physical 
basis of the validity of the test” is discussed. Mis- 
givings are expressed with regard to the case where 
inter-cup differences (India tea and China tea) are 
present, are randomised, but are not balanced between 
the treatment groups. Finally, the author compares 
randomisation at the “‘ design ”’ stage with a randomised 
decision at the “interpretive ’’ stage. The latter intro- 
duced to obtain precisely the desired level of significance. 


(R. J. Buehler) 


KISHEN, K. & SRIVASTAVA, J. N. 


(Dept. Agric., U.P. and Indian Coun. Agric. Res., 


9.2 (9.1) 
New Delhi) 


Mathematical theory of confounding in asymmetrical and symmetrical factorial designs 


—In English 


J. Ind. Soc. Agric. Statist. (1959) 11, 73-110 (10 references, 14 tables) 


The methods of finite geometry given by Bose & Kishen 
[Sankhya (1940) 5, 21-36] for solving the problem of 
confounding in the symmetrical factorial design 5”, 
when s is a prime or power of a prime, are extended to 
the construction of balanced confounded asymmetrical 
designs, by using curvilinear spaces or hyper-surfaces 
and suitably truncating the associated Euclidean geo- 
metry EG(m, s). 

A hyper-surface in EG(m, s) is represented by aj+ 
Ay fy(X%1) +42 fa(X2) +... +4mfn(%m) = 9, Fi) = Lajjx’, 
a; and a;;¢GF(s). Taking f,(x) = x", the authors prove 
three theorems: 

(i) If d is a divisor of s—1 = p"—1, then x* and 
x®-1~4 will give exactly [(s—1)/d]+1 distinct 
values when x varies from dy to a,_. 

(ii) The power-matrix S=(a{)(r, t =a, 1, ...s—1), 
where aj are non-zero elements of GF(s), is non- 
singular and has inverse S~* = (a; ",a$~"~'). 

(iii) Let y = f(x) = ayx+a,x7+...+a,_,x°"'. There 
exist a set of matrices such that as x is varies 
from a, to a,_,, only (k—1) distinct values of 
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y are obtained, so that including x = 0, we have 
k distinct levels. 


The authors use these theorems for getting a general 
method of constructing asymmetrical designs. The 
method is illustrated by constructing the following: 


(i) 3x3 x2 design in blocks of 6 plots, 
(ii) s* xq design in blocks of sq plots, 
(iii) s xq, Xq_ design in blocks of q,q, plots, 
(iv) 4x3x2x2 design in blocks of 12 plots, 
(Vv) 5; X52 xX... XS, designs in block of 5, x 53 x54%... 
plots, 
(vi) qx 27 and qxp? designs, q being any integer, 
and p odd prime power, 
(vii) (pyp2)” designs in blocks of pf? "pz ” plots. 


A connection between balancing in asymmetrical 
factorial experiments and balanced incomplete block 
designs is discussed. Some general methods of analysing 
partially confounded balanced designs are also discussed. 


(C. S. Ramakrishnan) 


9.1 (3.1) 


On the F-test in the intrablock analysis of a balanced incomplete block design—In English 


Sankhya (1960) 22, 279-284 (4 references) 


The paper aims at justifying the use of F-test for testing 
the significance of treatment effects in the intra-block 
analysis of a balanced incomplete block design. 

The author considers the random assignment of 
blocks of treatments of the design to actual experi- 
mental blocks and the further random assignment of 
treatments to plots within each experimental block. 
This randomisation induces a distribution for the ratio 
of sum of squares as given by the statistic (treatment 
SS)/(treatment SS+error SS). The first two moments 
of this distribution are obtained under the null hypo- 
thesis and compared with the first two moments of the 
Beta distribution B(4[v—1], 4[bk —b—v+1)). 

It is shown that for these two distributions, the means 
are always equal and that the variances are practically 
equal if b(k—1)>2; that is, if the number of experi- 
mental units is moderately large and the intra-block 


errors are homogeneous. It is pointed out that the 
situation is remarkably similar to that of randomised 
blocks and that the variance ratio distribution provides 
a reasonably accurate evaluation of the level of signifi- 
cance. 


(K. R. Shah) 


MORILLON, Y. (Research Institute of Iron Smelting, Paris) 


The theory of randomisation—In French 


9.1 (5.1) 


Rev. Statist. Appl. (1960) 8, 29-44 (2 references, 5 tables) 


This work on the theory of randomisation is based on 
previous work by Kempthorne [Design and analysis of 
experiments (1952) New York: Wiley, and ‘“ Random- 
isation theory of experimental inference”, J. Amer. 
Statist. Ass. (1955) 50, 946-967] and that done by 
Kempthorne & Wilk [‘ Fixed, mixed and random 
models ”’, J. Amer. Statist. Ass. (1955) 50, 1144-1167]. 

The author discusses four types of design which have 
as their common characteristic the utilisation of all the 
experimental units subject to the operation of random- 
isation. The study of the distribution of the ratio of 
mean-squares is simplified under these conditions as 
certain terms which occur in the decomposition of the 
sum-of-squares are not random in the null hypothesis. 

After stating that the designs he intends to study fall 
into four groups, 


(i) Completely randomised designs, 
(ii) Randomised complete blocks, 
(iii) Split-plot designs, and 
(iv) Latin squares 
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the author devotes the next section of the paper to a 
definition of additivity. The following section deals 
with the direct test for the null hypothesis that there is 
no difference between treatments and the fourth section 
is devoted to the basic mathematics necessary to the 
analysis of variance. 

Section five develops mathematically the designs 
discussed in the first section and is illustrated with 
several tables in the standard analysis of variance form. 

In conclusion the author shows conditions for a valid 
application of the randomisation test of treatment effects 
in connection with the four types of designs which are 
under discussion. It is also shown that randomisation 
tests are not justified for certain other hypotheses. In 
the case of the Latin square the various problems merit 
further study. 


(Mile. Turlot) 


9.1 (—-) 


On the construction of some partially balanced incomplete block designs with more than three 


associate classes—In English 


Bull. Calcutta Statist. Ass. (1959) 9, 87-92 (5 references) 


The author defines a partially balanced matrix with 
parameters «,, f, and 0,,,(p, q = 1, 2, ..., 8) as a matrix 
whose elements take integral values 1, 2, ..., s such that: 


(i) the integer p occurs a, times in each row and 
f, times in each column, 
(ii) there is a partially balanced association scheme, 
defined on the rows, and 
(iii) in any two rows j and j’ which are ith associates, 
integers p and q occur together in 0;,, columns 
if j and j’ are ith associates. 


Two balanced incomplete block designs N, and N, are 
defined to be associable if: 


(i) a one-one correspondence can be set up between 
the treatments as well as between the blocks of 
the two designs N, and N,, and 

(ii) if new blocks are formed by combining the 
corresponding blocks of the two designs, any 
two treatments i and j, i from N, and j from N,, 
occur together in these new blocks p times if 
i and j correspond and y times if they do not. 


The author shows that, if in a partially balanced 
matrix each integer is replaced by the incidence-matrix 
of one of a set of associable balanced incomplete block 
designs, the resulting matrix is the incidence matrix of a 
partially balanced incomplete block design. 


(K. R. Shah) 
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10.4 (1.5) 


The congestion time-limit distribution for a fully available group of trunks—Jn Russian 
Teor. Veroyat. Primen. (1960) 5, 246-252 (9 references, 2 tables) 


A fully available group of n trunks is considered under 
the assumption that a Poisson stream of calls with 
constant intensity 2 is serviced. The full availability 
group is a loss-system. The holding time is independent 
of the stream of calls and has an exponential distribution 
with a mean holding time equal to unity. 

Let €(t) = {€0(2), €,(0), ..., €,()} be a random vector, 
where ¢,(t) is the life time of the system in its « state 
(a = 0, 1, ..., m) during time interval [0, t]. The second 
moments of the random vector &(t) are determined as 
rational functions of 4. These results make it possible 
to apply integral and local limit theorems for practical 
purposes. 


(G. P. Basharin) 
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BLACKWELL, D. (Univ. of California, Berkeley) 
Infinite codes for memoryless channels—Jn English 


RE 


10.5 (-.-) 


Ann. Math. Statist. (1959) 30, 1242-1244 (1 reference) 


Consider a memoryless channel for which the input 
alphabet A and output alphabet B are both finite. For 
such a channel an infinite code for transmission at rate r 
is defined to consist of 


(i) a sequence of functions f, (n>0) where f, maps 
[rn]-tuples of zeros and ones into A, 
(ii) a non-decreasing sequence of positive integers 
M, satisfying M,,/n—>1 as n> oo, and 
(iii) a sequence of functions g, (n>0) where g, maps 
n-tuples with components in B into M,-tuples 
of zeros and ones. 


Suppose one is required to transmit an infinite sequence 
of zeros and ones in code across a channel, say (x,, x2, 
.s+y Xqo ++), and that the components of the message are 
made available in order, one component every r™! 
units of time. If the transmitter must send one letter 
every unit of time, the mth letter he sends is f,(x,, X2, -., 
Xtmj) Since only the first [rn] letters of the message are 
available to him. On the other hand, at the mth unit of 
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time, the receiver has received n letters, (11, V2, .--» Va) 
say, and estimates the first M, letters (x,, x2, ..., Xm,) 
of the original message to be g,(y,, ¥2, ---» ¥,)- For the 
rate r = 1, the above discussion is greatly simplified, 
and the author states and proves his theorem only for 
this case. 

The theorem of this paper is that if the capacity of 
the channel exceeds 1, then there exists an effective code. 
A code is said to be effective if for every message (x;, X2, 
vss) Xpo ++.) COnsisting of zeros and ones, only a finite 
number of the receiver’s estimates g1(1), ---. Sa(V1. «+> 
y,), ... are incorrect. The author proves the theorem 
by the use of Shannon’s exponential error bounds and 
states that while it would be desirable to extend the 
theorem to finite-state channels, the bounds are not yet 
known for general finite-state channels. 


(R. Pyke) 


BLACKWELL, D., BREIMAN, L., & THOMASIAN, A. J. (University of California, Berkeley) 


The capacity of a class of channels—In English 


10.5 (10.1) 


Ann. Math. Statist. (1959) 30, 1229-1241 (14 references) 


This paper generalises the definition of capacity of a 
discrete memoryless channel to the case of a class of 
such channels I. This capacity C = supg inf,R,(Q), 
where R,(Q) is the usual transmission note defined as 
the information content of the source with input 
probability Q minus the equivocation of one of the 
channels yinI’. Now R20 is an attainable note for the 
class I’ if there exists a sequence of codes (e"%, ¢,, 7) 
for the class with ¢,0. We say (e”*, ¢,. n) is a code for 
each n-extension of each channel in I if it is a one to 
one map from [e"*] sets of words of length n in the 
receiving alphabet into words of length n encoded in 
the transmitting alphabet such that the maximum 
probability of error using the code is of size ¢,,. 

The main result is that the supremum of the set of 
attainable notes for a class of discrete memoryless 
channels is equal to the capacity of that class and the 
method of proof gives an exponential error bound for 
any attainable note R<C. This result is demonstrated 
through a series of steps interesting per se. A basic 
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The capacity of a class of channels—Jn English 


inequality [see Feinstein, Foundations of Information 
Theory (1958) New York: McGraw Hill] gives the 
same bound for a single discrete channel for the maximum 
probability of error that Shannon [Information and 
Control (1957) 1, 6-25] had given for the average 
probability of error. This inequality permits a simple 
proof of the direct half of Shannon’s Theorem, that the 
capacity of a discrete memoryless channel is not greater 
than the supremum of attainable rates. 

The third theorem then gives an exponential bound on 
the error of a code for a single channel which depends 
only on the number of elements in the transmitting and 
receiving alphabets and C—R. This is then generalised 
to the case of a finite number L of channels in a class 
in lemma 3 and theorem 4. The underlying idea being 
to consider the mixture, with probability 1/L, of the 
L channels in the finite class as a channel and to replace 
the error bound ¢ by ¢/L. 

It is then shown that for fixed transmitting and 
receiving alphabets there is a large finite number of 


continued 


10.5 (10.1) 
continued 


Ann. Math. Statist. (1959) 30, 1229-1241 (14 references) 


channels such that any channel on these alphabets is 7 


close, in several senses, to one of this number. Lemma 5 
shows that if a channel has a sequence of codes (e”*, 
&,» ”) with ¢, =e "® for large n with B>O then this 
same sequence of codes can be used for all channels 
in a certain neighbourhood, determined by one of the 
measures of closeness in lemma 4, of the channel. This 
explains the attention to exponential error bounds. 

The direct half of the first theorem is proved by 
approximating the class I of channels by a certain 
finite set of channels from lemma 4, obtaining an 
exponential error bound for this finite set by theorem 4 
and using lemma 5 to show that such a code must work 
also for I’. The converse half is then proved with lemma 
6 by a variation of a result due to Fano. 


(S. C. Saunders) 
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CHERKASOFF, I. D. (Murmansk) 


10.1 (-.-) 


Criteria for co-ordinate homogeneity continuous Markoff processes—In Russian 
Teor. Veroyat. Primen. (1960) 5, 229-237 (10 references) 


The author considers a Markoff process with an in- 
finitesimal operator 
Re ary eae 
A, — A'(t, x) <r + BY (t, x) rae 
ox’ 0x'0x! 
where x = (x’, ..., x") is a point in Riemann space V, 
with a metric g;(¢, x). The necessary and sufficient 
conditions for the existence of transformation 
LE Pied Cane minor cg) 
which transforms the operator A, into the well-known 
operator 
62 
Ox" Ox 


A? = B''(t) 
is given. 
A statistical example is given in which these con- 


ditions are applied for establishing the density of the 
probability f(t, x; t, €) of a certain Markoff process. 


(I. D. Cherkasoff) 
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DAS, S. C. (Ravenshaw College, Cuttack, India) 


10.1 (-.-) 


A note on the structure of a certain stochastic model—In English 


Sankhya (1960) 22, 345-350 (3 references) 


The author considers the following stochastic model. 

There are (m—1) types of successes with the following 

properties: the occurrence of the rth type of success 

in any trial excludes for the next r trials the occurrence 

of any type of success. After this period the probability 

of having a success of the jth type is pa; (j = 1, 2, ..., 
m-1 

m—1) where >) a; =1 and p is the probability of 
j=1 

having any one of the m—1 types of successes. 

He identifies this as a Markoff chain consisting of the 
states Ep, E,, .... E,-; Where Ey is the state where a 
success can occur and £; is the state where a success 
cannot occur for the next m—/j trials. The system moves 
(OMe 10h, Las 2, ... 00 £,,, and from £, ‘to 
Pp eaeiOr pi —tdgv2,.4 5 Ml CE, = Eo). 

The author finds the stationary or the limiting prob- 
abilities for this Markoff chain and evaluates the mean 
and variance of the recurrence times for different states. 
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As an application the author mentions the study of 
accident data in situations where a worker getting 
involved in an accident and hospitalised for a certain 
time depending on the type of accident, is immune to 
further accidents during that period. 


(S. R. S. Varadhan) 
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DANSKIN, J. M. (Rutgers University, New Jersey, U.S.A.) 


The theory of reconnaissance—Jn English 


10.5 (0.0) 


WADC Technical Note 59-409 (1959) Aero. Res. Labs., Wright-Patterson Air Force Base 


iv+31 pp. (1 reference) 


A study is made of the optimum distribution of aerial 
reconnaissance effort against land targets in the presence 
of decoys. The model considered is one in which the 
reconnoitering forces allocate effort among various 
regions, their objective being the location of the targets, 
assuming that the side being reconnoitered is passive. 
The model is described in terms of photographic recon- 
naissance, but the theory is equally applicable to other 
methods of reconnaissance; visual, radar, or infra-red, 
for example, since no specific properties of photographs 
are used. The basic data assumed as given for the model 
are the probabilities p,; that an object of kind i gives a 
photograph of type j. These appear, vary and vanish 
as the level of reconnaissance changes. We shall be 
concerned with situations in which there is no specific 
military attacking system whose action is associated 
with the results of the reconnaissance. 

The results of the reconnaissance are indicated by 
pins stuck in a map, which is then regarded with great 
interest but without specific action by the appropriate 
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(Rutgers University, New Jersey, U.S.A.) 


commander. Associated with each map will be an 
uncertainty. We shall assume that the reconnaissance 
command is searching for information concerning missile 
installations and decoys, both of which are supposed 
to be of a certain specified type. We ‘shall suppose 
that the command has divided the area of interest into 
regions, in each of which there are just three possibilities: 
either there is nothing in the region; or there is a decoy 
in the region; or there is a missile installation in the 
region. Each pin indicates the apparent situation in 
the region by some device equivalent to the giving of 
three numbers, p,, p, and p3, the probabilities that the 
object is nothing, a decoy, or a missile installation. 
We have of course p;+p.+p3 = 1. We suppose that 
a target of kind 7 produces a photograph of type j 
during a reconnaissance with probability p,;;. The 
probability that a photograph of type j is produced is 
then q; = > PiPij: 

The information function of the theory of com- 
munication is chosen as the measure of effectiveness. 


continued 


10.5 (0.0) 
continued 


WADC Technical Note 59-409 (1959) Aero. Res. Labs., Wright-Patterson Air Force Base 


iv+31 pp. (1 reference) 


The information of a reconnaissance is defined as the 
change in the uncertainty of the map resulting from 
that reconnaissance. The expected information pro- 
duced by the reconnaissance is shown to be 


I(x) =e = PiPij log (4i;1 Pi) 


where x denotes the level of reconnaissance and c is a 
constant. It is proved that under certain reasonable 
assumptions the information function is increasing and 
first convex, then later concave. The problem is thus 
reduced to a familiar one in operations research, the 
maximisation of the sum of several convex-concave 
functions subject to linear side conditions. 


(H. L. Harter) 
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DOBRUSHIN, R. L. 10.1 (-.-) 
Properties of sample functions of stationary Gaussian processes—In Russian 
Teor. Veroyat. Primen. (1960) 5, 132-134 (4 references) 


Let ¢,(@), (| t|<0o) be a separable stationary Gaussian 
process with a continuous correlation function. Then 
the following alternative holds: either (i) for almost 
all w, sample functions of process €,(w) are continuous 
functions of ¢ or (ii) there exists a B>0 such that, for 
almost all w, sample functions €,(w) are such as 
? lim ¢,(@)— lim €,() 2B 

t>to t>to 
for any fy. In the second case almost all sample functions 
do not have points of first order discontinuity. 


(R. L. Dobrushin) 
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FINCH, P. D. (London School of Economics) 10.4 (1.5) 


A probability limit theorem with application to a generalisation of queueing 
theory—IJn English 


Acta Math. Acad. Sci. Hung. (1959) 10, 317-325 (3 references) 


The following single-server queueing system is con- 
sidered: the mth customer arrives at the instant t, 
and his service time is s,,n = 1, 2, .... If he finds the | 
server idle, then, before starting service he waits a time 
y,, Otherwise he commences service as soon as the 
(n—1)th customer finishes it. 

Let w, be the waiting time of the nth customer, then 


REE w,+u, if w,+u,>0 
Reaehh ener eater SO 


where u, =S,+1,—T+1, 4=1, 2,.... Let W(x) = 
Pr(w, <x) and W(x) = lim W,(x). It is proved that 
noo 


under certain conditions W(x) exists and satisfies an 
integral equation. In the case where the input process 
is of Poisson-type further discussions are given. 


(G. Bank6vi) 
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FINCH, P. D. (London School of Economics) 10.4 (1.5) 
Balking in the queueing system GI/M/1—IJn English 
Acta Math. Acad. Sci. Hung. (1959) 10, 241-247 (4 references) 


A single-server queueing system is considered with 
balking; that is to say a customer joins the queue con- 
sisting of mcustomers with probability 5,,; the sequence 
{b,,} is non-increasing, bj) = 1, b, = 0 (kK=N+1). 
A random variable n(f) is defined as follows: Nii, 
k =0, 1, ..., N+1 if at time ¢ there are k customers 
present. Let 7, = n(t,—0), where t,, denotes the instant 
of arrival of the nth customer. Let P, = lim Pr(y, = k) 


n> oo 


and Pé = lim Pr(n(t) = k) (k = 0, 1, ..., N+1). 
t-> co 
The author proves that under certain conditions con- 
cerning the input and the service processes, the limiting 


distributions P, and P; exist and are independent of the 
initial state; expressions for P, and P{ are given. 


(G. Bankévi) 
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FINCH, P. D. (London School of Economics) 10.4 (1.5) 
On the distribution of queue size in queueing problems—ZJn English 
Acta Math. Acad. Sci. Hung. (1959) 10, 327-336 ( 8 references) 


An m-server queueing system is discussed by the author i 
in this paper. A random variable y(t) is defined as . 
follows. ntji— ok. kK = 0, 1-2, ....5 if at time t there 
are k customers present. Let ¢, be the instant of arrival, 
s, the service time and w, the waiting time of the nth 
customer. Let 7, = n(t,—0) and €, = n(t,+w,+5,+9). 
The probabilities P,(t), O,(m), Rm), k = 0, 1, 2, ..., 
are defined by the equations 


P,(t)=Pr(n(t)=k), Q,(n) =Pr(n, =k), R,(n)=Pr(¢, =k). 
The author proves five limiting theorems concerning 
n(t), P,(t), Q,(m) and R,(n). 


(G. Bankévi) 


2/386 


as tis fs, 


FIRESCU, D. (University of Bucharest) 


Estimation functions for the fundamental 
chain with a finite set of states—Jn French 


10.2 (4.8) 
probabilities of a multiple stationary Markoff 


Bull. Math. Soc. Sci. Math. Phys. R.P.R. (1958) 50, 401-410 (5 references) 


Let us consider an s-multiple stationary Markoff chain 
with a finite set of states and suppose that the stochastic 
matrix of the fundamental probabilities of the given 
chain belongs to the regular positive case in the sense 
of Fréchet. 

First some theorems of multi-dimensional asymptotic 
normality are proved for the random variables corres- 
ponding to the number of occurrences of a succession 
of states in a sample of given size. Next it is shown that 
the considered estimation functions for the unknown 
fundamental probabilities of the chain are correct as 
well as Gaussian for the case when either all these 
probabilities are unknown or refer exclusively to one 
or more complete rows of their matrix: the respective 
estimation functions are asymptotically efficient. Ana- 
logous properties are also proved for the estimation 
functions for the inverse fundamental probabilities of 
the stationary chain of the family of the given chain. 


(R. Theodorescu) 
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GANI, J. (University of W. Australia, Nedlands) 


10.1 (-—.-) 


Time dependent results for a dam with ordered Poisson Inputs—IJn English 


Nature (1960) 188, 341-342 (6 references) 


The content of a dam at time ¢, when fed by an input 
X(t) and subject to a steady release, may be defined 
Z(t) = Z(t+6t)—dX(t)+(1—n)ot 

where 6t (O<n<1) is the time during ¢, t+ of that the 
dam is empty. The input is a Poisson process (with 
parameter 2) with respect to time but each input is of 
a constant size, «, and a, and always arriving in that 
order. The process Z(t) is Markoffian when the number 
of inputs up to ¢ is known. Using the device of the 
kronecker delta, the author writes down the transition 
distribution function of the dam content and shows 
that it satisfies differential-difference equations sirnilar 
to those used by Tackacs [Acta Math. Acad. Sci. Hung. 
(1955) 6, 101]. 

The Laplace-Stieljes transforms of the differential- 
difference equations lead to integral equations involving 
the probabilities of emptiness and which are not readily 
solved. Explicit results may be obtained however, by 
methods of direct probability. The author draws 
attention to some of his related work (to be published) 
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where he obtains probabilities of first emptiness 
recursively from a modified form of an equation used 
by D. G. Kendall[ J. R. Statist. Soc. B (1957) 19, 207-233]. 
A simpler method of evaluation using truncated generat- 
ing functions is also given [see also Math. Ann. (1958) 
136, 454]. 

Combinatorial arguments are used to establish the 
probabilities of first and second emptiness. Where 
%, = a, the transition distribution function reduces to 
a form given by Gani & Prabhu [J. Math. Mech. (1959) 
8, 658]. The method outlined in this work is equally 


applicable to the case of p>2 ordered inputs, a, ..., &,. 


(W. R. Buckland) 
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GRIFFING, B. 


(Div. of Plant Industry, CSIRO, Canberra) 


10.9 (2.8) 


Theoretical consequences of truncation selection based on the individual phenotype—JIn English 
Aust. J. Biol. Sci. (1960) 13, 307-343 (10 references) 


Theoretical consequences of truncation selection based 
on the individual phenotype are examined for the 
following cases of increasing genetic complexity: 


(i) an arbitrary number of alleles at a single locus, 
(ii) an arbitrary number of alleles at each of two 
linked loci, and 
(iii) a completely general genetic situation. 


The analyses are facilitated by generalising the 
concept of hereditary units to include not only the gene 
but also units of higher levels of organisation. 

Analyses based on the gamete as the basic unit of 
inheritance, but with a gene interpretation, permit a 
detailed examination of the consequences of selection 
and relaxation following selection for the two-locus 
case. It is shown that the immediate response to 
selection may be different from that predicted on the 
basis of gene analysis if an additive x additive type of 
epistasis occurs. However, due to the ‘‘ mutability ”’ 
of these higher order inheritance units, the population 
mean, on relaxation of selection, decays to that predicted 
by the gene analysis approach. 
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HANANIA, Mary I. 


(University of California, Berkeley) 


Finally, by using the individual as the unit of inherit- 
ance it is possible to make certain statements for very 
general genetic situations. These statements are formu- 
lated in terms of covariances among relatives. 


(B. Griffing) 


10.1 (10.3) 


A generalisation of the Bush-Mosteller model with some significance tests—Jn English 
Psychometrika (1959) 24, 53-68 (5 references, 2 tables) 


The authoress proposes a modification of the Bush- 
Mosteller learning model to test the hypothesis that 
reward and non-reward are equally effective in promoting 
learning in a continuous reinforcement situation. The 
hypothesis is reduced to a single parameter, although 
other parameters remain unspecified. 

The Bush-Mosteller model is explored and the modifi- 
cation is statistically developed with the proposed test 
having the property of being asymptotically locally 
most powerful among all tests of the same size and 
asymptotically similar. The test is illustrated in the 
testing of the single and alternate hypotheses. 


(G. L. Hines) 
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HANS, O. (Czechoslovak Academy of Sciences, Prague) 


10.0 (0.4) 


Measurability of extensions of continuous random transforms—In English 
Ann. Math. Statist. (1959) 30, 1152-1157 (8 references) 


In this paper the author gives probabilistic versions of 
the Tietze and Hahn-Banach extension theorems. Let 
(Q, S) be a measurable space. In the case of the first 
theorem mentioned above, it is shown that if V is a 
random transform of QxM into the reals for M, a 
closed subset of a separable metric space Y; and if V 
is continuous for each @ with | Vw, x) |<s(), then 
there is a random transform T of Qx X into the reals 
which agrees with V on Qx WM and is continuous for 
each w with | T(@, x) |<s(). 

In the case of the latter theorem, it is shown that 
if V is a random transform of Q~x M into the reals for 
M a linear manifold of a separable real normed linear 
space X, and if V is a bounded linear functional for 
each w with norm s(q@), then there is a random transform 
T of Qx X into the reals which agrees with V on Nx M, 
and is a bounded linear functional for each w with 
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HASMINSKY, R. Z. (Moscow) 


norm s(w). The second theorem remains true when X¥ 
is a non-separable Hilbert space and M is a subspace 
of X. The proofs of the theorems depend on a con- 
struction argument and thereafter a proof of measur- 
ability. 


(N. D. Ylvisaker) 


10.1 (0.1) 


Ergodic properties of recurrent diffusion processes and stabilisation of the solution to the 


Cauchy problem for parabolic equations—Jn Russian 


Teor. Veroyat. Primen (1960) 5, 196-214 (22 references) 


In this paper the existence of a unique invariant measure 
for Markoff processes satisfying some (in number 9) 
conditions is proved. This result is applied to obtain 
the asymptotic properties of the solution to the Cauchy 


f : 0 
problem for the parabolic equation ce = Lu when 


t>+0. It is established that these properties depend 
on properties of the solution of the external Dirichlet 
problem for the equations Lu = 0 and Lu= —1. The 
sufficient conditions for them, expressed in terms of the 


: : Ou 
behaviour of the coefficients in the equation Lu = re 


are given in the appendix. 


(R. Z. Hasminsky) 


IADRENKO, M. I. (Ukrainiane) 10.0 (—.-) 
Isotropic Gauss random fields of the Markoff type on a sphere—In Russian 
Dopoy. Acad. Nauk Ukrain. SSR (1959) 231-236 (7 references) 


The author describes the random fields, indicated in the 
heading, on a sphere S,, in (m+1)-dimensional space - 
and on a sphere S,, in Gilbert space possessing, in 
addition, some property of ‘“‘ Markoffness”. In the 
case of m = 2 this property reduces to the requirement 
that whatever the curve K dividing S, into two parts 
may be and whatsoever the points P, and P, separated 
by K, the random function €(P,) and &(P,) are in- 
dependent if the values €(P) on K are known. 


(M. I. Iadrenko) 
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KATZ, L. & PROCTOR, C. H. (Michigan State University) 10.1 (0.6) 


The concept of configuration of interpersonal relations in a group as a i 
time-dependent stochastic process—Jn English | 


Psychometrika (1959) 24, 317-327 (3 references, 7 tables) 


To explain changes in the sociometric configuration of 
a group to time, the problem arises of the extent to | 
which such changes may be viewed as the aggregation | 
of part-processes occurring at the level of two-person ‘| 
choice structures. A possible model is a Markoff chain ) 
in which three possible states are mutual choice, one-way 
choice, and indifference; one realisation for each pair 
of choosing individuals in the group. Choice data for 
an eighth-grade classroom are fitted to this model and 
are used to answer questions of constancy of transition : 
probabilities, order of the chain and sex differences. 

Within the limits of the number of cases available no | | 
evidence was found to suggest that the transition 


probabilities varied with time. For the two-month 
gaps there may be second-order dependence of the 
chain, but the chain of four-month gaps did not depart 
significantly from first-order dependence. No departures 
from the first-order Markoff chain model achieves 
statistical significance among the within-sex choices. | 


(R. E. Stoltz) 
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KATZ, M., Jr. & THOMASIAN, A. J. 


(Univ. of Chicago and Univ. of California, Berkeley) 


10.1 (1.6) 


_ An exponential bound for functions of a Markoff chain—Jn English 
Ann. Math. Statist. (1960) 31, 470-474 (9 references) 


Let X,, X2, X3, ... be a finite r-state (r=2) ergodic 
Markoff chain with arbitrary initial distribution and 
rxr stationary transition matrix P with p denoting the 
smallest non-zero element of P. Let p,, p, ..., p, be the 
stationary distribution for P. Assume further that there 
is only one ergodic class of states Ec R = {I, 2, ..., r}. 

Consider ny real-valued LENDS if deanctine on R. 


Let S,= aD IX), p= Y afl 


any mney). The authors use recurrent 
event theory in a manner analogous to Chung [Trans. 
Amer. Math. Soc. (1954) 76, 397-419] and Doblin [Bull. 
Soc. Math. France (1938) 52, 210-220] to prove the 
following theorem concerning the convergence of 
Te Ss (Oop: 

Let m be a positive integer and let e>0. Then 


Pr{|n~1S,—p| Ze for some n=m}<2Ae8”™ 
Pr{S,, 


Pr{S,,<n(u—«) for some n=m}< Ae B”™ 


and M= 


>n(u+e) for some n=m}<Ae 8’ 


where 


A = 8r/p"(1—e7 8”), B= p?"/2°M?r’. 
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KOLMOGOROFF, A. N. & ROZANOFF, Y. A. (Moscow) 


Since the transition matrix P is permitted to have zero 
entries the theorem can be applied equally to multiple 
Markoff chains; that is, to sums of the form 


sn = oe LX, Xxats Xp42) 


The article closes with an example involving a count- 
able state space Markoff chain, which otherwise satisfies 
the hypothesis of the theorem, and a real-valued function 
defined on the positive integers for which there exist 
no constants A and B such that 


Pr{S,2n(ute)}s Ae? 


Bounds for the one-sided inequalities are also given. 
Some possible applications for such bounds can be seen 
in Chernoff ‘“‘A measure of asymptotic efficiency for 
tests of a hypothesis based on the sum of observations ” 
[Ann. Math. Statist. (1952) 23, 493-507]. 


(W. L. Nicholson) 


10.1 (10.6) 


On a strong mixing condition for stationary random Gaussian processes—In Russian 
Teor. Veroyat. Primen. (1960) 5, 222-227 (9 references) 


In this paper the authors consider conditions which 
guaranteed the strong mixing of stationary random 
process €(t). It is proved, for example, that if the 
spectral density f(A) of the process (ft) is continuous 


and positive (parameter ¢ is discrete); or f(A) is positive — 


and uniformly continuous and for vie A 


ne 7 <f/@s = 


where parameter f is continuous, then strong mixing 
is the case. 


(Y. A. Rozanoff) 
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KAWATA, T. (Math. Dept., Tokyo Institute of Technology) 


10.6 (1.6) 


Some convergence theorems for stationary stochastic processes—In English 
Ann. Math. Statist. (1959) 30, 1192-1214 (13 references) 


Let &(1), — 00 <t<o be a continuous stationary process 
of the second order, in the wide sense, with means zero 
and spectral representations 


E(t) = is 


The first integral is a stochastic integral. The second 
is an ordinary integrai, with p(w) the covariance function 
of the &(t) process and F(A) the spectral distribution 
of the &(t) process. If L(s), a real-valued-function of 
bounded variation on each finite interval, has an L, 
(with respect to dF) Fourier transform k(x), that is 


fe [. ed L(s) — k(x) 


A 
then a process y(t) is defined by 


e“dZ(A),  p(u) = IP 


OO) 


eld F(A). 


2 
lim dF(x) =0 


A>—-© 
Boo 


n(t) = Ip é(t—s)dL(s). 


Here the integral is defined first over finite intervals 
and then a mean-square limit taken. The process y(t) 
is stationary with spectral distribution | k(x) |?dF(x). 


If in addition a real valued function f(t) is given such 


that | | f(t—s) | dL(s)< 0 for all real ¢ then set 


—-o 


X(t) = f(t)+ &(0) and define 


[° X(t—s)dL(s) = n()+ 


| aaah 


For a real valued function K(s), the author considers 
convergence theorems for n| X(t—s)K(ns)ds (as 
noo) similar to the classical theorems for the non- 
random case é(t)= 0. For example, Theorem 3.2 in 
the paper states that if f(s)[1+]|s |#]~' and [1+] s |]K(s) 
are in L,(— 0, 00), if f(t+u)—f(t) = O(u) as u-0, and if 
K(s) is bounded and K(s) = o[|s|~?] as s>0o then 


n Ih X(t—s)K(ns)ds— X(t) le K(s)ds j 5) (*) 
ee n 


— 


é 


Superimposing a random element on a Wiener-type 
formula of Bochner, the author obtains his theorem 4.2 
which states that if K(x) is absolutely continuous in 


Di 3o7 continued 
KAWATA T. (Math. Dept., Tokyo Institute of Technology) 10.6 (1.6) 
Some convergence theorems for stationary stochastic processes—In English continued 


Ann. Math. Statist. (1959) 30, 1192-1214 (13 references) 


every finite interval with x? | K(x) | uniformly bounded, 


T 
Thies Bae | | f(t) | dt is uniformly bounded and if 
T 


T 
lin Te * | fe dt 
-T 


T? 2 


exists for some & then 


lim | X(the7 *'aK(at)dt 
G20} 65 ‘ r s 
= lim x(per*at | K(t)dt 
T7?© BE —-T — 02 


the limits existing and being taken in the mean of 
order two. 


The author then considers the behaviour for large 


T of 
iia PN a erst 
| é(te a 
T 


pe ant 


which, if &(t) were not a process but an ordinary trigono- 
metric polynomial, would be the ordinary periodogram. 
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Under additional hypotheses on the &(¢) process which 
insure that it behaves like a Gaussian process, he proves 
that 

T 
lim &|J(T)-J(T’) |?- (1- 7) 2p*(x) = 0 
>To T 


T’ 


if x £0. The term 2p?(x) is replaced by 4p?(0) if © 


x = 0. Here p(x) is the density function (assumed to 
exist and to be bounded and continuous) of F(x). 
There is a similar theorem on the behaviour of co- 
variance [J(T), J(T’)] for large T and 7”. 


(R. Blumenthal) 
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LESLIE, P. H. (Bur. Animal Pop., Depart. Zoological Field Studies, Oxford) 
A note on some approximations to the variance in discrete-time stochastic models for 


biological systems—In English 
Biometrika (1960) 47, 196-197 (2 references) 


This paper is relevant to a description by Leslie [Bio- 
metrika (1958) 45, 16-31] of a discrete-time stochastic 
model for studying the properties of certain biological 
systems by numerical methods wherein certain empirical 
approximations were proposed to the variance of the 
*““number in existance ”’ at time f+1 conditional on the 
state of the system at time f. 

This current paper proposes a more general approxi- 
mation which is at least as good as the original and 
which has the advantage of being applicable to wider 
ranges for the values of the birth and death rates to be 
found in the mathematical description of the stochastic 
model. 

Such approximations save considerable machine time 
in the performance of, say, a numerical experiment with 
two interacting species since the necessity of including a 
logarithmic sub-routine in the programme is avoided. 


(J. G. Saw) 
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LEUND, C. S. (Moscow) 
On conditional Markoff processes—I/n Russian 


Teor. Veroyat. Primen. (1960) 5, 227-228 (4 references) 


The author considers a pair of random processes X,, 
Y, which jointly form a Markoff process Z, The 
conditional distribution of the process y,, for the condi- 
tion of known realisation of process X, during some 
time interval, is examined. Dynkin has proposed that 
if X, is a Markoff process, the conditional distribution 
of Y, will satisfy a functional equation similar to the 
known Kolmogoroff-Chapman equation. 

The author has proved this proposition, but details 
of the proof are omitted in this paper. 


(C. S. Leund) 
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LUM, Mary D. 


(Wright Air Develop. Center, Dayton, Ohio) 


10.5 (1.0) 


A comparison of an error-correcting eight-unit with other teleprinter codes for 


binary transmission and ternary reception—In English 


WADC Technical Note 58-231 (1959), Aero. Res. Labs., Wright-Patterson Air Force Base 


Vii + 36 pp. (2 references, 4 tables, 6 figures) 


Teleprinter information is transmitted electrically, using 
a sequence of “ marks ”’ and “spaces”. In true binary 
transmission, a mark represents a pulse sent on one 
frequency f,, while a space represents a pulse sent on 
another frequency f,. With this type of transmission, 
the decoder can make three possible decisions (ternary 
reception): firstly, accept pulse on /,, secondly, accept 
pulse on f,, or finally, accept neither the pulse on f, 
nor that on f,. This report is concerned with ternary 
reception of binary transmission. Binary reception can 
be considered a special case of this more general situation. 

As used in this report “‘ error-detection ”’ means the 
ability to detect at least one bit-error in a received 
character, “‘error-correction’’ means the ability to 
detect and correct at least one bit-error in a received 
character. 

A study is made of an eight-unit “four marks plus 
four spaces’’ binary code in comparison with the 
following six binary codes: a seven-unit “‘ four marks 
plus three spaces’’ code with shifts; a six-unit code; 
a five-unit code with shifts; a seven-unit parity check 
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McGREGOR, J. R. (University of Birmingham) 


An approximate test for serial correlation in polynomial regression—In English 


code; a six-unit parity check code with shifts; and an 
The common | 


eight-unit double parity check code. 
feature of these codes is that they can be made compatible 


with the existing 5-unit Baudot code since they all | 


contain at least 62 characters. The purpose is to deter- 


mine whether the proposed eight-unit code is optimum | 


(with respect to minimal susceptibility to character errors) 
for various bit-error rates. 


Formulas are derived for the probability of character | 
errors for the different codes as functions of the pro- | 


portion of bit-errors. 
probabilities versus proportion of bit-errors for the 


Tables and graphs of these 


range of 10 per cent. missing bits to 1 per cent. missing 


bits are given for two cases: firstly, missing bits only; 
and secondly, ten per cent. as many wrong bits as missing 
bits. An Appendix gives a list of the type of characters 
which are error-corrected or error-detected by each 
individual code. 

The results indicate that the eight-unit double parity 
check code is definitely the best of the seven codes 
studied. 


(Mary D. Lum) 


10.7 (2.7) 


Biometrika (1960) 47, 111-119 (10 references, 1 table) 


When the error terms in the regression model are serially 
correlated the usual least-square procedure fails.. Durbin 
& Watson [Biometrika (1950) 37, 409-428 and (1951) 
38, 159-178] developed bounds to the significance points 
of a test criterion d, defined as the ratio of the sum of 
squares of successive differences of residuals to the sum 
of squares of residuals from a fitted regression. 

In this paper, McGregor obtains a Pearson Type I 
approximation to the distribution of d and observes 
that the relative error between the true and approximate 
density functions is of order 1 /n* where n is the sample 
size. In addition for the case of regression on orthogonal 
polynomials, it is shown in an appendix that the exact 
distribution of d differs from that of d, (the upper bound 
given by Durbin & Watson) by order 1/n. 

A table is given comparing 1 per cent. and 5 per cent, 
points of d, with those of the approximate (Type I) 


distribution of d for various values of the sample size 
and orders of the polynomial regression. 

The proof of the results rests on the properties of a 
certain determinant and on an argument due to Daniels 
[Biometrika (1956) 43, 169-185]. 


(J. G. Saw) 


MALECOT, G. 10.1 (-.-) 


The problem of Brownian movement and its generalisations—In French 
Publ. Inst. Statist. Paris (1959) 7, 186-195 


In this paper the author describes two generalisations 
of the “‘ random walk ”’ process. 
The first consists of the supposition that the possible 
jumps, instead of being of constant amplitude, possess 
algebraic values X, with respective probabilities 74. 
The second consists of the admission that the times | 
of these jumps are random variables. | 


(Mile. Gervaise) | 


2/403 


MARYANOVICH, T. P. 10.4 (2.5) 


Generalisation of Erlang’s formulae in the case when the lines may be broken and 
renewed—IJn Russian 


Ukrain. Mat. Z. (1960) 12, 279-286 (3 references) 


A fully available group of n trunks is considered under 
the assumption that a Poisson stream of calls with 
constant intensity is serviced. The author considers a | 
loss-system: the lines may be broken in the moments ¢; i] 
and renewed in the time y; Pr{y<x} = G(x), Pr{t;.,—-— 
t,<x} = H(x), Pr{¢<x} = F(x), ¢ is the service time. 

Two following conditions are satisfied: 


ie [1—G(x)]dx<oo, ih [1— F(x)][1—A(x)]dx<oo. 


The stationary distribution p;;, where i is the number of 
busy lines, j the number of broken lines, and it+jSn 
is given. 


(B. Gnedenko) 


MASUYAMA, M. (Tokyo University) 10.1 (10.6) 
A response pattern analysis—In English 
Rep. Statist. Appl. Res. (JUSE) (1959) 6, 82-96 (10 references, 4 tables) 


The author presents three methods of singling out the 
true response pattern to one arrival of the meteoro- 
logical stressor with or without noise. In fact, he gives 
the methods to solve the single response pattern f(t) 
from the functional equation F(t) = f(t—a)+f(t—b) 
under the condition that F(t), a and b (a<b) are known 
and f(t) is one valued bounded function of finite range. 
He shows that the temporal or spatial n-method which 
is widely used in this field of research to find out the 
response pattern to the meteorological stressors is 
generally biased. In sections 3, 4, 5, 6, 7, he considers 
three methods to solve the simplest case in which only 
two signals arrive at t = aand t= b(>a). Insection 8, 
he extends the above idea to the two-dimensional case 
F(t, ty) = f(t, — ay, ty— a2) + ft, — by, t, —b9). 

Furthermore, in section 9, he treats the case which has 
many signals in one dimension. (See also No. 2/406). 


(S. Sakino) 
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MASUYAMA, M. (Tokyo University) 10.1 (0.4) 
Further studies on response pattern analysis—In English 
Rep. Statist. Appl. Res. (JUSE) (1959) 6, 97-99 (5 references) 


In this paper, the author tries to solve by an iterative 
method the difference equation with constant coefficients 


Co V(t) +c, Y(t—1) +c, Y(t—2)+...+G-; Yt-—k+1) 

+ Y(t—k) = G(t) 
where c,Y(t—a) is a one-valued bounded response 
pattern of finite range to the single signal or stressor 
of intensity c, at t = a and G(t) is the observed response, 


and where the Laplace transform of Y(t) is assumed to 
exist. (See also No. 2/405). 


(S. Sakino) 


NEWELL, C. F. (Brown University, Providence, R.I.) 
Queues for a fixed-cycle traffic light—In English 


10.4 (10.1) 


Ann. Math. Statist. (1960) 31, 589-597 (8 references, 1 figure) 


Analytic expressions are derived for the average queue 
length at a fixed-cycle traffic light, under equilibrium 
conditions. 

The model is a discrete one, in which arrivals and 
departures of cars at the traffic light occur only on a 
set of equally spaced time-points. Each cycle of the 
periodic light consists of r consecutive “red” time- 
points followed by g “green” points. At each point 
one car may arrive (probability «) or cars not arrive 
(probability 1—«); independently from point to point. 
Departures occur at green points only, one departure 
for each time-point so long as there is a car present; 
with no queue, a car may arrive and depart at the same 
green point. 

Queue lengths at corresponding time-points of 
successive cycles form a stationary Markoff chain. If 
q,, denotes the queue length just before the first red 
point of the xth cycle, and wu, the total number of 
arrivals during the xth cycle, the q, satisfy the recursion 
relation 

'G@x+1 = Max 1d as O}. 

This relation is used in Section 3 to find approxima- 
tions to &(q,) where the arrival rate is low. Starting 
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PEARCE, S. C. (East Malling Res. Sta., Maidstone, England) 


A method of studying manner of growth—ZIn English 
Biometrics (1960) 16, 1-6 (1 table, 1 figure) 


Questions about the growth rates of individuals in a 
growing organism can sometimes be answered by 
noting the relationship between the standard error of 
log(size) and the mean increment in log(size) of the 
organism. If x,(t) denotes the size of the ith individual 
at time ¢ then the following three standard errors are 
especially useful in exploring the laws of growth: 


(i) o,, that of log(size) at time f, that is to say of 
Li, = log x(0), 

(ii) o,9, that of the increment in log(size) between 
time 0 and time ¢, that is of L;,—Lio, 

(iii) o,, that of log(size) at time ¢ after adjustment 
by covariance on the corresponding values at 


; o 
time 0, that is of Li,— Po, — Lio» 
or) 
where fo, is the correlation between L;, and Lio. 
The author shows that only two of these are in- 


dependent. 
A study of graphs of the three standard errors when 
plotted against the mean difference between Lj, and 
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with an arbitrary distribution for q,, the sequence of 
distributions for q,, q3, etc., converges to the equilibrium 
distribution if the average number of arrivals per cycle 
is less than the maximum departure rate, that is, if 
a(r-+g)<g. The rate of convergence is very high if 
g—a(r+g)>[a(l—a«)(r+g)]?. An adequate approxima- 
tion to the distribution of g, may be obtained by finding 
the distribution of q, for some small n. Approximations 
for &(q,) are found in the limit r and g>0o, with r/g 
fixed, for the case u>1, where 


bh ig-orte)iirel@+2))-* 

In Section 4 use is made of generating functions to 
handle the case w<1. The generating function is first 
shown to be the quotient of two polynomials of known 
degrees, and these polynomials are in turn shown to 
be determined by the roots of z7—(1—a+az)"*% = 0.* 
For the case 4 approximately equal to unity, a discussion 
is given in Section 5. The paper ends with a comparison 
of results with some previously obtained by Webster 
*“ Traffic signal settings’ [Road Research Tech. Paper 
No. 39 (1958) Her Majesty’s Stationery Office, London]. 


(R. Hooke) 


10.9 (6.0) 


Lo sometimes suggests laws of growth. For example 


(i) if o, =o, from some time ¢ onward it implies 
that the initial size ceases to affect the growth 
rates, 

(ii) if the curve for o;9 continues to rise it means that 
the growth rates for different individuals vary 
and, 

(iii) if the curves for o, and o; are not straight it 
implies that the growth rates are not constant., 


These points are illustrated by two examples, the 
second of which deals with the growth rates of apple 
trees: graphs are given to illustrate the trials. 


(R. L. Chaddha) 
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RUNNENBURG, J. T. 


(Mathematical Centre, Amsterdam) 


10.4 (10.1) 


Probabilistic interpretation of some formulae in queueing theory—Jn English 
Bull. Int. Statist. Inst. (1960) 37, III 405-414 (13 references) 


It is customary to eee the variable z of a generating 


function A(z) = a Ppz" as “having no significance ”’ 


[see Feller, Ttcdiction to probability theory and its 
applications (1950) chapter 11 New York: Wiley], even 
if the quantities p, generated are probabilities. An out- 
standing feature of Van Dantzig’s ‘‘ Method of collective 
marks ” is that it recognises the possibility of inventing 
a “ significance’ for the variable z. Instead of using 
A(z) as a receptacle for probabilities which is devoid 
of intuitive meaning with respect to the problem in 
hand, one may try to interpret A(z) itself as a prob- 
ability in a slightly more complicated problem than the 
original; in this z denotes a probability, which may 
have any value between 0 and 1. 

It is shown by the paper, that the foregoing outline 
provides a method which has distinct advantages over 
other approaches in certain waiting-time problems. 
Three examples are dealt with; all originating from 
the one-server queueing theory with first-come, first- 
served queue discipline, in which the arrival intervals 
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SACKS, J. (Columbia University, N.Y.) 
Ergodicity of queues in series—Jn English 


are exponentially distributed, while the service times 
have an arbitrary distribution 
The examples are as follows: 


(i) Example A treats a formula due to Takacs [Acta 
Math. Acad. Sci. Hung. (1955) 6, 101-129], con- 
cerning the virtual waiting-time at an arbitrary 
moment. 

Example B gives a new derivation of the well- 
known formula for the Laplace-transform of the 
waiting-time distribution in the stationary situ- 
ation, due to Pollaczek [see page 156 of Kendall, 
D. G., J. R. Statist. Soc. B (1951) 13, 151-185]. 
Example C supplies intuitive content to a de- 
composition of the Pollaczek formula [see first 
footnote on page 208 of Kendall, D. G., J. R. 
Statist. Soc. B (1957) 19, 207-233]. 


(it) 


(iii) 


(J. T. Runnenburg) 


10.4 (10.1) 


Ann. Math. Statist. (1960) 31, 579-588 (5 references) 


Two queues are said to be in series if the output of the 
first queue is the input of the second queue. If W, is 
the waiting time of the nth individual in the first queue, 
and W* his waiting time in the second queue, then the 
queueing system is said to be ergodic if the joint dis- 
tribution of (W,,, W;) converges to a probability distribu- 
tion as n becomes large. 

In this paper, the conditions under which queues in 
series will be ergodic are derived; firstly for the case of 
two servers and then in the general case of s servers. 
The first significant result in the development is: 

PriZ.=1|Z,— 0} F(t) asin. 
where F is a two-dimensional distribution function whose 
variation over two-dimensional space may be less than 
one, that is, F may not be a probability distribution 
function. In the above, Z, = (W,, W, a). 

The second part of the argument is to show that 
under appropriate conditions, Z, is bounded in prob- 
ability as n becomes large, so that F will be a probability 
distribution. This is contained in a theorem: if 


&(g,)>max[é(R,), &(p;)] then F is a probability 
distribution. Here R, is the service time of individual n 
at the first server; p, is his service time at the second 
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server; and g,.; = T,+41—T,, Where t, is the time at 
which individual n enters the system. It is assumed 
that each of the three sequences {R,; nZ21}, {g,; n22}, 
{P,; n= 1} is a sequence of independent and identically 
distributed random variables, and that the three 
sequences are mutually independent. 

The necessity of the condition of the above theorem 
is shown, when first moments are assumed to exist, by 
the following: 


(i) If &(Ri)2 &(g2), then F(t) = 0, 

(ii) If &(p1)2 &(g2), then F(t) = 0. 

The question of ergodicity in the case of s servers is _ | 
handled in much the same way. Here let R; be the 


service time of individual n in the server o, o = 1, 2, 
was Let Ro, , =%.4,—1,, and let p, = &(R4), o = 0, 


1, 2, ..., s. Then the analog to the theorem is: if max 
1SoKs 
He<Mo then F is a probability distribution: where F 


is the s-dimensional analog of F in the first displayed 
expression. Finally, if max ,2WUpo then Fis identically 


1<o<s 


Zero. 


(R. M. Durstine) 
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SCOROCHOD, A. V. 


10.1 (.-) 


On differentiability of measures which correspond to Markoff processes—Jn Russian 
Teor. Veroyat. Primen. (1960) 5, 45-53 (5 references) 


The conditions for the differentiability of measures that 
correspond to Markoff processes are investigated and 
the density of one measure with respect to the other is 
calculated. 


(A. V. Scorochod) 
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SHISKIN, J. (Bureau of the Census, Washington, D.C.) 


10.0 (11.5) 


Decomposition of economic time series by electronic computers—In German 
Math. Tech. Wirtschaft. (1960) 3, 112-118 (5 references, 4 tables, 3 figures) 


The article is a German translation of Shiskin’s paper 
in Science (Vol. 128, pages 1539-1546). That paper, in 
turn, is a summary of several articles by Shiskin & 
Eisenpress on this topic: the use of computers in the 
decomposition of economic time series [see, for example, 
J. Amer. Statist. Ass. (1957) 52, 415-449]. These articles 
reflect the pioneer work carried out by Shiskin in this 
field during the last few years. 

This article in Science, as translated in Math. Tech. 
Wirtschaft first discusses the different kinds of economic 
fluctuations; the author next goes into the possibilities 
of applying computers in the decomposition of economic 
time series, into the trend (including cyclical components) 
and into seasonal and irregular components. 

In this paper, no details as to the method itself are 
given. The details of the rather elaborate and smooth 
mechanical method of decomposition as developed by 
Shiskin are set down in the other papers mentioned 
above. Broad emphasis, however, is given to tests of 


the “‘ goodness of fit”? of the method. Of the 10,000 
series so far decomposed by the method, 150 leading 
series were carefully tested in several ways to allow a 
more accurate evaluation of this mechanical method of © 
decomposition. 


(G. Bruckmann) 


Ki, icra A paar aN PLY Ri 
NET A oe aa 


i eS Shee Sle 


id asi Ky 


ait 
7 

4 \ 
si 

r 

~ ue 

i 
1 & 

, 


SIEGEL, S. (Pennsylvania State University) 


10.1 (2.5) 


Theoretical models of choice and strategy behaviour: stable behaviour in the 


two-choice uncertain outcome situation—ZJn English 


Psychometrika (1959) 24, 303-316 (22 references, 1 figure) 


The author points out that in the classical experimental 
situation, as first used by Humphreys, in studying two- 
choice uncertain outcome situations two widely used 
theoretical models yield somewhat different predictions. 
The Estes statistical learning model yields a prediction 
that subjects will learn to match their response ratios to 
the actual probabilities of occurrence of the events. The 
von Neumann and Morgenstern game-theoretic model 
would predict that a person will learn to maximise the 
expected frequency of correct predictions. To do this 
he will predict the more frequent event on all trials. 
The author purposes a resolution to the apparent 
paradox, a resolution which incorporates the utility- 
theory approach. More specifically, the author is 
suggesting that one should consider the utility of the 
person making a correct prediction of a given event. 
Three models are presented. Model III is of par- 
ticular interest in that it incorporates the Shannon in- 
formation measure and leads to somewhat different 
predictions than are yielded by either Models I or II, 
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SMITH, C. (Animal Breeding Res. Org., Edinburgh) 
Efficiency of animal testing schemes—Jn English 
Biometrics (1960) 16, 408-415 (3 references) 


Testing schemes or procedures are used in theimprovement 
of livestock. The amount of improvement obtained froma 
testing procedure depends on (a) optimum design with re- 
spect to predicted genetic gainand(b) optimum application 
to the whole population. Robertson [Biometrics (1957) 13, 
442-450] derived an expression for optimum design such 
that superiority of chosen families would be maximised. 

The author has considered three cases in which 
individuals, rather than families, are selected, and has 
derived prediction formulas for genetic gain when 

(i) only non-tested members of selected families are 
used for breeding, 
(ii) all members of selected families are used, 
(iii) only the tested animals are available. 


Variables in these formulas include the proportion p 
of groups selected, the genetic relationship r between 
members of the same group, the observed intra-class 
correlation ¢ for the character considered, the heritability 
ratio h?, the genetic standard deviation og, and the 
ratio K of number of animals tested to number of groups 
chosen. Expected genetic gains, using a wide series of 
testing situations for pigs, were calculated and are 
presented graphically. 
2/414 


An experiment is presented designed to test the hypo- 
thesis that asymptotic probability of a person’s predicting 
the occurrence of the more frequent event in a two- 
choice uncertain outcome situation is a function of the 
level of reinforcement present in the situation. It is 
hypothesised that the probability of predicting the more 
frequent event will tend toward unity as the rewards 
(positive utility) and costs (negative utility) of correct 
and incorrect predictions are increased. The data 
presented confirms the prediction. 


(R. E. Stoltz) 


10.9 (6.1) 


Genetic gain showed a nearly linear relationship to 
the logarithm of the fraction to be selected from the 
groups tested. Increases in predicted genetic gain were 
associated with decreases in ft, also with increases in 


size of family. Full-sib families were preferable if the 


fraction to be selected was small; half-sib families if the 
fraction was large. Selection of tested animals alone 
was advantageous if the fraction to be selected was small. 
The use of all members of tested families gave greater 
genetic gains than the use of only non-tested members. 

Selection based on family mean plus individual merit 
should add, in practice, very little to genetic improve- 
ment from selection on family mean alone. If selection 
is to be done simultaneously for two or more traits, 
it is appropriate to use a selection index. 

The individuals selected through the foregoing pro- 
cedures may be used directly to improve the whole 
population, or else used indirectly through a nucleus 
of breeding stock. The nucleus plan, though at first 
inferior, should ultimately lead to more genetic im- 
provement in the whole population than would the use 
of selected individuals directly on the whole population. 


(N. R. Thompson) 
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STRATONOVICH, R. L. (Moscow) 10,0 Ga) 
Conditional Markoff processes—Jn Russian 
Teor. Veroyat. Primen (1960) 5, 172-195 (4 references) 


In this paper the author gives relationships between the 
probabilities of conditional Markoff chains for neigh- 
bouring tests. The conditional probabilities at the end 
of the observation interval, the final probabilities, are 
satisfied with first-type equations corresponding to an : 
increase in the observation interval. The second-type | 
equations for the conditional probabilities within the 
observation interval are written in terms of the final 
probabilities. 

The following special cases are considered: Gaussian 
noise with independent values which becomes a delta- 
correlational process when the moments of time are 
compacted, and a continuous Markoff process. 

The related problem of the time-sign reversal of 
ordinary, a priori, Markoff processes is treated. 


(R. L. Stratonovich) 
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VINCZE, I. (Math. Inst., Hungarian Academy of Sciences, Budapest) 10.5 (6.4) 
On an interpretation of a concept of information-theory—Jn Hungarian 
Mat. Lapok (1959) 10, 255-266 (10 references) 


The J-divergence of a random variable €, with respect 
to a probability distribution ¢, is the expression 


1,8) = iE foo) log [fo /#(a) Jax 


— oe 


where f(x) and ¢(x) are the density functions of the 
variate € and the distribution ¢ respectively. In this 
paper the author gives the following interpretation of this 
notion: the probability distribution ¢ is interpreted as the | 
distribution of our interest and /4(€) as the information : 
contained in the value of € with respect to the distribu- 
tion of our interest. This interpretation is made plausible 
by starting from Shannon’s formula for the discrete 
case and passing to the limit. 

In case € and ¢ are normal with parameters po, 
dp and pl, dp respectively I4(€) = (u—Mo)7/(206). 


(K. Sarkadi) 
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VOGEL, W. (Universitadt Tiibingen, Germany) 


A random-walk problem and its application in the theory of sequential 


design—In German 
Arch. Math. (1960) 11, 310-320 (9 references) 

The following random-walk problem is treated: begin- 
ning at the time nm = 0 in the point k = 0, the random 
variable may jump at the times n = 1, 2, 3, ..., one 
unit in a positive or a negative direction; that is from 
k to k+1 or to k—1. The probability of a positive 
jump is p if k is positive, 1—q if k is negative, and 
r =(p+1—q)/2 if k =0. The corresponding values 
for a negative jump are: 1—p, q, and s = (¢+1—p)/2. 
The case p # 1 —q is of especial interest. 

If U, is the position k of the random variable at 
time n, the author discusses the asymptotic behaviour 
of the following expression: 


S\(P. ) =) PHU, <0)+4P1(U; = 0), 


that is, the probability that the position of the random 
variable will be negative at some of the times 1, ..., 7, 
and half the probability that it will be zero at some of 
these times. He proves several propositions on the 
asymptotic behaviour of this expression. 

These results can be applied to the following sequential 
design: given two possibilities to perform an experiment, 
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YUSHKEVIC, A. A. (Moscow) 


the outcome of which takes the values 1 (success) and 0 
only. The two possibilities differ in the probabilities 
of success which may be denoted by p and q respectively. 
A total of n experiments are to be performed according 
to the following sequential plan: the first possibility 


10.1 (5.7) 


must be used, if in the former experiments the number | 


of successes with the first possibility plus the number of 
failures with the second one is greater than the number 
of failures with the first plus successes with the second 
possibility. If however the second number exceeds the 
first one, the second possibility must be used; and if 
the two numbers are equal, either of the two possibilities 
must be used with probability one-half. If &(V;) is the 
expectation of the ith experiment, the risk function, by 
which the value of this sequential plan may be judged, 
can be written as R(p, q) = nmax(p, q)— XV). 
One can easily show, that this risk function is linearly 
connected with the expression S,(p, q) introduced in 
connection with the random walk discussed above, so 
that the results on asymptotic properties of S, can be 
applied to this risk function. 
(B. Schneider) 


10.1 (-.-) 


On the definition of a strong Markoff process—Jn Russian 
Teor. Veroyat. Primen. (1960) 5, 237-243 (3 references) 


The equivalence of the two definitions for strong Markoff 


Processes given in Dynkin’s book [The basis of the 
theory of Markoff processes (1959), Moscow] is proved. 


(A. A. Yushkevic) 
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ASKOVITZ, S. I. (Univ. Pennsylvania, Philadelphia and Albert Einstein Medical Center) 


11.3 (6.1) 


Graphic methods based upon properties of advancing centroids—In English 
J. Amer. Statist. Ass. (1959) 54, 668-673 (8 references, 6 figures) 


In an earlier paper [J. Amer. Statist. Ass. (1957) 52, 
13-17] the author presented a graphical method for 
obtaining the best fitting, i.e. least squares straight 
line to n points having equally spaced abscissae. In 
the present paper he gives: (i) two alternative solutions 
to the same problem; (ii) a construction for three 
unequally spaced points; (iii) a method which uses the 
best fitting straight line to n consecutive equally spaced 
points to obtain the best line given one additional 
point; (iv) a construction which locates the mean of a 
histogram. Proofs are given, or at least sketched. 


(i) The first alternative consists in repeated addition 
of new points as described under (iii) below. 
The second consists in starting from the second 
point from the left (rather than the first, as 
described in the earlier paper) and obtaining at 
the (n—2)th step the weighted centroid of the n 
points, with weights 0, 1, ..., 7—1; rather than 
1, 2, ...,n. Repeating with left and right inter- 
changed gives another weighted centroid, which, 
with the first, determines the least squares line. 
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BANERJEE, K. S. (State Statistical Bureau, Calcutta) 


(ii) 


(iii) 


(iv) 


Given unequally spaced points A, B, C, let Ao, 
Bo, Co be their projections on the x axis. Con- 
Struct A,, B,, C, alternately above and below 
A, B, C, so that AA, = ByCo, BB, = CyApo, and 
CC, = A By. Then the two intersection points 
of the two polygonal lines ABC and A,B,C, 
determine the least squares line. 

Let P and Q be the trisection points of the given 
best fitting line terminated at the first and nth 
abscissae, and let F be the additional point. 
Locate R on QF two-thirds of a spacing unit, 
horizontally to the right of Q. Then P and 
R determine the new best fitting line. 

The histogram is replaced by points centred on 
the tops of the bars. The mean is located by a 
method which involves the two weighted centroids 
described under (i). 


(R. J. Buehler) 


11.8 (-.-) 


A generalisation of Stuvel’s index number formulae—In English 


Econometrica (1960) 27, 676-678 


The purpose of this note is to generalise the derivation 
of the indexes given by Stuvel in his recent article, 
“New index number formulae’”’ [Econometrica (1957) 
25, 123-131]. 

After discussing Stuvel’s procedure for deriving the 
formulae, the author presents his generalisation which 
satisfies both the time-reversal and the factor-reversal 
tests. The basic identities, of which Stuvel took the 
arithmetic average, are now averaged with weights 
which incorporate an alternating function of price and 
quantity. 


(W. R. Buckland) 
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| 
FITZPATRICK, P. J. (Catholic Univ. of America, Washington, D.C.) 119: C3)s) 
_ Leading British statisticians of the nineteenth century—Jn English 


J. Amer. Statist. Ass. (1960) 55, 38-70 (59 references) 


The author reviews the lives and statistical contributions | 
of eleven outstanding British statisticians of the nine- : 
teenth century, namely, Playfair, Porter, Babbage, 
Farr, Guy, Newmarch, Jevons, Rawson, Galton, Giffin, | 
and Edgeworth. He states that perhaps the most | 
important contributions are those of Galton, who 
introduced the methods of correlation and regression 
analysis, the latter arising out of his studies on the 
relationship between the stature of an adult person and 
the average stature of the parents. 

In the conclusion, the author summarises and com- 
pares the training and professional activities of the men, 
the offices which they held, and the honours conferred 
upon them, together with their contributions to statistics. 


(S. P. H. Mandel) 
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GOLENKO, D. I. & SMIRYAGIN, V. P. (Computing Centre, Acad. Sci. USSR, Moscow) 11.7 (11.6) 
A source of random numbers which are equi-distributed in [0, 1]—Jn Russian 
Publ. Math. Inst. Hung. Acad. Sci. (1960) 5, 241-253 (2 references, 11 figures) 


In this paper the authors deal with an electronic device 
yielding random binary numbers; it makes use of 
occurrences of noise voltage which are greater or less 
than some fixed voltage level. The background from 
which these random numbers arise is of fully stochastic 
nature; thus, they might have some advantage over the 
pseudo-random numbers generated by any of the usual 
methods. The randomness of the obtained random 
numbers has been checked by several statistical tests. 
The results were satisfying; thus the device can be 
applied, for instance, in solving problems by the Monte 
Carlo method. 


(P. Medgyessy) 
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KENDALL, M. G. (London School of Economics) 


11.9 (-S, 


Studies in the history of probability and statistics. X. Where should the history of 


statistics begin ?—Jn English 
Biometrika (1960) 47, 447-448 


The author is concerned to point out that the isolated 
enumerations and inventories of the first sixteen hundred 
years of this present era do not constitute “‘ statistics ”’ 
in the sense that the word is now used. Statistics as 
we understand it only really began when the collectors 
of data began to reason about it and to seek to draw 
inferences which were not the immediate result of the 
straightforward tabulation of the data. The random 
element has been studied from earliest times. The 
collection of data such as that of the Book of Numbers 
or the Domesday Book is also a process of some 
antiquity, but as the author remarks, the marriage of 
the random element to the collection of data had to 
wait until well on into the nineteenth century. 

Kendall dates the real beginning of “ statistics ’? with 
the publication by John Graunt of Observations on 
Bills of Mortality in 1660 .v. 


(Florence N. David) 
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LOCKS, M. O. 


(University of California, Los Angeles) 


11.5 (7.1) | 


Automatic programming for automatic computers—In English 
J. Amer. Statist. Ass. (1959) 54, 744-754 (6 references, 2 tables, 2 charts) 


A description of the various stages of work and their 
attendant difficulties and complications when compiling 
programmes for automatic computers is the main concern 
of the author in the first part of this paper. Having 
described the advantages of automatic programmes, 
these include simplified methodology and the time 
saved by the repeated use of the same programme, he 
describes in detail an automatic programming system 
and enumerates the three basic types: namely the inter- 
pretative, assembly and compiling systems; these are 
briefly described and their merits compared. The 
availability of these three types for different computers 
is discussed: the question of excessively long running 
time and poor adaptability to many types of problems 
is then dealt with briefly and one method of solving 
this problem is given. An example is given of an auto- 
matically programmed analysis of variance F computa- 
tion using the Math-Matic compiler system for the 
Univac I and I. 


Hypothetical costs are estimated for computing the | 


Snedecor F in three different ways. These are as follows: 


(1) Desk calculator 


(ii) Computer with hand coding but without auto- | 


matic programming 


(iii) Computer with automatic coding and employing | 


a compiler. 


Having stated that, for simplicity, linear cost functions | 


will be assumed the remainder of this part of the paper 
gives the various costs and compares them: two tables 
and a chart are used to illustrate this section. 


The final part of this paper deals with the analysis of | 
the results and there is a summary which includes | 


prognostications as to the future of automatic pro- 
gramming and states that anticipated future develop- 
ments include automatic operating and debugging. 


(A. Booker) 
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MACKINNON, W. J. (University of Arizona, Tuscon) 


11.1 (2) | 


Compact table of twelve probability levels of the symmetric binomial cumulative distribution 


for sample sizes to 1000—In English 


J. Amer. Statist. Ass. (1959) 54, 164-172 (15 references, 1 table) 


A table of critical values for tests of the symmetric 
(p = 4) binomial cumulative distribution is presented. 
Compactness of the table obtained is accomplished by 
the use of an unusual table argument, d, defined as 
follows: let s and r be the frequencies of “‘ successes ” 
and of “failures”? in n trials such that s<4n and 
r23n; then the sum, r+s, is m and the difference, 
r—s, is d. The table gives the critical values of s at 
twelve two-tail probability levels (0-001, 0-01, 0-02, 
0-05, 0-10, 0-20, 0-30, 0:50, 0-70, 0-80, 0-90, and 0-95) 
for sample sizes to 1000. The entire table requires only 
those values of d from 1 through 105, thus the use of d 
rather than n as the argument reduces the size of the 
table considerably. 

Since the table has limits cn both sample sizes and 
probability levels, the author discusses other techniques 
for making tests involving the symmetric binomial 
distribution. In addition to giving references to some 
available tables, he gives the maximum error of the 
two-tail probability associated with both the normal 
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PERKAL, J. (University of Wroclaw) 


On the analysis of a set of characteristics—Jn English 


Zast. Mat. (1960) 5, 35-45 


After a short critical review of the advantages and 
disadvantages of known statistical methods used in 
analysis of factors (Hotelling, Spearman, Thurstone) 
the author proposes two new methods [Perkal, J., ““ Two 
new methods of analysing a collection of attributes ”’ 
Bull. Acad. Sci., Wil (1959) 7, 63-66: abstracted in this 
journal No. 1/324, 11.0] which can be summarised as 
follows: 

let x;; be the value of the ith attribute of the jth 
morviduales= 01, 2..,.., Kk, f =1, 2732.7, and: let the 
attributes be standardised over the set of individuals to 
the effect of having zero mean and unit variance. 

If all the attributes are non-negatively correlated, 


then m; = 2 )) x;; is the value of a common factor for 
toe 


the jth individual, while the residuals or indices w;; = 
x;;—m, contain the remaining part of information. In 
the contrary case a procedure for splitting the collection 
of attributes into groups of non-negatively correlated 
attributes is defined, which is based on the matrix of 
correlation coefficients between the attributes. The 
summarising factors are then computed for each group 
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approximation and the Camp-Paulson approximation 
when p = 3. It is noted that when n=5, the normal 


approximation has a smaller maximum error than the | 


Camp-Paulson approximation. 


The justification of the validity of a table using d | 


as the argument depends upon the monotonicity of the 


cumulative probability, with respect to s, for fixed d. | 


A proof of this is given as well as an explanation of 
the basis for the rules which appear in the instructions 
for the use of the table. 
obtain the table entries are given. 


(D. D. Grosvenor) 


Finally, the methods used to | 


11.0 (9.2) 


of attributes separately. A k-dimensional vector with | 
all its components equal to m, in the first case and with | 


equal components for each group of attributes in the | 
second case is assumed as the representation of the | 
common factor; and this procedure of computing it is | 


called by the author the “ single-vector ’’ method. 
A ‘‘multiple-vector’’ method is also defined. It 
consists in successive application of the above described 


single-vector method to the ever smaller groups of | 


indices, which, after suitable standardisation are regarded 
as new collections of attributes. The procedure stops 


when there are no more groups of non-negatively | 


correlated indices. 

In the way described above the factors are com- 
paratively easily and uniquely determined and may be 
immediately interpreted in natural sciences, although, 
unlike in Spearman’s and Thurstone’s methods, they 
must not be mutually orthogonal and uncorrelated with 
the residuals which appear when the computations are 
continued. 


(S. Zubrzycki) 
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